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1. INTRODUCTION 


The MINIVER computer program (ref. 1) Is an automated method of using the 
simplistic approach to aerodynamic heating. Prime benefits from using 
MINIVER computer code are the flexibility and economy resulting from the 
use of a point solution technique. This document describes the changes, 
modifications, and Inclusions which have been adapted to the current version 
of the MINIVER program. Extensive modifications were made to various sub- 
routines, and a complete new plot package added. This plot package is the 
Johnson Space Center DISSPLA Graphics System currently driven under an 
1110 EXEC 8 configuration. This document provides user instructions on 
executing the MINIVER program, gives a detailed description of the changes 
and modifications made, and provides an extensive description of the plot 
package. 
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2. DISCUSSION 


The version of MINIVER which was modified used the cumbersome method of 
allowing data to be I put by use of a data editor subroutine recognizing 
location numbers, rather than conventional Input techniques. As a result, 
any Increase In size of an Input array resulted In failure to execute the 
program. For this reason, the Input data method was changed and will be 
handled by FORTRAN NAMELIST. 

The old subroutines INPUT and INPUTA were completely removed, eliminating 
the need for storing all of the Input data within an array called W with 
fixed 1500-computer word locations. The user Is now free to change array 
sizes of the Input parameters, /1th little or no Impact in the editing and 
execution of the program. 

The main subroutine (H800) was extensively modified to handle storing of 
generated data In an 1110 UNIVAC-secured file, and a subroutine called 
STORED was written to handle this task. The contents of this secured file 
can be retrieved later for plotting purposes. 

The printed output was handled by the main program H800; this method was 
also changed. A subroutine called NEWOUT was written to allow the main pro- 
gram to be almost output free. 

The subroutine WRINP, which also handles the printing of Input data, was 
rewritten to accommodate the changes made with the Inclusion of the new 
Input technique. 

Another subroutine changed was SETMUP. This routine sets up thick-skin 
parameters and writes out the Input. 

Changes in subroutines TRANS and OPTMYZ were minor and are Indicated through 
comments cards within the program listing. 
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Once calculations for a body point have been completed, the user has the 
option of plotting the generated data or storing It for later processing. 
A routine to perform this task was written, and a complete description Is 
given In section 10 of this manual. The main program makes a call to the 
routine which drives the plot package. 


3. INPUT DESCRIPTION 


MINIVER Input Is basically the same as before modification, but the method 
of getting the data In the program has been changed. The old method used 
the ENCODE, DECODE technique; the new version uses FORTRAN NAMELIST. The 
computer name for NAMELIST Is DATALO. 

Table I shows the Input parameters by group, with units and symbols. Two 
new columns have been added: program name (the name used In the program) 

and array size (the number of computer words currently assigned to that 
particular parameter). The array size Is assumed to be a single value when 
none Is shown. 

Table II shows the new parameters Introduced Into the modified version, 
not classified under groups or symbols. 
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parameters by group 






ABLE I. - Continued 
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ENT1 Number of FF $ table entries 






TABLE I. - Continued 
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TABLE 1. - Continued 



HFAC q/q factor 








TABLE I. - Continued 
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TABLE I. - Concluded 



TSINKZ Inner wall sink temperature table 



TABLE II. - NEW INPUT PARAMETERS 


Program name 

Description 

Array size 

ICASE 

Case number 


LN6PLT 

Flag used for plotting. >C long plots 

<0 short plots 


IHCOPY 

Flag used for plotting. >0 no plot generation 

<0 make plots 


IPLTWE 

Flag used for plotting while In Demand 
mode. >0 plot while on Demand 
<0 no plot 


IFLGOT 

! 

Flag used for printed output while on Demand 
mode. >0 print output data 
<0 no print 


BDYPNT 

Body point ID. (Alpha input, 24 characters) 

4 

TRANME 

Trajectory ID. (Alpha input, 36 characters) 

6 

IBPNUM 

Body point number 


MAXTME 

Max time used per run ( in seconds) 
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4. PROGRAM EXECUTION 


There are several ways to execute MINIVER; the most convenient Is In demand 
mode. Upon sign-on, type on device: 

GSTART ES35-N0651 6*L0REN . RUNMIN 

Contents of RUNMIN, Including the overlay, are shown In figure 1. Notice 
the statements: 

0ASG.T 9. 

If the user wishes to plot the data generated during the run at a later date, he 
must save the contents of unit 9 at end of run, as follows: 

0COPY.I 9, FILENAME. ELEMENTNAME 


Close attention should be paid to the @ADD statements, as these files con- 
tain the Input data. The user must generate his own Input files. It Is 
obvious that the user will not get plots on the device with a @START command; 
therefore, If plots are desired, the IPLTWE (I Want Plots While Executing) 
must be set to 1. Plots will be In microfilm. 

Another way to start execution is through the @ADD statement 
@ADD ES35-N06516*L0REN. RUNMIN 

Print output will appear on the screen if the IFLGOT flag Is set to 1. The 
user may want to skip this procedure; if so, set IFLGOT to 0. 

A user familiar with demand terminal operation should be able to choose the 
better method. 
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1 *6RUN,R/R 170LOX, E/3207, ES32-L78771, 20, 50 NODECK LORENZO 

2 1 0FREE TPF0. 

3S0FREE 9. 

4*6ASG,T 9.,F/|/TRK/500 




33* IN DUNSTN 

34 * SEC PHEXPN1S, (FLOU1 ) 

35* IN PNEXPN 

36* SEC DOUNID1S, CFLOU1 ) 

37* IN DOUNID 


*&fr,3**SEC PHIDit# CFLOlli ) 
39*111 PH ID 

- 4G*SEG CRSFLU1*, CHAIN 
41* IN CRSFLU 
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Figure 1. - Concluded 
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6. OUTPUT LISTING 


The listing corresponds to the Input shown In section 5. 
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7. PLOT PACKAGE DESCRIPTION 


MINIVER plots are generated using the DISSPLA utility routines as currently 
operational on the UNIVAC 1110 EXEC 8 operating system at the Johnson Space 
Center. All DISSPLA routines were not used to accomplish this task. 

Section 7.1 shows the subroutines used and gives a short explanation of the 
usage of each. For more detailed information, see the DISSPLA manuals 
(refs. 2 and 3). Several subroutines were written to drive these routines, 
PU3TL0 (the main driver), HNDYLO, DRAWLO, DSMMLO, and PPL0T1. Description and 
usage of these subroutines is given in section 7.2. A detailed flowchart is 
shown in section 7.3. Output plots are shown in section 7.4. 

7.1 USAGE OF DISSPLA ROUTINES 

AXSPLT - To obtain rounded axis scaling parameters for units/inch 
(i.e., GRAPH type) axis 

CALL AXSPLT( AMIN, AMAX.ORIG, STEP ,AXIS) 

AMIN - Least data value 

AMAX - Greatest data value 

ORIG - Returns rounded axis origin 

STEP - Returns rounded step size 

AXIS - Returns rounded minimum axis 

BASALF - To obtain desired alphabet 
CALL BASALF(' STANDARD') 

CURVE - To draw a curve 

CALL CURVE{ XARAY ,YARAY ,NPNTS , IMARK) 

XARAY - Array containing X-values 
YARAY - Array containing Y-values 
NPNTS - Number of points to be plotted 


IMARK - Frequency of marker symbol 

> 0 Points connected with symbols 
= 0 Points connected. No symbols 
< 0 Points not connected 

DISSPLA (automatically) allows up to 14 different symbols 

ENDPL - To end a plot and create a new physical page 
CALL ENDPL (I PLOT) 

I PLOT is the plot number 
> 0 Summary on printer and plot 
= 0 No summary on either 
< 0 Summary on plot only 

ERASE - Erase screen contents 
CALL ERASE 

ERTRAN - To fetch current date and time 
CALL ERTRAN (9, DATE .TIME) 

The number 9 is required by the system. The other two arguments are 
self-explanatory. 

GRACE - To set margin around the subplot area (beyond which curves will be 
scissored) to an arbitrary value. Default is 0.5 inch 

CALL GRACE(GRACEM) 

GRACEM - Width of grace margin around subplot area in inches 

GRAPH - To set up linear axis specified in units/inch 
CALL GRAPH (XORIG.XSTEP.YORIG.YSTEP) 
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XORIG - Value of X at the axis origin 
XSTEP - X step size In units/inch 
YORIG - Value of Y at the axis origin 
YSTEP - Y step size in units/ inch 

GRID - To draw a grid in the subplotted area 
CALL GRID( I XGRID , IYGRID) 

IXGRID - Number of grid lines per X-axis step 
IYGRID - Number of grid lines per Y-axis step 

HEIGHT - To set character height. Default is 0.14 inch 
CALL HEIGHT(HITE) 

HITE - Character height in inches 

INTNO - To plot an integer {in inches) from the physical origin 
CALL INTN0{ INUM, XPOS ,YP0S) 

INUM - Integer to be plotted as a string of digits 
XPOS - X-coordinate in inches 
YPOS - Y-coordinate in inches 

IQWAIT - An interrupt to allow user viewing of display 
CALL IOWAIT ( IARG) 

IARG - Number of seconds of image to remain on screen 

LEGEND - Identifies the curves on a plot by their markers as provided by 
DISSPLA. The text for the legend must be supplied in an array 
which has been packed with the routine called LINES. The sequence 
of the lines should correspond to the order of use of the markers. 

CALL LEGEND ( IPKRAY .NLINES, XPOS, YPOS) 
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IPKRAY - Name of array containing packed lines 
NUNES - Total number of curves to be identified 

XPOS - X-distance from physical origin to lower left corner of legend, 
in Inches 

YPOS - Y-distance from physical origin to lower left corner of legend, 
in inches 

LINES - To pack a line of text 
CALL LINES (LSTRNG, IPKRAY, I LINE) 

LSTRNG - Character string terminated by $ 

IPKRAY - Pack array to receive the line 
INLINE - Sequence number of line being packed 

MESSAG - To plot a message string, in inches from physical origin 
CALL MESSAG (LMESS , I MESS , XPOS , YPOS ) 

LMESS - Characters to be written 
IMESS - Number of characters in LMESS 

XPOS - X-distance from physical origin to start of message, in inches 
YPOS - Y-distance from physical origin to start of message, in inches 

MIXALF - Refer to OISSPLA manuals (refs. 2 and 3) 

NOCHEK - Suppress listing of points out of range. The default option is for 
the point out of range listed on the printer. 

CALL NOCHEK 

TEKEGM - DISSPLA interface with 1110 UNI VAC 
CALL TEKEGM(480) 

System expects the number 480 
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XINTAX - Integer numbering on X-axis 
CALL XINTAX 

YINTAX - Integer numbering of Y-axis 
CALL YINTAX 

YAXAN6 - Angle labels on Y-axis 
CALL YAXANG(ANGLE) 

ANGLE - Angle from horizontal, in degrees 

7.2 MINIVER PLOT ROUTINES USAGE 

PLOTLO - The driver routine for the plot package 

CALL PLOTLO (TIME .NHFLAG ,ARI DEF.ATRE ,TZ ,ZZ , VZ .ALFAOT .DELTAT , ITHICK ,LNGPLT , 
TIN.MAXTME, DEVICE) 

Argument in the call statement as defined in section 3. 

Subroutines required: DRAWLO.DSMMLO.PPLOTl .HNDYLO 
Libraries required : LOCALIB.DISSPLA.PLOTIO 

DRAWLO - A utility routine which collects DISSPLA routines that are called 
several times. The routine was written to avoid numerous calls to 
the same procedures. 

CALL DRAWLO ( XO , XD , XL , YO , YD , YL , XARRAY , YARRAY .NPOINT , I FLAG , I PASS ) 


XO 

- X-origin in inches 

XD 

- X delta increments in inches 

XL 

- X-axis range in inches 

YO 

- Y origin in inches 

YD 

- Y delta increments in inches 

YL 

- Y-axis range in inches 

XARRAY 

- X array to be plotted 
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YARRAY - Y array to be plotted 

NPOINT - Number of X-Y pairs 

IFLAY - To Indicate closing of a given frame 

I FLAG a 1 One plot per frame 

IFLAG » 2 Two plots per frame 

I FLAG * 3 First plot of a set of two 

IPASS - Not used 

Subroutines required: None 

Libraries required : L0CALIB,DISSPLA,PL0T10 

DSMMLO - Find minimum and maximum of an array 
CALL DSMMLO (NP, A, YMIN.YMAX) 

NP - Number of points 
A - Array name 
YMIN - Minimum value 
YMAX - Maximum value 
Subroutines required: None 
Libraries required : None 

HNDYLO - A utility routine to eliminate the calling of the same routines in 
the plot driver 

CALL HNDYLO ( XARRAY , YARRAY ,NP,I LABL .IPASS) 

XARRAY - X array to be plotted 

YARRAY - Y array to be plotted 

NP - Number of x.y pairs 

ILABL - Y axis label (36 characters maximum) 

IPASS - Not used 
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Subroutines required: PPL0T1 ,DSMML0,DRAWL0 
Libraries required : LOCALIB.DISSPLA.PLOTIO 

PPL0T1 - To place a label on the y-axls 
CALL PPL0T1 ( ILABL , ICASE ) 

ILABL - Label to be used (36 characters maximum) 
ICASE - Case number 
Subroutines required: None 
Libraries required : L0CALIB,DISSPLA,PL0T10 
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7.3 FLOWCHARTS 



ZJOl 


SUBRWTINE PLOTU>(TIMR,NHrLAti,ANir>W,ATRE,Tt, 
tt,VZ,Ai/AOT,DELTAT, ITIIICK, UI , ;J'LT,TIH,MA*TKK t 
MTICB) 


JI 


orrnn norm, tranne, device 


377 


□ 


FAJUMKTDl J1 
i_ 


17, J2 500, J3 32 


it; 


TKX rOLLOUIM AM ATS AM DIMENSIONED J2| 

TIME1, 00X1, A ELI, QC1, QC1T,TK»71, TV1, 

Hi, cpi, ami, alpiiai, aiozi, hpacei, 

Tl». 9T(a),TT(a),VZ(5O),TZ95O),ALPAOT(50) 
W( 50) , D«LTAT( 50) ,tDrm(4) ,TKA10U(6) , 

nup(aoo).uaB.( ao.t). ipak(ao) 


L 


ALL PLOT LAB 113 CONTAIN SIX COKFUTW WORD. 
( 1) Altltud#(Ft.) 

( 2) Valoolt/Ut/aao) 

( 3) Alpha (dagraaa) 

( i) Laminar Fa o tor# 

( 5) Daflaatlon Faotora 
( 6) Local Maoh Nuabar 
( 7) Rajnoldf Nuabcra * 

( I) Rafao 

( 9) Local Fraaaura(lbf/ft/a*i) 

(10) Cp 

(11) Q Dot (Itu/ft-Sq-aao) 

(12) T Vail (Hag F) 

(13) Control fturfnoa Daf. (Dag.) 


J " 




/ tr A3KC i'N. 110 

\ 


1 YM 

• .4 

IT IN 11 : TIN 


Q- 


• ( A 


OMAX 

KJUTMI 

DELTA 

MAXTMX/10 



ORIGINAL pAnp 

0F «*>« Wauty 


CALL TEXEGM(AlO) 

CALL BASAIPfL/CJTD') 
CALL MIXALP( 'STANDARD') 


comorAabs/ label 


cokmAibojo/ BDTPNT.TRANME, TCASR, IHCOPI 


CCHOi/AIDAT/ XOR, ESTEP, LAX IS 


COMCM/ARRAT/UCI ,TIME1 ,BRA1 , ALPHA1 , AELZ1 , 


HFACZ1 , XMKi , RH1 ,CF1 ,TV1 ,TWF1 , 


QT,TT,TV2,QC1T 

1 

CCMO/fcAX/QM'VETAM, AKLZM, H7ACZM, PKM, CFFM, * 


TRKFM 

i 

C(MOI/lPLOTr/J1lCD, n, NT, NT2, NFF, FCT 


COMM m,IBUF f r?U 




INITIALIZE MC PS TERMINAL 
INTOtrACt, AND BRANCH TO 

proper puce 


INITIALIZE MICROFILM INTERFACE 
AND BRANCH TO PROPER PLACE 



CAM. AXSPLT(O.D,OKAX,N.O,ECR, ESTEP, EAEIS) 
CALL IINDIL0(TZ,ZZ,NT,1,1) 

| CALL liNUYt 0(T2,VZ,NT,2,1) 


j 


TKS 


N ALPHA ^ 0 x 


NO 


CALL IWI'TL0(r,ALPA<»T,NT,3,1) 


40 
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LBBiD(xrAX,>,a.},*.a) 

mmo(tTiiiR,M,6.a) 



HU 

call 


CALL KIDTL0(T0«’.,fK1,R t 9,)) 
CALL KK2YL0(T IKK 1 , CF1 , IT , 1 0, J) 


IHCOn* 0 


C CALL Dltn(O) 



UGTLT s 0 


V*° 
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u 
f *> 

^ TTIUE 

j^CALL j 
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■ I 

r CAM. imi»TUl(TIKKI,|l !.l, JT,7, j) 
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AMDEF 1 


CALL MHDILC(TPai,DaTAT,n,11,)l 


CALL HKBU0{TD«1,nW1,rT,ia,)) 
CALL HKCyL0(TI>Z1,TKEri,XT l ia,J) 
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CALL KDCOfT 
CALL RA5I 
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T .*? /mi > o.ooi 


r. , 
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< ITHICK^ 0 


^ TA'JI 
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| SBBMVTINE BUNU>(XO, XO, XL, TO, TO, Tl, XAWUT. TAJUUT , NFOINT, lELAfl, IE ASS) 
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1 
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; 
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iBiir(2oo),BDtmU), 
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z± 


common m, our, iru 

cowon A) now/ SDrm.TNANME, icase, ntcorr 



CALL X I)IT AX 

call m:a 

CALL TANO(O.O) 



CALL ENCPL(O) 
CALL I0UARO0) 




C tALL XlllTAX 1 

_g.ALL TANO (Q.O) j 


! 


i 


CALL CNAfHUO, XD.YO, YD) 

CALL CNID<2,2) 

CALL EWKA»(N, DATE, TIME) 

call iixicirr(o.io) 
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CALL MKtQIIT(O.ia) 
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CALL EXIT 
END 


It SUBROUTINE AIR62 (HGH.T.RHO.P.U) 

SIC H-ALTITUDE 

3*C N-LENCTH Of TABLE (SET TO 12 FOR 1959 MODEL! 

4IC HB-ALTITUDE BASE 

BtC TAB -BASE TEMPERATURE 

SIC GLMB* TEMPERATURE GRADIENT 

TIC RMOB-BASE DENSITY 

BtC T» TEMPERATURE 

9»C RHO*DENSTTY 

lSlC P -PRESSURE 

11»C U-UELOC1TV OF SOUND 

late C-FRAUITY 

13IC CIUEN H,N,HB( 1 )TO HB(N ) ,TMB( 1 ) TO TMB(N),GLMB( 1 ) TO 

14IC GLMB(N >,RHOB(l ) TO THOB(N) 

IStC COMPUTE T,RHO,P,U,G 

IttC 

l?t DIMENSION HB(22).TMB(22>.CLMB<22>.RH0B(22> 

IttC DIMENSION HB(22>,TMB<22>,GLHB<22>.RH0B(22> 

19t DATA (HB(I). 1*1. 22)/*., 36*89. 24,65616. 8.1*4985. 88. 154199. 48. 

2« l 217*6*3. 68.2M131. 24, 2S9186. 36. 291153. 22. 323**2. 64, 354753. 29. 

ait 23864D6. 18, 48S78D.86.512B45. 95,543215. 24, 6*5268. 39.728243. 79, 

22 t 3939894 .72,1234645.71.1 52*799 .26.1798726.44. 2*68776 .31/ 

231 DATA <TNB< I ), 1*1,22 1/518.67, 21389.97, 41 1 .57,2*487.17,454.77, 

241 12*325. 17.374. 17.469. 17,649. 17, 1729.17, 1999. 17,2179.17,2431 .17, 

26i 22791.17.3295.17,3889.17.4357.17,4663.17,4861.17/ 

26t DATA <CLHB( I ), 1*1,22 )/-3.S662E-3,0. ,5. 4864E-4, 1 .S362E-3.0. ,-1.0973 

271 lE-3,-2.1946E-3.*.,1.6953E-3.2.8343E-3.5.6867E-3.1.1444E-2.8.635E-3 

28 > 2,S.?744E-3,4.*61E-3,2.927E-3,2.381 IE-3,2.*152E-3, 1 .63S5E-3, 1.1*1 IE 

29t 3-3.7. 3298E-4.*./ 

3*1 DATA <RHOB( 1 ), 1*1,22 )/2.3769E-3.7.*612E-4, 1 .7*82E-4,2.5661E-S, 

31 t 12.7698E-6,1.473SE-6.4.8719E-7,3.8826E-8.6.15*8E-9,9.6511E-1*. 

32t 21 .9*71E-1*,4.7266E-11 .3.5624E-12.2.2488E-12, 1 .SS92E- 12,8 .4345E-13, 

33» 33.*346E-13.6.956E-14,1.26*8E-14,3.*599E-15,9.0*31E-16,3.*463E-16/ 

34 t DATA N/22/, CO/32. 174*485/. RE/29898908./ 

3St 12 H*(RE*HGR)/(RE+HCR) 

361 C*C0<((RE/<RE+HGMmt2> 

37t DO 1 1-2. N 

39t If (H-HB( I ) >2, 1,1 

39 i 1 CONTINUE 

4«t GO TO 5 

4U 2 IF(CLMBd-l) >3.4.3 

421 4 RMO-RHOBd-1 )«EXP(-(H-HB(I-1 ))*G0/(1716.4827*TMB(I-1 )) ) 

431 T-TMBCI-1) 

441 GO TO 6 

4Si 3 T-TMBd-1 )*CLMB(I-l )*(H-HB(1-1 ) ) 

46t RHO*RHOB( I-l >*EXPt-d .*♦00/ ( 1716.4827* CLMB d-1 ) ) »AL0C(T/TMB(I-1 ) 

471 11) 

48* CO TO 6 

49* 5 T*TNB(N> 

5** RHO-RHOB(N)SEXP(-(H-HI(N) )XGO/( 1716 . 4827*TMB<N ) ) > 

61* URITE (6, 1*0* ) H.HCM 

52* 1*** FORMAT (2E14.5) 

53* 6 U*SORT ( 1 . 4* ( 1716.4827XT ) ) 

54* P*RHO*( 1716.48278T ) 

55* RETURN 

"6* END 
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ATAS4 


II SUBROUTINE ATAS4 (H,T, DENS,P, SPS > 

81 REAL LAB 

3< DIAENSION PB( 14 ),ZI( 19 >,PK(6,S >,RH0K<6,3 ),TK(6,5 >,TAB< 

41 114 ), LAB( 14 ) # TB< 14 ) 

St DATA PBASE /6.231991759E-5/ 

St DATA ZI /19832. 1 , 17853.3,28999. ,49999. ,83994 

?t X. ,9.E4, 1 .E5, 1 • 1ES, 1 .2E5, 1 . 5E5, 1 .6E5, 1 .7E5, 1 . 9E5,2.3E5, 
81 X3.E5,4.E5,5.E5,6.E5, 7.E5/ 

9* DATA PIC /I .6871582E-2, 

lit X-1.1425176E-4,-1.36l2327E'9,7.3624145E-14, 

lit X-l . 0899315E-17, 3.3946432E-22, -7.9919777E-2, 

181 X-8 . I 946438E -5 , -5 . 5S22383E - 9 . 3.111 6969E- 1 3 , 

131 X-l. 66878276-17,3. 8319351E-22,9.8414277E-1, 

14t X-2. 697691 7E-4, 8. 5227541E-9, -3. 962826 3E- 13, 

15 1 X1.9146471E-17,-1 .9264318E-22, 1 . 14118495E1, 

IS t X-4. 1 1 497477E-4, 1 . 336648S5E-8, -3.595 1897SE- 13, 

17t XS • 199972S4E-18, -2.899558946-23,9.99324461, 

181 X-2.S8298177E-4,3.76139346E-9,-4.29887236E-14, 

19t XI. 691821486-19,-1. 925889276-25/ 

891 DATA RHOK /I .33921 17E-2, 

81 t X-8. 85 9296 46 -5, -4.2143856E-9,5 .95175576-13, 

82t X-3.9744789E-17, 7. 87712736-22,1.26671226-1, 

83t X-l .33731 476-4, 2. 966737 IE-9, 2. 3396 1 896-13, 

841 X-3.2562583E-17, 7.99352896-22,9.27512666-1, 

25* X-l .43496796-4,-2.82717366-9,4.74898926-14, 

2St XI .88632466-18, -4.278241 IE-23/ 

27t DATA TK /2 .966787762, -6.77319916-3, 

88t X8. 4619895E-7,-l .7894949E-19, 1 • 14S145E-14, 

89t X-2. 4898788E-19,2.6892151E2, 4.39753S2E-3, 

39t X-8.9159672E-7, -2 .8929791E-1 1,5 .97248566*15, 

311 X-l . 14993726-19,3.796455762, -3 .28S8965E-2, 2 .96456366-6, 
32 t X-4. 3283944E-1 1,-5. 7S972426- 17,8. 2924583E-21, 

33 1 X2. 944798E 1,2. 97698384E -2 , - 8 . 63938789E -7 , 

341 XI .663924176-1 1 , -9. 39976 18S6- 17,-4.999951986 -22, 

35l X-4.988659S362, 3.921372816-2,-4 .9S189691E-7, 

36t X-3. 262198546-12, 9. 66659364E-17,-4. 788442796-22/ 

37t DATA TAB / 189.65,219.65,269.65,369.65, 

38t X969.65, 1119.65,1219.65, 1359.65, 1559.65,1839.65, 

39t X2169. 65, 2429. 65, 2599. 65, 2799. 65/ 

49t DATA LAB / 3. 6-3, 5. 6-3, 19.E-3,29.E-3, 15.E-3, 

411 X19.E-3,7.E-3,5.E-3,4.E-3,3.3E-3,2.6E-3, 1 .7E-3, 1 • IE-3, 

42t X1.1E-3/ 

43i DATA PB /. 172244361E-4, .315971712E-5, 

441 X.774389897E-6, .2659771 11E-6, .535849383E-7, 

45t X.391284945E-7, .29591 1 1 17E-7, . 178715656E-7, 

46 i X. 739258 171E-8. .2995731 16E-8, .439456696E-9, 

47 t X. 1 17315489E-9, • 379198961E-19, • 1281 15339E-19/ 

48 t Z-H/3. 289833 

49t IF (Z> 5,9,6 

Site N LESS THAN 9 SET TO 9 

51t 5 Z-9.9 

52 t GO TO 9 

53* 6 IF (Z-799999.) 9,9,8 

54 tC H GREATER THAN 7B9999 PETERS SET TO 799999 
55 * 8 Z-799999. 

56* 9 N-9 

57* 19 N-A41 

58 1 IF (Z-ZKH)) 49,49,19 

591 49 IF (N-5 >29,29,65 

69* 29 Z2-Z8Z 

Sit Z3-Z28Z 

681 Z4-Z28Z2 

631 Z5-Z3SZ2 

64 I C STD COAPUTATION FOR TEAPERATURE 
SSi 65 IF (N-5) 79,79,89 
SilC H LESS THAN 83994. AETERS 

67* 79 T- TK ( l ,N )4TK(2,N )t2*TK < 3,N )IZ2*TK( 4,N )tZ3*TK ( 5,N )tZ4 + 

68 1 1TK(6,N)8Z5 

69 t T -Ttl.8 

791 GO TO 119 

711 89 IF (N-6) 85,82,85 

72 i C T FOR 83994. THRU 99999. 

73t 82 T-189.6S81.8 

74* GO TO 119 

75tC T FOR GREATER THAN 99999 

76t 85 A-A-6 

77l T«TAB(A)4LAB(A)f (Z-ZICN-l ) ) 

78t T -Til. 8 
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ATAS4 
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ATAS4 
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ATAS4 

ATAS4 

ATAS4 

ATAS4 
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ATAS4 

ATAS4 

ATAS4 

ATAS4 

ATAS4 

ATAS4 

ATAS4 

ATAS4 

ATAS4 

ATAS4 

ATAS4 
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ATAS4 

ATAS4 
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ATAS4 
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ATAS4 

ATAS4 

ATAS4 

ATAS4 

ATAS4 
ATAS4 
ATAS4 
ATAS4 
ATAS4 
ATAS4 
ATAS4 
ATAS4 
ATAf 4 
ATAS4 
ATAS4 
ATAS4 
ATAS4 
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791 119 IF (N-5 > 139.139.1S9 

99 i C H LESS THAN 83994 DETERS ATI4S4 

111 139 CON* 19.9 ATNS4 

•21 IF (N-3 > 13S.13S.132 ATNS4 

•3t 132 CON-9. 8066SE-4 ATNS4 

•41 13S P -CON »EXP<PK<1,N>*PIC(2,H>*Z*PICC3.N)*Z2*PK(4.N>8Z34 ATNS4 

•Si 1PKIS.N >>Z4+PK(£,N)t25 > ATNS4 

Ml P -PF208.9 ATNS4 

•7l GO TO 399 ATNS4 

•Si 1S9 IF (N-6 >160, 169, 179 ATNS4 

89i 169 P*PIASE8EXP< <-1.373381523E12*(Z-83094. >>/( 189.65* ATNS4 

99l 1 <6344869. *Z ><<6344869. *83994. )>> ATNS4 

•II P -PJ208.9 ATNS4 

92 ■ GO TO 399 ATNS4 

931 179 N*N-6 ATNS4 

941 P*EXP< ALOGCPKH > >♦< 1 . 37339 1523E 12/ ( LHB(R >*(6344869. +Z )* ATNS4 

951 1 (6344869. »Z1 (N-l ) > > )«AL0C(TN8(N >/(TNB<n (♦(LNStHlUZ-ZI (N-l >>>>)) ATNS4 

96i P *Pt298.9 ATNS4 

971 399 IF (N-3> 316,316.320 

981 316 DENS* < 1 . 16799729 )(EXP( RHOK( 1 ,N )*RH0K(2,N KZ+RHOK ( 3.N >( ATNS4 

991 1Z2»RH0IC(4.N)*Z3*RH0K(S,N>*Z4«-RH0IC(6.N>*ZS> ATNS4 

1991 DENS *DENSt. 09194 ATNS4 

19H GO TO 369 

1921 329 DENS*S.82S871E-4*P/T 

1931 369 SPS*49.92U77*S0RT(T> 

1941 RETURN ATNS4 

1951 END ATNS4 
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SUBROUTINE BINTRP< X, V, I TABLE, ZZ ) 


It 
2iC 
3tC 
4tC 

5lC THIS SUBROUTINE CONPUTES PRANDTL N0.,LEUIS NO., AND HSUBD 
SlC 


7iC 

8tC 


91 


COnnON/ERRFLG/NERROR 

101 


CONNON/FFNUN/MSUBD, ALEU 

Ilf 


DINE NS I ON XJ( 190 ), YI (26 ) 

12t 


EQUIVALENCE (A, ZI) 

13t 


IF ( IFLAG.EQ. 1 ) GO TO 9 

14! 


DO 1 J-1,40 

IS! 


DO 1 1-1,26 

16t 

1 

Z( I, J >-2I< J, I ) 

17f 

9 

IF ( ITABLE.EQ. 1 ) GO TO 10 

18i 


IF(ITABLE.E0.2)G0 TO 20 

19s 


IFCV.LT.l.E-6) Y-l.E-6 

20! 


N-66 

21! 


NW-90 

221 


N- 19 

231 


NN-26 

24! 


IF ( (Y .LT. YI (N ) ) .OR. 

25t 


NN-6S 

26t 


G n TO 30 

27! 

10 

N-l 

28i 


IF ( V.LT. 1 .E-8 ) V-l.E-8 

29! 


NN-40 

301 


N-l 

31! 


NN-10 

32! 


IF ( (Y .LT. YI (N ) ) .OR. 

33! 


HN-0 

34t 


GO TO 30 

3St 

20 

N-41 

36! 


NT1-65 

37i 


N-ll 

38! 


NN-18 

39t 


IF ( (V .LT. VI (N ) ) .OR. ( 

401 


NN-40 

41! 

30 

N-NM 

42t 


DO 40 I -N, NN 

43> 


IF ( Y.LT. YI ( I ) ) GO TO 5< 

44t 

40 

CONTINUE 

45! 


I -NN 

46! 

50 

N-N^l 

47 « 


DO 60 J-N, NN 

48! 


IF C X .LT. XJ(J) ) GO TO 

49 1 

60 

CONTINUE 

50i 


K-NN-NN 

51! 


GO TO 80 

52 1 

70 

K-J-HN 

53i 

80 

ZNIN- Z(I-1,K-1 >*(<X-XJ(. 

54! 

1 )) 

55 1 


ZUAX-ZU.K-l >♦< ( X-XJ ( J-l 

56! 


zz*zniN.((v-vi(i-i))/(vji 

57! 


IF ( ( [TABLE .EG. 1 ) . AND. ( ZZ 

58 ! 


IF ( ( I TABLE . £0 . 2 ) . AND. (ZZ . 

59t 


IF(<ITABLE.E0.1).AND.(ZZ, 

•0! 


IF< (ITABLE. EO. 3 >. AND. (ZZ. 

61! 


IF((ITABLE.EQ.3).AND.(X.( 

62i 


IFLAOl 

63i 


IF (ITABLE .EO. 1) HSUBD* 

641 


IF (ITABLE .EO. 2) ALEU-2 

65i 


RETURN 


A ( 1 040 ) , 21(40,26) 


.LT. XJ(N) )) CO TO 300 


(X .LT. XJ(fl) )) GO TO 310 


.it. xjcnn ) go to 3ie 


70 


BINT00I5 
BINT0020 
BINT 0025 
BINT0030 
BINT0035 
BINT0040 
BINT0045 
BINT0050 

BINT0055 
BINT0060 
BINT0065 
BINT 0070 
BINT0075 
BINT0080 
BINT0085 
BINT0090 
B1NT0095 
BINT0100 
BINT0105 
BINT0119 
BINT01 15 
BINT0120 


BINT0130 

BINT0135 

BINT0140 


0. nzz-0. 

,5))ZZ-.S 


BINT0150 
BINT9155 
BINT0160 
“i ,K-1BINT0165 
BINT0170 
) BINT0175 
BINT0188 


24000 . ).AND. (ZZ.GT. .734 ) )ZZ- .734 


66 i 300 URITE (6,900 ) X,V 

671 900 FORHAT ( 72H0 8888? NEGATIVE TENP OR PRESS RATIO LESS THAN 10-6 

68 t 1UBR. BINTRP . T-,E12.6,6H P/P0-,E12.6> 

69 t NERROR-1 

70t RETURN 

71 t 310 URITE(6,919) X,V 

72 I 910 FORNAT ( 73N0 ttttt NEGATIUE TENP OR DENSITY RATIO LESS THAN 10 
73! 1SUBR. BINTRP. T-,E12.6,1SH DENSITY RATI0-,E12.6 > 

741 NERROR-1 

76i RETURN 


BINT0185 
BINT0190 
IN SBINT0195 


B I NT 021 0 
8 INBINT0215 
BINT6220 


H-16 


i 


76i 
77i 
711 
791 
tit 
9H 
Mi 
83 1 
•41 
•St 
Mt 
•71 
Sti 
89t 
tit 
91 1 
92t 
93t 
94t 
95 1 
96t 
97t 
98t 
99: 
ltii 
ltlt 
182i 
183t 
194t 
185i 
186t 
197t 
188i 
199t 
Hit 


1241 


1411 

142t 

1431 

1441 

M5i 

146t 

M7t 

148t 

149t 

1S9I 

1S1I 

1521 

153* 


DATA ( X J ( J ) , J • 1 , 98 > /9 . . 5489 . , 5768 . . 6 1 29 . , 6489 . , 6849 . , 7299 . , 7569 . , 8 1 NT9239 

17929. . 8289. .8649. .9899. .9368. . 9729. . 10989. . 19449. . 19899. . 1 1 169. . t 1BINT0235 

2529 . . 1 1 889 . . 1 2249 . . 12699 . . 1 2969 . . 1 3329 . . 1 3689 . . 1 4849 . . 1 4408 . , M769B 1 NT9240 
3. # 15129. , 15489. , 15849. , 16290. , 16569. , 16929. , 17289 . , 17649. , 18099. , BIHT024S 

4 1 9899 . . 23499 . . 36999 . . 9 . . 3699 . . 4599 . . 5499 . . 6399 . . 7299 . . 8 1 99 . . 9999 . , B I NT02S0 

59999 . . 1 9899 . . 1 1 799 . . 1 2699 . . 1 3598 . . 1 4 499 . . 1 5399 . . 1 6299 . . 1 7 1 89 . . 1 800B I NT0255 

69. . 18909. . 19899. #29799. #21699. .22590. ,23499. ,25290. , 0.0 , 3609. ,BINT0269 

7 4599. , BINT0265 

8 5499 . , 6399 . , 7290 . , 8 1 99 . , 9900 . , 9909 . , 1 0890 . , 1 l 709 . , 1 2699 . , 1 3590 . B I NT0279 

9, 14400. , 15390. , 16200. , 17100. , 18009.. 18990., 19890., 20799. , 21690. ,BINT0275 

1 22509 . . 23400 . . 25200 . / B I NT0280 
DATA (VI(I), 1-1,26)/. 9S8E-8, .958E-6, .958E-5 , .958E-4 , .9S8E-3 BINT0285 

1 , . 958E-2 , . 958E -1 ,.958 , 9.58, 19.16 , 1 .E-6, 1 .E-4, 1 .E-3, 1 .BINT0290 

2E-2, l.E-1 , 1. , 10., 100. , 1 .E-6 , l.E-4 , l.E-3 , t.E-2, 1.BINT0295 

3E-1 ,1. , 19. , 109. / BINT0300 

DATA (A(K),K-1, 290 ) / 9., 1.8EO, 3.82E+3, 6.9E43, 

1 8.2EO , 19.6E43 , 11.43E43 , 12.0EO , 12.23E43 , 12.26E + 3 , BINT0310 

2 39S12.45E+3 , 9. , 1.S6E+3 . 1.621EO , 1.779E43 , 2.138E43 ,BINT0315 

3 2.857E43 , 4.122E43 , S.98EO , 8.18E43 , 10.11EO , 11.3SE43 , BINT0329 

4 11.96EO , 12.23€>3 , 12.35EO . 12.49 E^3 , 12.43EO , 12.44£4|JNT0325 

S3 , 23B12.45E+3 , 0. , 1.522EO , 1.553E^3 , 1.697E^3, 1 . 727EO, B INT0330 
61 .97E*3 , 2.432E43, 3.219E43 , 4.422EO , 6.037E43 ,7.869EO ,9.5BINT0335 

762E+3, 19.83EO, 11.6E43 ,12.01EO , 12.22E43 , 12.33E43 , 12. 38EBINT9340 

80 , 12.42E43 , 12 . 43EO , 12.44E43 , 19112. 45EO , 9. ,1.411E*3, BINT934S 

91.489EO , 1.531E43 , 1.579E^3, 1.664EO , 1.819E43, 2.092EO ,2.SBINT03S0 
143EO , 3.235EO , 4.21SE43, S.480EO . 6.945E43 , 8.432EO , 9.BINT03SS 

2745E^3 , 10 . 748EO , 11.428EO , 11.85E43 , 12.99EO , 12.23EO .BINT0360 


3 12.32EO , 12.37E43 

4^3 , 11112. 45EO , 9., 1 . 002E + 3 
5 1 . 490E43 , 1 . 593EO , 1.690E43 
3.683EO , 

10.994E+3 
12.22E43 


6 2 . 979E^3 

7 9 . 1 45E^3 

8 12.12E43 

9 12.40E43 


12.40EO , 12.41E43 , 12.42E43 , 4412.44EBINT036S 
1.256E43 , 1 . 397EO , 1.472E43, BINT0370 
1.847EO , 2.091E43 , 2.4S8E43, BIHT037S 
4.580E43 , 5.649E43 , 6.831EO , 8 . 039E43, B INT0380 

. 10.82E43. 1 1 . 36C+3 , 11.72E^3. 11.96E+3, B1NT0385 
12.39E43 , 12. 34EO . 12.36E43, 12.39E43 , B1NT0399 


12.42E43.12.43E43 , 3I12.44E43 . 5I12.45EO/ 


DATA ( A(K ),K*201 # 280 )/ 0., .461E + 3. .7198E43. 


B1NT9395 


Hit 

1 

.9597EO. 

1 .184EO 


1 . 322EO, 

1 .414E.3 

0 

1 .488E43, 

1.563E43, 

B I NT, 405 

112i 

2 

1.658EO. 

1.791E43 


1.978E43, 

2.236E43 


2.585E43 , 

3.041E43 

, B IMT0410 

113i 

3 

3 . 61 1E+3, 

4.304E43 


5.102E+3. 

5.989E43 


6.922E43 , 

7.86E43 

, BIMT0415 

1141 

4 

8.754E43, 

9.550E43 


10.23E43, 

10.78EO 


11.21E43 . 

11.53E43 

, BIHT0420 

U5t 

5 

11.76E43, 

U.95E43 


12.08EO, 

12.18E43 


12.24E43 # 

12.29E43 

, BIMT0425 

116t 

6 

12.33E43, 

12.35E43 


12.38E43, 

12.39E43 


3I12.45E43, 

0. 

,BINT0430 

1171 

7 

1.660E42, 

2.869E42 


4.483E42, 

6.378E42 


8 324E42 , 

1.009E43 

, B1NT0435 

118t 

81 

, . 1 154E+3, 

1.269E43 


1 . 364E+3, 

1.452E43 


1.543E43 , 

1.659E43 

, BINT0440 

1191 

9 

1.784E43, 

1.954E43 


2 . 169E + 3, 

2.437E43 


2 . 769EO , 

3.168E43 

0 

1291 

1 

3.636E43, 

4.176E43 


4.772E43, 

5.431E43 


6.125E43 , 

6.849E43 

, BINT0450 

121t 

2 

7.554E43, 

8.238E43 


8.992E+3, 

9.483E43 


10.01E43 , 

19.45E43 

, BINT0455 

122i 

3 

10. 83E+3, 

11.14E43 


11.49E43, 

11.60E43 


11.77E43 , 

U.99E43 

, BINT9460 

123< 

43*12. «SEO/ 







BINT0465 


DATA ( A(K ),K-281 » 360 )/ 0. 


5.466E41, 9. 779E ♦ 1 , 


125t 

1 

1.611E42, 

2.467E42 

0 

3.537E42 # 

4.777E42 , 

6.118E42 

, 7.475E42 

, BINT0475 

126t 

2 

8.783E42, 

9.973E42 

0 

1.106E43, 

1.207E43 , 

1.393E43 

, 1.401E43 

, BINT9480 

127i 

3 

1.505E43, 

1.622E43 

0 

1.7S6E43, 

1.914E43 , 

2.097E43 

, 2.314E43 

, BINT0485 

128i 

4 

2.56SE43, 

2.849E43 

0 

3.174E43. 

3.538E43 # 

3.937E43 

, 4.373E43 

, 0INT0490 

129t 

5 

4.840E43, 

5.341E43 

0 

5.842E43. 

6.365E43 , 

6.899E43 

, 7.417E43 

,BINT0495 

I30t 

6 

7.922E43, 

8.397E43 

0 

8.872E43, 

9.286E43 , 

11.60E43 

,2tl2.45£43 

,BINT0500 

13H 

7 

9. 

17.30 

0 

31.65 

53.20 

83.2 

,1.236E42 

, BINT0S05 

132: 

8 

1 . 743E42, 

2.359E42 

0 

3.069E42. 

3.85E42 , 

4.71E42 

, 5.59E42 

, BINT0S19 

133t 

9 

6.28E42 , 

7.089E42 

0 

8.291E42, 

9.182E42 # 

1.007E43 

, 1.096E43 

,BINT0515 

I34t 

1 

1.189E43, 

1.285E43 

0 

1 .388E43, 

1.501E43 , 

1.624E43 

, 1.760E43 

, BINT0520 

135t 

2 

1 . 910E-3, 

2.076E43 

0 

2.259E43, 

2.459E43 , 

2.689E43 

, 2.921E43 

, BINT0525 

136t 

3 

3.179C43, 

3.46EO 

0 

3.752E43, 

4.061E43 , 

4.387E43 

, 4.732E43 

, BIHT9530 

137t 

4 

5.98E43 , 

6.80E43 

0 

19.5E43 , 

12.45E43/ 



BIHT0535 

138t 

DATA (AIK), 

K-361,640)/ 

0., 13.88, 

22.48, 




139t 

1 

37.72 , 

59.42 

0 

88.47 

1.257E42 , 

1.713E42 

, 2.252E42 

, BINT0545 

1491 

2 

2.865E42, 

3.482E42 

0 

4.267E42 , 

5.926E42 , 

5.817E2 

, 6.612E42 

, BINT0550 


3 7.414E42, 8.227E42 , 9.031E42 , 9.879E42 , 1.074E43 , 

4 1.260E43, 1 . 360E ♦ 3 , 1.470E43 . 1.599EO , 1.720E43 , 

5 2.016E+3, 2.187E43 , 2.372E43 , 2.568E43 , 2.784E43 , 

6 3.255E43, 3.513E43 , 3.785E43 , 4.07EO , 5.6E43 

7 12. 45E+3, 211.35 , 1.125, 1.09 , .93 , .73 , .62 


1.164E+3 , B1NT9555 
1 . 863E+3 , B1MT0560 
3.91E43 , BINT0565 


10.5E43 
.56 


, BINT0570 


8 

.52 , 

. 509, 39t .50, 

1.35 , 

1.335, 

1.172, 

1.147, 

1.143, 

1.051, 

BINT0S80 

9 

.786, 

.699, .580, 

.549, 

.520, 

.505, 

281.50, 

1.35 , 

1.346, 

BINT9585 

1 

1.233, 

1.157, 1.154, 

1.130, 

.998, 

.750, 

.619, 

.578, 

.541, 

BINT9599 

2 

.521, 

.599,270.59, 

1.35 , 

1.355, 

1.355, 

1.213, 

1.167, 

1.144, 

BINT9595 

3 

1.114, 

.970, .762, 

.628, 

.588, 

.574, 

.535, 

.519, 

.509, 

BINT9600 

4251.50, 

1.35 , 1.355, 

1.391, 

1.327, 

1.227, 

1.180, 

1.160, 

1.195, 

BINT0605 

5 

.982, 

.808, .669 , 

.695 , 

.587 , 

.570 , 

.530 , 

.518 , 

.597 , 

BINT0619 

6231.59, 

201.35,1.404 , 

1.405, 

1.320, 

1.225, 

1.196, 

1.158, 

1.115, 

BINT061S 




{<- 1 / 


1541 71. #19 . .885 , .746 , .652 , .661 , .576 . 

1551 8 .567 . .S62.S9t.S6/ 

1561 DATA (A(K ).K*64x, 1666 >/ t. 35. 1.35. 1.467, 1.437, 1.424. 

1571 1 1.347. 1.26S, 1.215, 1.167. 1.131, 1.671, .975 . .866 . .756 , 

1561 2 .669. .618 . .584 , .558 . .546 , .526 . ~ ‘ 

1591 32S1. 35, 1.468 ,1.445 , 1.464, 1.449, 1.464, 


k 

i7it 


.567, . 536, *515 , 


.Sit, .562 , 181.5, 
1.327, 1.268, 1.235, 


16# i 

4 

1.158. 1.118, 

1.059, .979 

• 

.893 

6 

.805 , .726 

, .662 . .618 , 

16H 

5 

.585 

. .564 , 

.531 , .525 

6 

.518 

, 

161. 5#, .77 

, .772 , .732 , 

162i 

6 

.726 

. .71# , 

.681 , .68# 

6 

.69# 

6 

.7## , .716 

, .718 , .725, 

163i 

7 

.731 

.27B.734, 

.77 , .772 

6 

.75# 

6 

.735 , .734 

, .717 . .683 , 

1641 

8 

.683 

, . 694 , 

.702 , .712 

6 

.72# 

6 

.726 . .73# 

,268. 734,. 77 , 

165t 

9 

.772 

6.763 . 

.739 , .734 

6 

.729 

6 

.708 , .632 

, . 686 , . 697 , 

1661 

1 

.764 

. .712 , 

.719 , .729 

,261.734 

, .77 , .772 

, .77# , .747 , 

1671 

2 

.738 

6 .734 , 

.725 , .704 


684 


.691 , .7, .708, .716, .725, 

1681 

3 


261.734, 

.77 , .773 

6 

.773 


.762 , .744 

, .738 , .732 , 

160t 

4 

.723 

6 .765 , 

.687 , .689 

6 

• 7#3 


.71# , .721 

.729 ,258.734 

17#i 

5 

.77 

. .772 , 

.774, .773 

6 

.756 


.743 , .739 

, .732 . .724 , 

1711 

6 

.769 

, .693 , 

.684 , .696 

6 

.71# 


.72, .727,248 

734,. 77,. 772, 

1721 

7 

.775 

6 .777 , 

.772 , .758 

6 

.747 


.739 , .732 

, .724 , .715 , 

1731 

8 

.70# 

6 .69# , 

.696, .693 

6 

.707 


.714 , .725 

,228.734/ 

1741 

DOT# 

<*(K),K«1##1,1#4#)/ . 

77 

.772, 

.775, .775, 

.775, .771, 


1 .764, .753 , .744 , .746 , 728 , .72# , .1 

2 .687 , .687, .695 , .764 , .71# , .716 , 

3168.734/ 

END 


11 , .BVI » «IV1 i 

.722 , .787, .738, 


8INT#62# 
BINT 0625 

8INT663S 
8INT664# 
BINT#645 
81NT#65# 
BINT#655 
B 1NT666# 
81NT6665 
81NT#67# 
BINT0675 
•INT0680 

BINT#69# 
B INT0695 

B1NT#7## 

B1NT071# 

BINT#715 

8 I NT# 785 
81NT673# 
8INT8735 
8INT874# 


8-18 


1 

2 

3 

4 
6 

5 
7 
I 
9 

19 

11 

12 

13 

14 

15 
IS 

17 

18 
19 
21 
21 
22 

23 

24 

25 

26 

27 

28 
29 
39 

31 

32 

33 

34 

35 
38 

37 

38 

39 
49 

41 

42 

43 

44 

45 
48 

47 

48 

49 
59 

51 

52 

53 

54 

55 

56 

57 


SUBROUTINE CHEEUVtNSTEP ) 

DIMENSION ALPHA (25 ). 5(25). DTENP( 109 ).C<25 >.ZPLUS(25 ).ZMINUS<25 ). 
1 A(2S ) . B ( 25 ) . C ( 25 ) , 0(25) 

DIMENSION TSAUE ( 25 ) . RHO ( 25 >,CP(25 ).COND< 25 ).ZZ(S9 > 

COMMON/FLUFLD/ PHIL( 16 ) . TU.MU. REST < 34 > 

COMMON/ THICK/ DTEO. TDOT # TEO. OC. OR. HCOUT. ENISS. HREC. 

1 DT.CTIME, PL. DTO. DTI. TE02 

C0MN0N/MANDV/INAX.DX(2S>.T<2S>.NPFLAG<2S>. NTIUZ. TSVNK. ENISIN. HCIU 
1 .TCSF.TIU2(59 I .HCIUZ( 59 )» TCAS2( 59 ).TSINK2(58 ) . TUI , TS1NK . TGAS.HCIN. 

2 IFIRST.DELT 

FORMATC 1H9. 13X. 36H!SSNEUTON-RAPHSON FAILED TO CONUERGE. 

1 4SH ON A UALUE FOR THE OUTER WALL TEMPERATURE AT. F18.3. 

2 UH SECONDS!!! ) 

2991 FORMATC /1H .SHK<1>*. E19.3.8H DX(l)*,E9.3. 2X, 4MEPS*. Ell. 3. 2X. 

1 6HNSUBC- . E9.3. 2X. 5HTREC-. E9.3. 2X. 6HHUALL* . E9.3. /1H . 
2SHTREC-. E 19.3. 2X. 6HTUALL*. E9.3. 2X, SHSTEP*. 12) 

3999 FORMAT( 1H9. 13X. 36H!!!NEUT0N-RAPHS0N FAILED TO CONUERGE. 

1 45H ON A UALUE FOR THE INNER UALL TEMPERATURE AT. F19.3. 

2 11H SECONDS! St ) 

3991 FORMAT (1H .8HKUMAX)-. E9.3. 2X. 9HDX(IMAX)-. E9.3. 2X. 6MEPSIN*. 
1 E9.3. 2X. 3MHC*. E19.3. 2X. 6HTSINK*. E9.3. 2X. 5HTGAS*. E19.3. 2 
2X. 5HSTEP-. 12) 

4999 FORMAT (/IH .29X. 17HLAST TEMPERATURES/ ( 19X. 5E29.5/)) 

COMMON/ERRF LG/NE RROR 
COMMON/FLAG/ IDEAL 
DATA STEUIE/ .4769E-12/ 


SET INITIAL CONDITIONS. 


IF (IFIRST .EQ. 0) GO TJ 159 
DT • DELT 

159 IF (NSTEP .NE. 1 ) GO TO 300 
TUOLD • TU 
DO 175 I ■ 1. IMAX 

175 TSAUE < I ) ■ T C I ) 

CALL STOCK ( I MAX . MPF LAG . T . RHO . CP . COND ) 

IF ( IDEAL.GT .0 ) CO TO 176 

CALL MOLIER(HU. PL. 2. TU. 22. SS. RR. GR ) 

IF< IDEAL.GT .0 > CO TO 176 

CALL NOLIERCHREC. PL. 0. TREC. Z2. SS. RR. GR ) 
IF ( IDEAL .GT. 0 ) GO TO 176 
GO TO 180 

176 HU* . 24STU 
TREC-HREC/.24 

189 HC1N-HCIU 

TGAS-TGSFM59.7 
TSINK-TSVNKM59.7 
IF (NTIUZ . LE . 0 ) GO TO 903 
1 F ( HC I UZ ( 1 ) . LE . 9 • ) GO TO 901 
CALL TBLINCCTIME.T1UZ. HC1N.HCIU2.2.Z2. NTIU2 ) 
991 IF (ABS(TGAS2( 1 ) ).LE. .0001 ) GO TO 902 

CALL TBLIN<CTinE.TlUZ.TG.TGA$Z.Z.ZZ. NTIUZ) 
TGAS-TGM59.7 

902 IF ( ABS ( TSINKZ (1)1- LE • - 0001 ) CO TO 903 

CALL TBLIN(CTIME. TIU2.TS. TSINKZ. Z.2Z. NTIUZ) 


88 

69 

89 

81 

82 

83 

84 
65 
88 

87 

88 
89 
79 

71 

72 

73 

74 

75 

76 

77 

78 


TSINK-TSM59.7 
993 CONTINUE 

OC • HCOUT ! (HREC - HU) 

OR • EMISS ! STEUIE ! TU! t 4 

00 • -OR ♦ OC 

01 • ENISIN ! STEUIE ! TUI!!4 - HCIN S (TCAS -TUI) - ENISIN ! 
1 STEUIE ! TS INKS! 4 


STEP ONE OF RUNGE-KUTTA INTEGRATION. COMPUTE OUTER AND INNER UALL 
MEAT FLUXES USING INITIAL TEMPERATURES AT TIMEd) OR RESULTS OF 
LAST SERIES OF CALCULATIONS AT TIME (N ) . 


399 CONTINUE 


COMPUTE COEFFICIENTS OF TRI-DIAGONAL MATRIX OF INTERIOR 
TEMPERATURES AT TIHE(NM/2) USING RESULTS OF PREUIOUS CALCULATIONS 


8-iy 


79i CO) • RHOO) S CP<1 > t 0X0) 

Ms GOUT • CO) 

III ZPLUSO) • 2.ICONDO )BC0ND(2 )/(C0ND(2 )tDX( 1 ) ♦ CONDO )*DX<2) ) 

Ml ZNINUSO >• I. 

831 FOUT • ZPLUSO )/G(l> 

141 NO) • I. 

•Si IO) • 1. 4 FOUT * (DT/2.) 

Mi CO) • -FOUT 8 (DT/2. ) 

•71 DO) « TSAUEO) ♦ 00 / GOUT t DT/2. 

Ml MINUS • MAX - 1 

•Bl IFIINAX.E0.2) CO TO 411 

8*1 DO 4M I • 2, MINUS 

111 ZNINUSU) • 2. • CONDO >*C0NDd-l ) / (C0ND( 1-1 )IDXO) 4 CONDO X 

82i 1 DXO-D) 

I3i ZPLUSO) • 2. t C0ND(I)*C0NDO4tl / (CONDO*I)tOX( I ) 4 CONDO)* 

941 1 DX(I41)| 

8St CO) • RHOO) * CPd) * DX(I) 

86i A(I) • - (DT/2.) t ZNINUSO )/G(I ) 

871 DO) • 1. 4 (DT/2. ) t (ZPLUSO) 4 ZNINUSO))/ CO) 

881 CO) • -(DT/2.) t ZPLUSO )/C( I ) 

881 DO) • TSAUEO) 

lMi 481 CONTINUE 

llll 411 ZNINUS(MAX) • 2. * C0ND( INAX )tC0ND( MAX-1) / (CONDI INAX-1 )*DX( INA 
182i IX) 4 CONDI INAX )*DX( MAX-1)) 

1131 ZPLUSI MAX ) • «. 

1141 GONAX) • RHOOHAX) • CPI MAX) * DXOHAX) 

18S < CIN • C(INAX) 

1161 FIN • ZNINUSI MAX ) / GONAX) 

1171 AONAX) • -FIN * (DT/2.) 

llll UMAX) • 1. 4 FIN « (DT/2.) 

1181 CONAX) • 8. 

llll D(INAX) • TSAVEl MAX ) - DT/2. ( OI/CIN 

line 

U2iC 

113IC SOLVE FOR NEW TENPERATURES BY GAUSSIAN SUMMATION METHOD. 

114IC 

USlC 

llll CALL GAUSS (A. I, C. D. ALPHA. S. T. MAX) 

117IC 

llliC 

UHC COMPUTE THE DERIVATIVES ASSOCIATED UITM NEU TEMPERATURES. 

IMiC 

line 

1221 CALL DRIVELIDTENP, C. ZPLUS, ZNINUS. 00. 01. T. NSTEP, MAX) 

123IC 
184 1 C 

12S>C AT STEPS 1.2. AND 3 OF RUNGE-KUTTA, USE LATEST TEMPERATURES TO 

12SIC COMPUTE NEU OUTER UALL TEMPERATURE. 

187IC 

line 

188IC IF NSTEP • 4, RUHCE -KUTTA IS COMPLETED. CALCULATE FINAL NODE 
138IC TENPERATURES. UALL TEMPERATURES. AND HEAT FLUXES 

lane 

132IC 

1331 IF (NSTEP .EO. 4) CALL FINALTI TSAVE. DTENP. MAX, 4, DT. T) 

1341 TNEUT1 -Til) 

1351 LIN • 2 * MAX 4 1 

13SIC 

137IC 

13HC AT STEP 4. USE TEMPERATURE (N ) PLUS DERIVATIVE COMPUTED AT STEP 3. 

i3nc 

148IC 

1411 IF (NSTEP .EO. 3) TNEUT1 • TSAUEO) 4 DTENP! LIN ) t DT 

1421 TERN • 2. I CONDI 1 >/ DX(1) 

1431 AOUT • ENISS * STEVIE 

1441 HCOUTI • HCOUT t ( HREC-HU )/ (TREC-TU ) 

1451 BOUT • HCOUTI 4 TERN 

1451 COUT • HCOUTI t TREC 4 TERN * TNEUT1 

147IC 

i4nc 

149IC SOLVE FOR NEU OUTER UALL TEMPERATURE IV NEUTON-RAPHSON METHOD. 

is«ie 

line 

1681 CALL NEUTIAOUT. BOUT, COUT, TU. I CANT ) 

153IC 
164 1 C 

15SIC IF CONVERGENCE NOT OBTAINED, ABANDON CASE. 

1SCIC 
157iC 


8 W 


isai 

159«C 

lMtC 

line 

162«C 
163iC 
164i 
166' 
!Mt 
167s 
IMt 
1691 
1781 
1711 
1721 
173IC 
174IC 
17IIC 
1721 
177i 
178i 
1781 
1881 
18H 
1821 
183s 
1841 
186i 
1881 
1871 
188IC 
189tC 
196tC 
191 tC 
192IC 
193IC 
194* 
195: 
1961 
1971 
1981 
199t 
2881 
28H 
282s 
283t 
284t 
285i 
2861 
287« 
288i 
2891 
218i 
21H 
212i 
213t 
214t 


IF < ICANT .NE. 8) GO TO 988 

8FPEAT PROCESS TO COMPUTE NEU INNER UALL TEMPERATURE. 


TNEUT2 • T(IRAX) 

LIN • 3 t IRAK 

IF (NSTEP .EQ. 3) TNEUT2 - TSAUE(IRAX) ♦ DTERP(LIN) S DT 
TERM2 • 2. I COND(IRAX)/ DX(IRAX) 

AIN — ERISIN t STEU1E 
BIN •- ( HCIN ♦ TERN2) 

CIN • AIN t TSINKU4 - HCINSTGAS - TERR2ITNEUT2 
CALL NEUT( AIN, BIN, CIN, TUI, ICANT) 

IF (ICANT .NE. 8) GO TO 918 

CORRECT MEAT FLUXES FOR NEU UALL TEMPERATURES 


IF (IDEAL .GT. 8) GO TC 418 

CALL NOLIER (HU,PL,2,TU,Z2,SS,RR,GR ) 

IF (IDEAL .GT. 8) GO TO 418 
GO TO 428 
418 MU • .24 t TU 
428 CONTINUE 

OC • MCOUT t (MREC - MU) 

OR • ERISS t STEUIE t TUBS4 

00 • -OR ♦ OC 

01 • ERISIN B STEUIE S TUIBB4 - HCIN t (TGAS -TUI) - ERISIN S 
1 STEUIE * TSINXBS4 

588 TDOT • ( TU - TUOLD )/DT 


CORPUTE RAXINUR STABLE TINE AS FUNCTION OF TME INNER AND UALL 
EQUILIBRIUM TEMPERATURES. 


IF (NSTEP .EO. 4) CALL STABLE(FIN, FOUT, CIN, GOUT, 01, TREC ) 

IF(NERROR.EO. 1 ) RETURN 

IF ( IFIRST .NE. 1) GO TO 458 

CALL STABLE (FIN, FOUT, GIN# GOUT, 01# TREC) 

IF (NERROR.EO. 1 ) RETURN 
DO 425 I • 1, IRAX 
425 T ( I ) • TSAUE(I) 

IFIRST • 8 
458 CONTINUE 
RETURN 

988 URITE(6,2888) CTIME 
GO TO 928 

918 UR1TE (6 3888) CTIME 

928 URITE(S!2881 ) COND(l), DX(1)# EM1SS, HCOUT, HREC, HU, TREC, TU, 
l NSTEP 

UR1TE(6, 3881) COND(IRAX), DX(IMAX), ERISIN, HCIN, TSINIC, TGAS, 
l NSTEP 

URITE (6,4888) (T(I># 1 • 1# IMAX) 

NERROR-1 

RETURN 

END 




II 
II 
3IC 
4iC 
SIC 
HC 
?IC 
• IC 
91 
191 
111 
III 
13iC 
14i 
lSi 
III 
17i 
lit 
19i 
I9i 
III 
III 
23i 
I4i 
2*t 
III 
17 1 
Ml 
291 
391 
3H 
3It 
33t 
34IC 
3SiC 
39iC 
37IC 
31 1 
39t 
491 
41t 
421 
431 
441 
45 1 C 
45 1 C 


471 

48t 

49t 

59t 


SUBROUTINE CRSFLU(CFFLG,ELt ,ELT,EL,ELNBDA, DSUBO, CORNR, UE, PE, 
1 RHOE , PU , ALPHA , UDOT , UU , XNACHU , WX ) 

CONSTANT WIDTH RECTANGLE - IDEAL GAS 
CONSTANT WIDTH RECTANGLE - REAL GAS 
SHARP EDGE DELTA - IDEAL GAS 
SHARP EDGE DELTA - REAL GAS 

IF< ALPHA.LE.9. ) CO TO 98 
NCFFLC*CFFLG*.99l 
ALPHR- ALPHAS. 91 74S33 
IFCNCFFLG-2 ) 1,2,19 

1 WBAR-C .7450. 1 4IC0RNR/DSUB0 >f . 447BTANC ALPHR ) 

GO TO 5 

2 IF (UDOT .LE ) UDOT-5.31 
ENN-XHACHUISINC ALPHR ) 

IF(ENN.LT. I • ) CO TO 3 
DUDXD • SORT ( 2 . t ( PE -PU ) /RHOE ) 

GO TO 4 

3 DUDXD- 3. BUUtSINC ALPHR ) 

4 UBAR-2. SUDOTSDUDXD/UE 
IF<NCFFLC.GT.2> GO TO 12 

5 UX-DSUBO/WBAR 
UV-EL/UX 

IFCWV.LT. l.E-5) GO TO 98 
IF CWV.GT .6. ) GO TO 5 
ELL-9.5K 1 . -EXPC -2. BUV ) )SUX 
ELT- ( 1 . -EXP( -UV ) )fUX 
GO TO 99 
5 ELL-9. SSWX 
ELT-WX 
GO TO 99 

SHARP EDGE DELTA (SWEEP ANGLE LANBDA ) 


19 IF (NCFFLG-3 )1 1 ,11,2 

11 BETA-1. 4. 333tTAN( ALPHR )BTAN<ELHBDA1 .9174533 ) 
GO TO 13 

12 BETA- 1 . ♦WBAR1TAN( ELRBDAB .9174533 > 

13 ELL-EL/BETA 
ELT.2.BEL/CBETAM. ) 

GO TO 99 


98 ELL-EL 
ELT-EL 

99 RETURN 
END 




II SUBROUTINE DETRALIRN.EL.PH1.ENCL.ENCT) 

ll DIMENSION A(S> 

3i CONAON/FLUFLO/TINF .PINF, RMOINF . AINF.RMOE.PE, TE.HE, XNUC. UE.RMOO.PO. 

41 I TO. HO. XAUO.RHOU.Ty.HU.XHUU.RHOS.TS.HS.XNUS. RMOR.TR, MR. XNUR.U.XAC 

Si 2, REl, HINT. XNUINF. PR, PU. DU. TU.XNUU. XN4CHU.WU.hu, ALPHA. 

Si 3RHOSTL.HSTL.TSTL,XRUSU.RHOSTT.HSTT,TSTT.XRUSTT.HRECL.HRECT.S 

7l t.CINF.CAAAU, CANAS, CAHAE.GANAO.GU.CAHAC.PRL.PRT 

Si CONNOH/ERPrlC/NCRROR 

01 DATA (A(J).J-l.fi)/ I .03754, .0043770. - .6107E-4. .904SI, .01107, 

101 1 -.04MSBE-4/ 

111 REL'RHOESUEtEL/XNUC 

111 1 THETA* l . S70796-PHH .01714S3 

131 IF (PHI .CT .89. ) FTR* 1 . 

141 IF (PHI .07.00. 1 CO TO 2 

IS t C*t ./< 1 . 4tXAACHUtV2 I 

IS • D*(l.-C>a(THETA<«2-THCTAtSIN(4*THETA)/2.«(l.-C0S(4.tTHETA) >/0) 

171 144.tC«(THETAtta-THETA*SIN(2.tTHETA)4(l.-C0S(2.«THETA)»/2. ) 

111 FTR*2.tTHETA*SINi THETA M( ( 1 .-C >*COS(THETA )I*2*C >/$ORT(D ) 

101 2 F3* ( HO-HW )/ ( HRECl-HU > 

101 CALL FAVRID(RN.ENFI 

111 PRL* .71 

III PRT*.71 

23 i ENCL*ENF>FTRtF3 

I4i 10 sun*0. 

Ml IFIPMI.LT. 25. » n*l 

2Si IF (PHI . GC .25 . I N*4 

271 N*H*2 

Ml K*0 

291 DO 11 I*A.N 

301 SUN*SUA*A( I >tPHItt(K I 

3H 11 K*K*1 

32 i F*SUA 

33i RR*RH0E'2.37?E-3 

341 CALL BINTRP(TE,RR,1,HD) 

3Si IF (HERROR.HE .0 I RETURN 

301 ENCT* . 933*REL*t . 8*< 1 .037SF )* ( 1 . ♦ .S8*HD/M0 >«XAUE/(EUPRT*» .607 I 

371 S*PRL».007 

3Si 90 RETURN 
301 END 
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693 

Ill 


III 


131 


Hi 

199 

111 

192 
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992 

lit 
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tit 


III 


23i 

994 

14 • 


til 

995 

211 


27i 
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2li 


3li 


31i 

996 
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33» 
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36» 

111! 

37i 


3li 


32 1 


4|« 

913 

411 


42t 

912 

43i 


441 
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471 

911 

481 
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59i 
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SUBROUTINE DINT(XX.XT,Vl,VTl,2Z.ZT,NL.Nl,V2.VT*.V3,VT3.Il,VE,VT4> 
DOUBLE INTERPOLATION SUBROUTINE. IF IL*1. THE SECOND INDEPENDENT 
VARIABLE IS NOT CONSTANT WITH THE FIRST. (THE RANGE OF THE SECOND 
IS NOT THE SANE AT EACH VALUE OF THE FIRST. * 

OINENSION ZT(NL >, XT(NL,NL >, VTI < NL, NL >. VT2( NL.NL >»VT3(NL,NL I 

DINENSION VT4INL.NL I 

IF(2Z-ZT(1 >>BR2,M1,B#1 

IF ( IL.EO. 1 ) CO TO BB3 

IF ( XX . LT.XTI 1, 1 > ICO TO BA2 

DO DM 1*1. NL 

L*l 

11*1-1 

IF ( 22 -2T ( I > I9A2.8B4.8M 


RETURN 

RATIP*<ZZ-ZT(LL>>/(2T(L)-ZT(Ll)> 

IF<tL>llll. till. 1113 
DO M4 J-l.NL 
LN*J 
LLN-J-1 

IF(XX-XT(J,L> IMS, 9B6.M4 

CONTINUE 

GO TO M2 

RATIO* I XX-XTI ILN.L > >✓( XT( LN.L >-XT( L.~'.L > > 

VI *YT1 (ILN.L )*RATIO«( VTI ( LN. L >-VTl ( LLN, L > > 
V2*VT2( LLN.L >*RAT 10* (VT2( LN.L )-VT2( LLN, L 1 1 
V3*VT3( LLN.L >*RATIO*< VT3( LN.L >-VT3( LLN, L 1 ) 
VE*VT4( LLN.L >*RAT10*(VT4( LN.L )-VT4( LLN, L > 1 
RETURN 
Vl*VTl(LN,L> 

Y2*VT2( LN,L ) 

V3*VT3(LN.L> 

VE-VT4UN.L) 

RETURN 

DO 913 J-l.NL 
LN* J 
LLN*J-1 

IF (XX-XT(J, 0)91 1.912. 913 

CONTINUE 

CO TO M2 

VI *VT1 < LN, LL >-RATIP»l VTI (LN, LL >-VTl ( LN.L * ) 
V2*VT2< LN.LL >-RATIP*( YT2( LN.LL >-VT2( LN.L ) > 
V3*VT3( LN, LL >-RATIP*( VT3( LN, LL )-VT3(LN,L > ) 
VE*VT4( LN.LL >-RATIP*( VT4( LN.LL >-VT4( LN.L > > 
RETURN 

RATIO* CXX-XTI LLN.L > >/( XT( LN, L l-XTI LLN.L > > 
Y4-VTI ( LLN.LL > ♦RATIO* (VTI (LN.LL >-YTl ( LLN.LL > > 
VS*VT1 ( LLN.L >♦ RAT 10* (VTI ( LN, L >-VTl ( LLN.L > ) 
Y1*Y4-RAT1P*(V4-VS> 

V4*VTe( LLN.LL >*RATIO*( YT2( LN.LL >-VTZ( LLN.LL > > 
YS* VT2( LLN.L >*RATIO*( YT2( LN.L l-YT2(LLN,L > > 
V2*Y4-RATIPt( Y4-YS ) 

Y4*YT3( LLN.LL >«RATIO* < VT3(LN. LL >-YT3( LLN, LL > > 
Y5*VT3(LLN.L>«RATIO*(VT3(LN.L>-YT3( LLN.L* > 
V3*V4-RATlPt(Y4-VS I 

V4*YT4( LLN.LL >♦RAT10*( YT4(LN.LL >-YT4( LLN.LL > > 
YS*YT4(LLN,L >*RATI0* < VT4( LN.L >-VT4( LLN, L > > 
YE*V4-RATIP*( Y4-VS > 

RETURN 

DO 91B J* 1 ,NL 

IF(XX.LT.XT<1,L> >CO TO BR2 

LN* J 

LLN-J-1 

IF(XX-XT(J.L> >916,917, 918 
CONTINUE 
GO TO M2 
VS-YTKLN.L > 

V7*YT2(lN.l > 

Y9*VT3( LN, L > 

Y#*VT4(LN,L > 

CO TO 1114 

RATIO*(XX-XT( LLN.L > >'(XT(LN, L >-XT ( LLN.L > > 

VS* VTI (LLN.L >«RATIO*< VTI (LN, L >-VTl ( LLN.L > > 
V7*VT2( LLN. L MRATIO* ( YT2( LN. L >-YT2( LLN, L > > 
Y9*VT3(LLN,L >tRATIO*(VT3( LN.L >-VT3( LLN.L > > 
VB*YT4( LLN. L>tRATIO*(VT4(LN,L>-VT4( LLN.L >> 


Sfc P40 

M>Ii QUai. 


H-24 


8838855828883885382888 


I 1114 DO 021 J-l.AL 
LA-J 
LLAO-1 

IF (XX-XT( J. LL ) >919,929,921 
921 COATINUC 
GO TO 892 
929 V4-YTKLA.LU 
V9*VT2(UI»LL I 
Vt-VT3<LA,LL> 

YA-VT4<LA,LL> 

GO TO 928 

919 RATIO- <XX-XT< LLP, LL>)/»XTi LA, LL>-XT< LLP, LL> > 
Y4-VT1 CLLP,LL >*RATIOiC YT1 ( IP, LL >-YTl ( LLP, LL I > 
Y9-VT2C LLP,LL )*RATIOi( VT2C LP, LL )-YT2(LLP,LL ) I 
Y2-VT3<LLP,LL J^RATIOt C VT3( LP, LL >-YT3< LLP, LL > > 
YA-YT4CLLP,LL>*RATI01CVT4<LP,LU-YT4<LtP,LL>) 
92S Yt-Y4-RATIPt<Y4VS> 

Y2-Y9-RATIPt<YG-Y7> 

V3-Y*-RATlPt(Yi-Y9) 

YC • YA- RAT IPt C YA- Y9 ) 

RCTURH 

CND 


8-2 f) 







il SUBROUTINE DINTHXX.XT.Yl.VTl.ZZ.ZT.NL.NL) 

2«C DOUBLE INTERPOLATION SUBROUTINE. 

3t DINENSION ZT(NL).XT<NL,NL> YTMNL.NL) 

41 1F(Z2-ZT<1> >892.891,891 

S« 801 IF (XX-XT(l.D) 885.893,803 
6t 885 Yl-12.7 
7» RETURN 

•t e»? V1-.2222E+39 
9< RETURN 

181 L83 DO 888 1*1. NL 

111 L-I 

12i LL-I-1 

13i IF(ZZ-ZT(I I >982.884.889 

141 808 CONTINUE 

1S> GO TO 382 

16i 992 TTIP- <ZZ-ZT(LL ) >/(ZTCL >-ZT(LL > ) 

1?« JC TO 4 113 

181 894 r 984 J-l.Nl 

19i wl-J 

29t LLN-J-i 

211 IF< XX-XTC J.L > >985,996.984 

22 l 984 CONTINUE 
23i GO TO 892 

241 98S RATIC*(XX-XT(LLn,L>>/(XT(LN.L)-XT(LLn.L>> 

25 < V1*YT1 (LLN, L )-»RATIO*< YT1 ( LN.L >-YTl (LLH.L ) > 

26t RETURN 

27« 986 V1-YT1(LN,L> 

281 RETURN 

291 1113 IF <XX-XT< 1 , L > ) 89S.1115.111S 
38i 1115 DO 918 J-l.NL 
3l« LN-J 

32i LLH-J-1 

33* IF<XX-XT<J.L> >916.917,918 

34t 918 CONTINUE 

3S« CO TO 892 

36« 917 YS-YTl(LN.L) 

371 GO TO 1114 

381 916 RATIO*(XX-XT( LLN, L > )/(XT(LN,L )-XT(LLN,L > > 

39i Y5-YT1 (LLN, L >*RATIO*< YTl (LN, L )-VTl ( LLN, L ) > 

49i 1114 DO 921 J-l.HL 
41 i LN-J 

42l LLH-J-1 

43i IF (XX-XTC J, LL ) >919,929,921 

441 921 CONTINUE 

4S> GO TO 892 

46l 928 V4-VT1 (LN.LL > 

47 t GO TO 925 

481 919 RATIO- <XX-XT< LLN.LL ))/< XT < LN, LL > -XT C LLN.LL ) ) 

491 V4-VT1 ( LLN.LL HRATIOTC YT1 ( LN.LL >-YTl ( LLN.LL > > 

S8i 925 Y1-V4-RATIR*(>^-V5) 

SI i RETURN 

S2t END 


8-26 


SUBROUTINE DOUNID< ANGLE ) 
REAL NUSTRfl, HFRSTH 
REAL IIDNSTfl 
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SUBROUTINE DRIUEL ( DTEPIP, G, ZPLUS, ZPIINUS, QO, 01, T, NSTEP, IflAX) 
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SUBROUTINE FINALT ( TSAUE, DTEUP, IPIAX, NCOL 
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28 
21 


SUBROUTINE FLOUCFF, ALFA, BETA ) 

REAL MU, ND , ML , NE 
DIMENSION FF ( 9 ) , ALFA ( 9 ) , NF ( 9 ) 

CONNON/F LUF LD/ T I NF , P I NF , RHO I NF , A I NF , RHOE , PE , TE , HE , KNUE , UE , RHOO , PO , 
1 TO , HO , KNUO , RHOU , TU , HU , XNUU , PHOS , T S , HS , XNUS . RHOR , TR , MR , XNUR , U , XNE 
2, REL , HINF , XNU I NF , PR, PU, DU , TU , XNUU , XNACHU , UU , HU , ALPHA, 

3RH0STL , HSTL , TSTL , XNUSTL , RMOSTT , HSTT , TSTT , XNUSTT , HRECL , HRECT , S 
t, G INF, GAN AU, GAMA S, GAMAE, G AMAO. GU, GANAG, PRl,PRT 
CONNON/DNSTRN/TD, PD, DD, ND, VD, i>D, HD , AD, GANAD 
CONNON/ FRSTN/TFRSTN, PFRSTN, DFRSTN, MU. UFRSTN, SU, HFRSTN, AU 
COMMON /TOTAL/' HTOTAL , PTOTAL , TTOTAL , RMOTOT , CTOTAL 
COMMON/CVL/ H$L,PSL,RHOSL, XNUSL, TSL,USL 
COMMON/STAC/PSTAG, RHOSTG, TSTG, XMUSTG 
COMMON/FLAG/ IDEAL 
COMMON/ERRFLG/NERROR 
COMMOM/HFLAG/ NHFLAG 
DATA PI/3. 1415926S4/ 

IDEAL-8 
GR-32.2BS3.3S 
DO 8 K-1,9 
8 NF(K)-FF(K >4.881 
NINF-.24BTINF 

CALL HANSEN(XMUINF,PINF,TINF) 

IC-1 

GINF-1.4 

TU-TINF 

PU-PINF 

DU -RHO INF 

MU-U/AINF 

UU-V 

HU-HINF 

AU-AINF 

GAMAU-GINF 

XNUU -XMU INF 

PFRSTM-PU 

TFRSTM-TU 

DFRSTM-DU 

UFRSTM-UU 

HFRSTN-HU 

1 IF(MU . LT . 1 . > GO TO 28 
IFCNF (K )-36 >2,4,6 

2 CALL PCSUCMU, ALFAOC ), 1 , OFT, BETA ) 

IF(NERR0R.NE.8> RETURN 

IF (OFT- 1 . ) 18,28,5 

4 CALL PCSU(MU,ALFA(IC>,2,0FT,BETA> 

IF (NERROR.NE .8 > RETURN 

I F ( OF T - 1 . >18,28,5 

5 BETA-98. 

GO TO 10 

6 BETA- ALFA (1C) 

10 CALL DUNSTM( BETA* PI / 180) 

1F( IDEAL-1 >30,100,30 

IF MACH NO IS LESS THAN ONE. USE MOD. NEWTONIAN AND LAST SLOPE 


20 IF ( K . LT .2 ) GO TO 160 
HO-HTOTAL 

PO-PTOTAL 

RHOO -RMOTOT 

TO-TTOTAL 

GAMAO-GTOTAL 

CALL MANSEN< XNUO,PO, TO ) 

22 I F ( NF ( 1C 42 > • EG . 29 ) GO TO 26 
BETAP-ALFACIC41 )tPI/180 
IC-IC42 

27 IFOC.CT.7) GO TO 25 
IF< ALFA(K )-0. >25.25,22 
26 BETAP-(ALFA(IC4l )-ALFA0C42 ) )tPI/188 
IFCBETAP.LT. 0. > BETAP-0. 

IC-IC43 
CO TO 27 

25 PE-PU4(PO-PU )f SINCBETAP >tt2 

CALL MOLIER(HE,PE, l,TE,ZE,SD,RHOE,GAMAE > 
IF ( IDEAL .GT. 0 ) GO TO 168 
IF(HE.GT.MO) GO TO 23 
IF (BETAP .CE. PI/2 ) GO TO 23 
UE2-CH0-HE >150183. 



H- V 


Til UE- SORT CUE 2) 

791 AE*S0RTCGAnAESCR>2E<TE> 

Ml XAE-UE/AE 

lit CALL HANSEN (XIHJE, PE. TE> 

III CO TO 24 

131 23 HE’HO 

141 PE*PO 

lit TE-TO 

Ml RHOE'RHOO 

I7t XNUE-XHUO 

III GANAE'GAHAO 

Ml UE*9. 

Ml XNE-9. 

91 l 24 HSL*HE 
921 PSL-PE 

931 TSL-TE 

941 XHUSL*XNUE 

951 OSL*UE 

Oil RHOSL-RHOE 

97l HD*HU 

Oil PD*PU 

991 TD-TU 

INI DD-DU 

1911 ND-NU 

1921 UD-UU 

1931 SD-SU 

1941 AD'AU 

1951 CANAD'CANAU 

1961 CARASL'CANAE 

1971 PSTAG-PO 

19li RHOSTC-RMOO 

1991 TSTC*TO 

ii9i xnustg-xnuo 

1111 GO TO 52 

U2IC 

1 1 3 i C DETERMINE LOCAL FLOU CONDITIONS 

114IC 

HSl 39 IF(NF (K ) .GT . 38 ) GO TO 43 
llll XFCNFCK+l i.EO.14) GO TO 38 

1171 IF (NFOC+l )-16 > 31,32.36 

llll 39 IFCNFCICM ) -19) 46,34,37 
1191 31 XTA8LE-3 

1891 GO TO 33 

18H 32 ITA9LE«4 

1221 33 CALL PCSUCNU.ALFACIC+l ), ITAILE.OFT.PC ) 

1131 IF(NERR0R.NE.9) RETURN 

1241 8ETAP«ALFA<KM>«P1'189 

18Si IF(0FT-1.)35. 29,34 

1261 3S IFCPC.LT. 9.) PC-9. 

1271 PL-PU*PC8DU8UU*UU/2. 

1281 GO TO 49 

129> 34 BETAP-ALFACKM )*PI/189 

1391 PL-PU+CPTOTAL-PU HSINC1ETAP »82 

1311 CO TO 49 

1321 37 PSAU-PTOTAL 

1331 TSU-TTOTAL 

134* RSU-RHOTOT 

1351 SSU-SD 

1391 CSU-GTOTAL 

137i HDSU-HD 

1311 PDSU-PD 

139t TDSU-TD 

1491 DDSU-DD 

14H CDSU-CANAD 

1421 XNDSV-ND 

1431 UDSU-UD 

1441 ANGLE-PI '2 

1451 CALL DUNSTNC ANGLE) 

1461 PSTAG-PTOTAL 

1471 PTOTAL-PSAV 

1491 RHOTOT -RSU 

1491 TTOTAL-TSU 

1591 SD-SSU 

1511 GTOTAL-CSU 

1621 HD-HDSU 

1631 PD-PDSU 

1541 TD-TDSU 
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155i M>*DDSU 

154i GAHAD-CDSU 

1571 HD-XHDSW 

15tl OD-UDSU 

ISO' BETAP-ALFA(KM)*PI/184 

lMt PL-PU*<PSTAG -PU>tSIN<BETAP>*»2 

1411 CO TO 44 

1621 38 PC-ALFA<KM> 

163i BETAP-4. 

1641 PL-PU+PC*DU*UUtWU/2. 

1651 44 CALL AOLIERCML.PL, l.TL.ZL.SD.DL. CANAL) 

1661 1F( I DEAL. GT. 4) GO TO 144 

167« ZTL-ZLBTL 

1681 AL'SQRT 1 GAAALtCRtZTL ) 

1681 IF<ML.CT.MTOTAL> CO TO 41 

1741 UL2-<HT0TAL-HL)«S4143. 

17H UL*S0RT(UL2> 

1721 RL*UL/AL 

1731 CO TO 42 

1741 41 HL«HTOTAL 

1781 PL-PTOTAL 

1741 TL*TTOTAL 

1771 DL-RMOTOT 

1781 OARAL'CTOTAL 

1781 Ul-UU*C09(BETAP)M2 

1841 IF (8ETAP .EO. 4.) WL-4. 

18H RL»UL/AL 

188* 00 TO 42 

1S3IC 

1841C CALCULATE 04LI0UE SHOCK PRESSURE 

186 1 C 

1861 46 WLT*WU«C0S(BETASPI/184> 

1871 IF<ALFACK*1 KCT.BETA) CO TO 43 

1681 ANGL*PI/184*(BETA-ALFA<K*in 

188' UL«ULT/COS(ANCL> 

1841 HL*HT0TAL-UL«UL/54143. 

1811 CALL N0LIER(HL,PL.3,TL.ZL.SD.DL.GAI1AL> 

182' AL*SORT(CAnALtCRtZL>TL ) 

183' PIL-VL/AL 

1841 BETAP-ALFA<KM XPI/184 

185' 42 IF(NF(K+2).NE.29) GO TO 44 

186' CALL PnEXPN(ALFA<K»2>.HT0TAL,HL,AL.SD. PL. TL.UL.DL.NL. CANAL) 

1871 IF < IDEAL.CT.4 ) CO TO 144 

188' I F ( NERROR . GT . 4 ) CO TO 999 

199' BETAP*(ALFA(K<1 l-ALFACK+2 > XPl/184 

244) K*KO 

241' CO TO 45 

242' 44 K>K«2 

243' 45 IF (K .GT .9 ) GO TO 54 

2441 IF(NL. LE .4. > CO TO 54 

245' IF (NF (K > .LE . 4 ) GO TO 54 

246< PU-PL 

247i TU«TL 

248' DU*DL 

249' HU-HL 

214' AU*AL 

211' OU*UL 

212' NU-NL 

213' CAHAU-CANAL 

214' PFRSTR-PU 

21S' TFRSTA-TU 

216' Df RSTN-DU 

217' UFRSTN-OU 

218' HFRSTN-HU 

2191 GO TO 1 

224 ' C 

221 'C SWEPT CYLINDER STAC LIHE PROPS 

222 1 C 

223) 43 CALL HOLIER'HSL.PSL. 3.TSL. ZSL.SD.RHOSL.GAHASL > 

2241 IF ( IDEAL.CT.4 > CO TO 144 

225' CALL HANSEN < XNUSL.PSL.TSL ) 

226' USL-OU*COS(BETA*PI/184) 

227i PE-PSL 

228* HE-HSL 

229' TE*TSL 

234' RHOE*RHOSL 

2311 XNUE*XNUSL 
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mi 

is: 

•XI 

•Mi 

•371 

•Ml 

•Mi 

Mi 

•41 1 
2481 
8431 
8441 
84Si 
84<i 
847t 
848t 
848i 
8S8i 
86 l« 
852 s 
8531 
854t 
85St 
85€i 
8571 
8581 
8591 
258! 
25H 
252i 
253i 
254t 
255i 
2551 
2571C 
258 1 C 
259 1 C 
278i 
271« 
272i 
273i 
8741 
275i 
2751 
277i 
278! 
279i 
288: 
281! 
882t 
283i 
284i 
28Si 
285i 
287i 
288 1 C 
289 !C 
898 1 C 
89H 
298« 
293i 
294i 
295 1 C 
295iC 
297 1 C 
2981 
8991 
388i 
MU 
3821 
M3i 
M4t 
3851 
3861 
M7« 
3881 
3891 


AC*$QRT(CAIIASLiCRtTSL ) 

UC-VSL 

XHE-UE/AE 

GAHAE-GANASL 

PSAU-PTOTAL 

TSU-TTOTAL 

RSU-RMOTOT 

GSV-CTOTAL 

MDSU-ND 

PDSV-PD 

TDSU-TO 

DDSU-DD 

GDSU-GANAD 

XNDSU-ND 

UDSU*UD 

ANGLE-PI/2 

CALL DUNSTH( ANGLE) 

PSTAG-PTOTAL 

CALL NOL1ER(NTOTAL,PSTAC,8,TSTG,ZSTC,SSTG,RNOSTG.GARSTG) 
1F( IDEAL • NE • 9 ) GO TO 188 
CALL HANSEN(XNUSTC,PSTAG,TSTC) 

HSL-HE 

PTOTAL-PSAU 

RHOTOT-RSU 

TTOTAL-TSU 

GTOTAL-GSU 

HD-HDSU 

PD-PDSU 

TD-TDSV 

DD-DDSU 

GAMD-GDSU 

ND-XHDSU 

UD-UDSU 

BETAP-9ETAfPl/188 
GO TO 51 

DEFINE BOUNDARY LAVER EDGE PROPERTIES 

58 PE -PL 
TE-TL 
RHOE-DL 
HE -ML 
UE-UL 
AE-AL 
RE-NL 
XNE-HE 
GANAE- CANAL 

CALL HANSEN ( XNUE # PE , TE ) 

51 PO-PTOTAL 
HO-HTOTAL 
RHOO-RHOTOT 
TO-TTOTAL 
GANAO-GTOTAL 

CALL HANSEN <XNUO,PO. TO) 

52 CALL HANSEN <XNUU,PU.TU> 

XHACHU-NU 

DEFINE WALL PROPERTIES 

CALL NOLIER<HU # PE,2,TU,ZU.SU,RHOU,GU> 

IF ( IDEAL . GT . 8 ) GO TO 188 
CALL HANSEN (XNUU, PE ,TU) 

ALPHA- IE TAP /P I 1 1 88 

ECKERT REFERENCE ENTHALPY PROPERTIES - REAL GAS 

HRECL*HE«8.85SUEfUE/58183. 

HRECT-HE#9.88tUEtUE/S8193. 

HSTL-8. 51 (HE^HU H9* 221 (HRECL-HE ) 

HSTT *8. 51 (NE ♦HU )^8.22f (HRECT-HE ) 

CALL NOLIER(H5TL,PC, 8, TSTL, ZS # SS # RHOSTL, GANAS ) 

CALL NOL I ER ( HSTT , PE , 8 , TSTT , ZS . SS , RHOSTT . G ANAS ) 

IF ( IDEAL. EG. 1 ) GO TO 188 
CALL HANSEN (XNUSTL, PE ,TSTL) 

CALL HANSEN (XHUSTT, PE, TSTT) 

IF(NHFLAG.EO.S) GO TO 78 
IF(NHFLAC.E0.7) GO TO 78 
GO TO 999 
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319*C 


311 tC 
312*C 

DEFINE RMO R , HU R PROPERTIES 

3131 7# 

CALL HOLIER<MO.PE#9#TSP#ZSP.SSP, RHOSP . GAHASP ) 

3141 

IF ( 1DEAL.GT .9 ) GO TO 199 

315i 

CALL HANSEN (XHUSP. PE # TSP ) 

31Si 

RHSP-RHOSPtXHUSP 

317i 

RHE-RHOEIXHUE 

318t 

RHEE -RffSPt ( 1 . 85- . 858 ( RHSP/RHE ) ) 

319i 

RHU-RHOUtXHUU 

32ti 
321 *C 

RORHUR-RHEE* C 1 . 6- . 58 C RHEE/RHU ) ) 

322 * C 
323 tC 

ITERATE FOR RHO R 

324i 

1*1 

32Si 

RHP-RORHUR/( 1 .E-13SPE ) 

325t 

IF ( RHP-4 . )69» 69.61 

327t 59 

TRGESS-2.21E4tRHPtt< -2. 47 ) 

328* 

GO TO 62 

329* 51 

TRGESS • 2 . E7tRHPtt (-7.39) 

339* 52 

IF ( TRGESS. GT . 8999.) TRGESS- 8999. 

331* 

CALL HOLIER ( HR , PE . 2 . TRGESS .ZRG.SRG, RHORG , GAHAG 

332« 

IF ( 1DEAL.GT.9 ) GO TO 199 

333* 

CALL HANSEN (XHURC. PE . TRGESS) 

334* 

RAGE S S • RHORG t XNURG 

335* 

ERR-ABSt ( RNGESS-RORNUR )/RORNUR ) 

336* 

IFCERR.LT. .91 ) GO TO 65 

337* 

RESID-RORHUR-RNGESS 

338* 

TRGP2-TRGESSI1 .95 

339* 

TRGN2-TRGESS19.95 

349* 

CALL H0LIER(HR,PE.2, TRGP2.ZRG.SRG.RH0RP2.GAHAG 

341* 

IF < IDEAL. GT. 9) GO TO 199 

342* 

C.^LL HOL I ER < HR . PE . 2 . TRGN2 . ZRG . SRC . RH0RH2 , GAHAG 

343* 

P( 1DEAL.GT. 9 ) GO TO 199 

344* 

TERIOR- ( RHORP2-RNORH2 )/ C . ltTRGESS ) 

345* 

CALL HANSEN (XHURP2. PE. TRGP2) 

346* 

CALL HANSEN (XHURH2. PE. TRGN2) 

347* 

DERIUH- ( XHURP2-XNURH2 )/( . ltTRGESS ) 


DERIU— XHURGSDERIUR-RMORCtDERIUH 

RATE-RCSID/DCRIO 

TRGESS- TRGESS-RATE 

351* 

I-IM 

358* 

2F(I.LT. 59) GO TO 52 

353* 65 

RHOR-RHORG 

354* 

XHUR-XHURG 

355* 

TR-TRCESS 

356* 

357*C 

GO TO 999 

358 iC 
359*C 

IDEAL GAS FLOU PROPERTIES 

389* 199 

K-l 

361* 

TU-TINF 

368* 

PU-PINF 

363* 

DU-RHOINF 

364* 

HU-U/AINF 

365* 

UU-U 

366* 

HU-HINF 

367* 

AU-AINF 

368* 

GANAU-GINF 

369* 

PFRSTH-PU 

379* 

TFRSTH-TU 

371* 

DFRSTH-DU 

372* 

JFRSTH-UU 

373* 

MFRSTH-MU 

374* 191 

IF (HU.LT . 1 • ) GO TO 159 

375* 

lF(NF(K)-36 M92.193.196 

376* 192 

I TABLE -1 

377* 

GO TO 194 

378* 193 

ITABLE-2 

379* 194 

CALL PCSU < HU . ALF A U > , I TABLE . OF T . BETA ) 

389* 

IF ( NERR0R.NE.9) RETURN 

381* 

IF (OFT-1 . >197.29. 195 

382* 195 

BETA-99. 

383* 

GO TO 197 

384* 195 

BETA-ALFA(K) 

3«Si 1*7 CALL DOUNIDI BCTAtP 1/ 18t ) 
MfiiC 
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397 IC DCTERNINE LOCAL FLOW CONDITIONS 
3ttiC 

399! IF(NF(K ).CT .31 ) 60 TO 111 

399! IF(NF (KM ) .EO. 14 ) 60 TO 121 

3911 IF<NF(KM >-l!> 198.199,123 

392t 123 IF(NF(KM )-18 ) 124.111.12S 

3931 199 ITABLE-3 

3941 GO TO 119 

3911 199 ITABLE-4 

3911 119 CALL PCSU(NU, ALFACKM ). I TABLE. OFT. PC ) 

3971 IF ( NERROP . NE . 9 ) RETURN 

399t BETAP-ALFAOCM )BPI/199 

3991 IFtOFT-l. )U2,29.lll 

499i 111 BETAP-ALFA(K4l)tPI/lB9 

4911 PL*PU4(PT0TAL-PU)tSIN(BETAP)Xt2 

4921 GO TO 113 

493t 112 IFCPC.LT.l.) PC-9. 

494 I Pt -PU4PC IDUXUUtUU/2 . 

4911 GO TO 113 

4911 181 PSAU-PTOTAL 

4971 TSU-TTOTAL 

4991 RSU-RMOTOT 

4991 ANGLE -PI/R 

4191 CALL OOUNID( ANGLE) 

41 1 « PSTAV PTOTAL 

4181 PTOTAL- PS AU 

4131 RHOTOT-RSU 

4141 TTOTAL-TSU 

41H BETAP-ALFACKM )tPI/lB9 

41H PL • PU4 C PST AG -PU ) t S I N < BE T AP ) 1 12 

4171 GO TO 113 

418! 121 PC-ALFACKM) 

4191 BETAP-9. 

4891 PL-PU4PCtDUSUUtUU/2. 

4811 113 IF (PL .CT. PTOTAL) PL-PTOTAL 

4821 NL-SORT(S.t( (PTOTAL/PL )St(2./7. >-l . ) ) 

4831 HT0TAL-9.24STINF+USU/S9193. 

424i TT-MT0TAL/.24 

425 1 TL-TTtl PL/PTOTAL )BI . 281 

4811 DL -PL/(GRtTL) 

427i AL-49.92tS0RT(TL) 

4881 HL-9.24ITL 

429f CANAL- 1.4 

4391 IF(HL.GT.HTOTAL) GO TO 114 

431 « UL-S0RT( (HTOTAL-HL ) 159 193. ) 

4321 GO TO US 

433! 114 HL-HTOTAL 

4341 PL-PTOTAL 

4311 UL-UUtCOSCBETAP )1S2 

431! IF (BETAP .EO. 9.) UL-9. 

437! NL-UL/AL 

438i GO TO US 

439!C 

449 iC CALCULATE OBLIQUE SHOCK PRESSURE 

441 IC 

4421 124 VIT -UUICOSC BETA1PI/1B9 ) 

4431 IF (ALFACKM ).GT .BETA ) GO TO 116 

4441 ANGL-PI '1891 ( BETA-ALFA(K + 1 ) ) 

44H UL-ULT/COS(ANGL) 

441i MTOTAL- .24BTINF ♦UtU/59193. 

4471 ML-HT0TAL-UHUL/S9193. 

441 1 IF (HI «LT .9. ) GO TO 111 

4491 TL-HL/.24 

459* PL-PD1(TL/TD)BB3.S 

411* DL-PL/(GRtTL) 

452* AL • 49 . 92BS0RT ( Tl ) 

453* NL-UL/AL 

4541 CANAL- 1 . 4 

455i BETAP-ALFACKM )XP1/189 

451* US IF(NF(K*2).NE.29) GO TO 117 

417* CALL PNID(ALFACK42),ML,AL.PL.TL,UL.DL,Nl) 

458* BETAP- ( ALFA (KM >-ALFA(K*2 ) )tPl/189 

459* K-KO 

419* GO TO 118 

4111 117 K-K+2 

412* 118 IFCK .GT.9 ) GO TO 129 

413! IFCNL.LE.9.) GO TO 129 


RI 
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4 

i 

4 


IF(Hr<n.tE.9> 00 TO 129 

FU-Fl 

TU«Tl 

DU-DL 

4ttt 

HU •Ml 

4*91 

AU-Al 

479l 

uu-ui 

471 • 

HU-HL 

4?tl 

GAMU-GAHAl 

4731 

PFRSTH*PU 

4741 

TFRSTH-TU 

4791 

OFHSTH-DU 

47ii 

UFRSTH-UU 

4771 

HF28TH-HU 

4791 

60 TO 191 

479IC 
499 iC 

MCH NUH8ER LESS THAN 0 H£( IDEAL GAS 

491 iC 
4911 1S9 

IFCC.LT.2) CO TO 189 

493* 

HO-HTOTAl 

4941 

PO-PTOTAL 

49ii 

RHOO-RMOTOT 

4991 

TO-TTOTAl 

497i 

CAMO-CTOTAl 

499* 

CO TO 162 

4991 169 

Po*PUf<i.*.2snusmj)tt3.s 

4991 

HO-MUfUUtUU/59183. 

49H 

TO-MO/.24 

4991 

GAMO* 1.4 

4931 162 

CALI HANSEN (XHUO,FO, TO ) 

4941 16S 

IF (NF (K^2 ) .EQ.29 ) GO TO 166 

499i 

•CTAP-ALFAOCU )IPI/188 

4991 

K-K42 

497i 167 

IFCK.CT.7) GO TO 168 

4991 

IF ( ALFACK )-9. >168,168,169 

4991 166 

IE TAP- ( ALFA(KM >-ALFA<K42 ) UP1/189 

9991 

IF < 8ETAP. LT .9. ) IE TAP- I. 

99H 

K-KO 

9981 

GO TO 167 

9931 169 

PE-PUMPO-PU MSINC IETAP )SS2 

9941 

T0-M0/.24 

999i 

GAHAO* 1 .4 

996t 

TE-T0S(PE/P0)tt.286 

9971 

RMOC*PE/<GRITE> 

99«> 

RHOO *P0/ ( 32 . 2IS3 . 3SIT0 ) 

S#9i 

HE-9.24STE 

9191 

GAME-1.4 

9111 

IF(HE.CT.HO) GO TO 163 

S12i 

IF (IETAP .CE. PI/2) GO TO 163 

9131 

UC*SQRT ( (MO-HE >159193. ) 

9141 

AC • 49 . 92SS0RT ( TE ) 

919i 

XNE-UE/AC 

9161 

CALL NANSEN (XHUE, PE. TE ) 

9171 

GO TO 164 

S19I 163 

PE -PO 

919t 

ME -NO 

929t 

TE*TO 

S2H 

RMOC-RMOO 

i 

R 

XHUE-XHUO 

UE-9. 

UC*UUIC0S( IETAP >112 

9291 164 

XHC-9. 

HSL-HC 

9«7i 

FSL-PE 

928i 

TSL-TE 

9291 

XHUSl-XNUE 

6391 

USL-UC 

63H 

RMOSL-RMOC 

832i 

HD* HU 

933i 

FD-FU 

8341 

TD-TU 

6391 

DD-DU 

6391 

HD* HU 

8371 

UD-UU 

8381 

SD-SU 

839i 

AD-AU 
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S4fi 


GAAAD-CAflAU 

(41 1 


GARASL-1.4 

w: 


PSTAQ-PO 


rhostc-rhoo 

•441 


TSTG-TO 

•49i 


xrustc-xruo 

•4§i 

•4?iC 


CO TO 12! 

•48tC 
•49 iC 
SMi 


SUEPT CVLINDC9 STAC LINE PROPS 

!!• 

UT-USC0S(9ETMP1/189) 

nn 


VN-9. 

5 »• 


JMVT*LT,VC) W-S09T<UDIUO-UTtun 

•53« 


EAN»UN/AD 

••41 


TSL-TOt C 1 . ♦ . 21EHNSEAN 1 

•••« 


PSL-PDf<TSl/TD)tt3.5 

•Ml 


HSL-.24ITSL 

M7i 


RNOSL-PSL/CCRiTSL) 

•Mi 


GANASL-1.4 

M9t 


USl-UUtCOS ( IETMPI / 1 M ) 

M#> 


CPU HANSEN < XRUSL, r$L, T$L ) 

MU 


PE-PSL 

•Mi 


HE -Mil 

M3i 


TE-TSL 

M4i 


RHOE-RHOSL 

MS* 


GARAE-GARASL 

•Mi 


XRUE-XNUSL 

M7i 


UE-USL 

•Mi 


XRE-UT/AD 

569i 


PSAU-PTOTAL 



TSU-TTOTAL 

•7H 


RSU-RMOTOT 

•79i 


ANCLE -P 1/2 

•73i 


CALL DOUNID< ANCLE > 

•74i 


PSTAG-PTOTAL 

575 « 


HT0TAL»9.24tTlNf ♦UIU/59I93 . 

•79 1 


TSTC-MT0TAL/.24 

•77i 


RMOSTG-PSTAC/ ( GRtTSTG ) 

•79i 


CALL HANSEN < XRU5TG, PST AG, TSTG ) 

•79« 


PTOTAL-PSAO 

M9i 


TTOTAL-TSW 

MU 


RHOTOT-RSU 

Mti 


9CTAP-9ETAfPI/199 

m3!c 


CO TO 119 

•me 

•me 


DEFINE 90UNDARV LAVE A EDGE PROPERTIES 

••71 

129 

PE-PL 

•m 


TE-TL 

•m 


RMOC-DL 

•Mi 


ME* ML 

nn 


AC-AL 

•9f! 


GAME* CAM L 

•B3i 


UE-UL 

•941 


RE -ML 

•m 


XRE-RE 

•m 


CALL HANSEN (XRUE, PE, TE ) 

•97i 

119 

HO-MTOTAL 

•m 


PO-PTOTAL 

•m 


T0-M0/.24 

9M> 


RMOO-PO/<GRITO ) 

MU 


GARAO-1 .4 

•Mi 


CALL HANSEN < XRUO, PO, TO ) 

M3i 

M4iC 


CALL HANSEN(XRUU,PU,TU) 

•me 

•MiC 


DEFINE UAIL PROPERTIES 

M7i 

121 

HU* 9. 24 try 

•Mi 


GU-1.4 

•m 


XIIACHU-NU 

•19i 


RMOU-PE/IGRITU) 

•III 


CALL HANSEN (XRUU, PE, TU) 

•!•' 

•13iC 


ALPHA-9ETAP/PIS199 

•14IC 

•me 


ECKERT REFERENCE ENTHALPY PROPERTIES - 


'-39 




RO f >WUR-RfCCS< 1 .6#-#.##t( RNEE/RNU > ) 
ITERATE FOR RHO R AND WU R 


RNP-RORNUR/ < 1 . E- t 3IPE ) 

IF (RHP-4 • >13# # 1 3# # 131 
13# TRGESS-2.2lC4tRHPtl<-2.47) 

GO TO 132 

131 TRGCSS*2 . E71RRPtt (-7.30) 

I 132 I F ( TRC r SS . OT . 1 •### . > TRCESS* 1####. 
RM0RC*PE/ (CR1TRGCSS ) 

CALL NANSEN ( XRURC # PE # TRG£ SS ) 

RACE SS • RMORC * XRURC 
ERR*ABS< ( RRGESS-RCRRUR )/RORHUR ) 
IFCERR.lt.. # 11 GO TO 13S 
RE$ID # RORRUR-RRGESS 
TRCP2-TRC£SSI1.#S 
TRGR2*TRG£$SS# . 9S 
DCRIUR«-PC/( GRtTRCESSITRCCSS ) 

CALL HAASE N ( XRURP2 , PE , TRGP2 ) 

CALL MAhSEN ( XRURR2 , PE , TRCW2 > 
DERIUR- ( XAURP2-XHURR2 )/( . 1 tTRGESS ) 
DERIU--XRURCtDERlUR-RNORClDERIUR 
RATE-RESID/DERIY 
TRCESS-TRGESS-RATE 
I-1M 

IF ( 1 . LT . S# ) GO TO 132 
i 135 RHOR-RMORG 
I XRUR- XRURC 

i TR-TRCESS 

#7# I HR-.24STR 

#711 GAWAGM.4 

#72i 99# RETURN 

#71 
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THIS SUBROUTINE CALCULATES VISCOSITY BASED ON HANSEN 
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SUBROUTINE HATRESt HATL, ON, RHOP1, CPU, DEL, TDOT, TU, THAT, CPI1AT, NHTL ) 
DIMENSION CPHAT<10), THAT (10), XX (10) 

IF ( PIATL . GT . 0 ) GO TO 10 
T-TU-459.7 
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SUBROUTINE NOLIER<H # P#NOPT#T#Z,S#RHO# GANNA) 


11 
2!_ 

3tC NOPT-9 LOOK UP PROPS BASED ON P AND H 

4IC NOPT-1 LOOK UP PROPS BASED ON P AND S 

5«C NOPT-2 LOOK UP PROPS BASED ON P AND T 

6tC NOPT-3 LOOK UP PROPS BASED ON H AND S 

?! DINENSION FLP<33,29),HZ(33,29),TT<33,29)#ZT(33#28)#GANE<33#29)# 

•I IENTR0<33#29 >#FLP0(669 )#HZ0(669 >#TT0<669 )#ZT0(669 ).CANE0(669 ) # 

9l 2ENTR00 ( 66# ) # FLPZ ( 29 ) # HTBL < 33 > # ENTROU (29,33) # FLPU (29,33) 

lit EQUIUALENCE(FLPO,FLP>, (HZO.HZ ), <TTO # TT )# (ZTO,ZT )# (GANEO,GAN£ ), 

111 1 ( ENTROO, ENTRO ) 

121 DATA P0,CP0R,M0#S0R,G,R,ALE,RT0,CP/21 16. ,3.48158, 117.346#23.6, 

131 132 .2# S3 .35# 2 *392585 #33. 795# .23869/ 

141 CONNON/FLAG/ IDEAL# IFP 

1SI IDEAL-9 

181 Z-1.8 

17 I GANNA-1. 4 

181 IFdZ.EQ.33) GO TO 5 

19t DO 2 K- 1 ,29 

29i LL-339CK-1 ) 

21 1 DO 2 L-l#33 

22i FLPO( LL*L )-FLPZ(K ) 

23i IF(K.EQ.l) GO TO 2 

241 MZO(LL+L)-MZO(L) 

25t 2 CONTINUE 

26t IZ-33 

27t JZ-29 

28i DO 1 1-1# IZ 

29i HT1LC I )-HZ( I# 1 ) 

39t DO 1 J-1#JZ 

3H ENTROU( J,I )-ENTRO( I, JZ-J+1 ) 

321 1 FLPU( J# I )-FLPC I, JZ-JM ) 

33« 5 IF (NOPT . EO.3 ) CO TO 49 

341 PL-AL0G19(P/2116. ) ♦ 19. 

351 IF ( NOPT- 1 ) 19,28,39 

36t 19 IFCH.LT.199. ) GO TO 199 

37« CALL DINT(H,MZ,Z,ZT,PL,FLPZ,IZ,JZ,T,TT,S,ENTR0,9,GANNA,GANE> 

38t IF (Z-.2222E* 39)11, 12,12 

391 11 RHO- P/( 32. 2S53. 3SBZST ) 

49i GO TO 59 

411 199 IF (H. LT.9. ) GO TO 12 

42* T-H/CP 

431 RH0-P/(G*RIT ) 

44i 191 S- (CP0RIAL0G19(H/M0 )-AL0G19(P/P0 ) )*ALE*S0R 

45i GO TO 59 

461 12 IFP- i 

47t IDEAL-1 

481 GO TO 59 

49i 29 IF(PL.LT.FLPZ(t ) ) GO TO 22 

591 DO 21 J-2# JZ 

511 IF(PL-FLPZ(J>> 283,292,21 

52 « 21 CONTINUE 

53i 292 IF ( S-ENTRO ( 2# J ) ) 295,294,294 

541 293 IF ( S-ENTRO (2# J-l ) ) 295, 294,294 

55i 294 CALL DINT (S,ENTFO,H,HZ,PL, FLPZ, IZ# JZ#T,TT,Z#ZT, 1,GANNA#GANE ) 

56t IF ( M- . 2222E09 ) 1 1 , 22 , 22 

57* 295 M-RT0tCP0Rtl9It((S-$0R+ALESAL0G19(P/P0))/(CP0RSALE)) 

681 T-M/CP 

69l RHO-P/ (GtRtT ) 

69* GO TO 59 

611 22 IFP- 1 

92 i IDEAL-1 

631 GO TO 59 

641 39 IF(T, LT .419. ) GO TO 381 

65 1 CALL DINT(T #TT,H,HZ,PL#FLPZ, IZ, JZ, Z#ZT,S,£NTRO, 1# GANNA, CANE ) 

66i IF (H- . 2222E09 )1 1 ,31,31 

67* 391 IF(T. LT .9. ) GO TO 31 

68t M-CPtT 

69* RMO-P/(GtRtT ) 

79t GO TO 191 

7H 31 IFP-I 

72* IDEAL-1 

731 GO TO 58 

741 49 IFCH.LT.199. ) GO TO 49i 

76 1 CALL DINTKS, ENTROU, PL, FLPU,H#HT|L#JZ#IZ) 

79i IF ( PL- . 2222E+39 >41 , 42, 42 

771 41 IFCPL.GT. 12.69897) GO TO 44 

78t P- 19. BKPL-19. )S21 16. 


R-44 




791 

89t 

8H 

82 1 
83« 
841 
8Si 
881 
87t 
88i 
89i 
98i 
9lt 
92t 
93i 
94i 
9St 
96i 
971 
98t 
99t 
I8#t 
1911 
182t 
193t 
1941 
18S t 
186i 
1871 
188t 
199t 
1191 
HU 
U2t 
1131 
1141 
1161 
116i 
117i 
U8i 
119t 
129i 
1211 
1221 
123i 
124t 
125i 
126i 
1271 
128t 
129t 
138i 
131i 
132i 
133i 
1341 
1351 
136t 
137i 
138t 
139i 
1491 
141t 
1421 
143t 
1441 
145t 
146i 
1471 
1481 
1491 
1591 
1511 
1521 
153i 
1541 
1551 
1561 
1571 
158i 


GO TO 19 

491 IFlH.tT.O. ) GO TO <? 

ALP- ( SOR-S >/ALE^CPOR*ALOGI0<H/HO ) 

P-PO*10.MALP 
T-H/CP 

RHO-P/<G«RtT) 

GO TO 59 
42 IFP-1 
49 IDEAL- 1 
GO TO 59 
44 IDEAL-2 
59 RETURN 

DATA (FLPZ<J),J-l,20> /6. 39193, 6.69897. 7.0. 7.39193. 7.69897. 8. 
19. 8.39193. 8.69897. 9.9. 9.39193. 9.69897. 10.9. 19.39193. 19.698 
297. 11.9. 11.39193. 11.69897. 12.9. 12.30193.12.69897/ 

DATA (TTOC L ).L- 1.294 )/ 0.. 419.. 834. . 1616. .2339. .3932. .3478. . 

13848.. 4925.. 4175.. 4329.. 4468.. 4690.. 4880.. 5148.. 5739.. 6238.. 6565., 

26836.. 7995.. 7155.. 7279.. 7380.. 7693.. 7911.. 8199.. 8262.. 8419.. 8568., 

38723.. 8999.. 9149.. 9450., 9.. 419.. 834. . 1616. .2339. .3932. .3517. . 

43997.. 4139.. 4399.. 4469.. 4699.. 4759.. 4940.. 5260.. 5760.. 6399.. 6795., 

56984.. 7195.. 7359.. 7499.. 7596.. 7927.. 8172.. 8374.. 8559.. 8712.. 8878., 

69954.. 9250.. 9499.. 9780., 0.. 419.. 834. , 1616. ,2339. .3932. ,3546. . 

73964 . . 4224 . . 4399 . . 4563 . . 4798 . . 4865 . . 5960 . . 5346 . . 5775 . . 6359 . . 6795 . , 

87992. . 7305.. 7495.. 7648.. 7772.. 8125.. 8388.. 8595.. 8784.. 8950.. 9135., 

99324 . . 9529 . . 9759 . . 1 9930 . , 9 . . 419.. 834 . . 1 6 1 6 . . 2339 . , 3032 . , 357 1 . , 

1 4035 . . 4295 . . 4599 . . 4669 . . 4830 . . 4999 ..5180.. 5450 .. 5899 .. 6450 . . 6870 . , 

27218.. 7450.. 7640.. 7790.. 7929.. 8323.. 8604.. 8827.. 9918.. 9216.. 9495., 

39594.. 9780.. 19039.. 10359., 9., 419.. 834. , 1616 . ,2339. ,3032. , 3586. . 

44100. . 4419.. 4639.. 4815.. 4985.. 5150.. 5340.. 5594.. 5999.. 6498.. 6995., 

57353 . . 7620 . . 7849 . . 8915 . . 8 1 72 . . 8597 . . 8901 . . 9150 . . 9378 . . 9576 . . 9765 . , 

69963. . 19180.. 19429.. 10739.. 0., 419., 834. , 1616. ,2339. ,3932. , 3596. , 

74159 .. 4485 . . 4735 . . 4930 ..5100.. 5270 .. 5480 .. 5724 . . 6950 . . 6589 . . 7059 . , 

87470.. 7769.. 8999.. 8190.. 8334.. 8811.. 9144.. 9410.. 9648.. 9846., 

919960.. 10289.. 10509.. 10740.. 1 1040.. 9.. 419.. 834.. 1616.. 2339.. 3932,/ 
DATA ( TTO ( L ) , L - 205 , 499 ) /36 1 8 . . 4298 • , 4555 . , 4829 . , 5040 . , 5249 . , 5 4 1 9 . . 

15619.. 5872.. 6195.. 6650.. 7165.. 7596.. 7908.. 8150.. 8355.. 8532.. 9036., 

29414.. 9657.. 9927.. 10160.. 10380.. 10590.. 19826.. 11 100.. 11499., 0., 

3419., 834 . , 16 1 6 . , 2339 . , 3032 . , 3639 . , 4259 . , 4655 . , 4969 . , 5292 . , 5498 . . 

45649 . . 5849 . . 6966 .. 6350 .. 6790 . . 7250 . . 7722 . . 8995 . . 8360 . . 8600 . . 8725 . , 

59378. . 9756 . . 10960. . 10359. . 19580. . 10899. . 11049.. 1 1309. . 11580. , 

611889., 9., 419., 834., 1616., 2339., 3032., 3640., 4300., 4759., 5969., 

75323 . . 5559 .. 5758 .. 5990 . . 6228 . . 6505 . . 6990 ..7319.. 7839 ..8215.. 8529 . , 

88760. . 9999. .9639. . 10949. . 10399. . 10660. . 19920. . 11180. . 11420. , 

91 1680.. 11979.. 12299.. 0., 419., 834 .. 1616. ,2339. , 3032 ., 3654 ., 4314 . , 

1 4809 .. 5 1 50 .. 5445 .. 5683 .. 59 10 .. 6159 .. 6394 .. 6659 .. 7945 .. 7485 .. 7956 . , 

28350 . . 8680 . . 8955 . . 9198 . . 9864 . . 10350 . . 10699 . . 1 1909 . . 1 1299 . . 1 1 559 . , 

31 1810.. 12199.. 12389.. 12799.. 0.. 419., 834. , 1616. , 2339. , 3032 ., 3654 . , 

44328.. 4889.. 5295.. 5598.. 5880.. 6119.. 6355.. 6628.. 6909.. 7259.. 7669., 

581 18.. 8559.. 8910.. 9219.. 9468.. 19239.. 19740.. 11 150.. 11 480.. 11790., 

612080. . 12370. . 12659. . 12989. . 13399. .0. .419. .834. .1616. .2339. .3932. , 

73654 .. 4416.. 4945 . . 5395 . . 5724 ..6915.. 6285 . . 6540 . . 6808 ..7190.. 7450 . , 

87819.. 8279.. 8695.. 9085.. 9415.. 9792.. 10519.. 11070. .11500.. 11879., 

912200. . 12519. .12829. .13119. .13459. .13800., 0., 419., 834. ,1616./ 
DATA ( TTO ( L ) , L • 40 1 , 590 ) /2339 . , 3932 . , 3654 . , 4455 . . 4989 . , 5475 . , 5839 . , 

16150.. 6459.. 6735.. 7992.. 7390.. 7640.. 8000.. 8433.. 8848.. 9245.. 9600., 

29900.. 19890.. 11 400.. 11889.. 12289.. 12639.. 12970.. 13289.. 13619., 

313969.. 14310., 0., 419., 834. , 1616. ,2339. , 3932. , 3654. ,4469. ,5910. , 

45548 . . 5976 . . 6348 . . 6655 . . 6950 ..7281.. 7575 ..7910.. 8275 . . 8658 ..9100., 

59496.. 9865.. 10188.. 11209.. 11880.. 12429.. 12820.. 13279.. 13609., 

613969. . 14319. .14699.. 15089. .0. .419., 834. , 1616. ,2339. , 3032 ,3654., 

74483 . . 5925 . . 5589 . . 6080 . . 6474 . . 6839 . . 7150 . . 7488 . . 7800 ..8140.. 8500 . , 

88892 .. 9300 . . 9799 . . 1 9070 . . 1 0404 . . 1 1 509 . . 1 2269 . . 1,» 49 . , 1 3300 . . 

913770. . 14149. .14500.. 14889. .15270. .15669., 0., 419., 834. ,1616.. 

12339.. 3932 . . 3654 . . 4497 . . 5038 . . 5675 ..6169.. 6589 . . 6975 . . 7340 . . 7686 . , 

28098 . . 8359 . . 8740 ..9198.. 9525 . . 9939 . . 10310 . . 10683 ..11840.. 12650 . , 

313280.. 13810.. 14270.. 14709.. 15100.. 15490.. 15890.. 16290., 0., 419., 
4 834., 1616.. 2339., 3032., 3654., 4505., 5959., 5725., 6266., 6750., 7188., 

57566 . . 7956 . . 8390 . . 8670 . . 9050 . . 9450 . . 9859 . . 1 9269 . . 1 9660 ..11952., 

612299.. 13199.. 13919.. 14590.. 15910.. 15489.. 15930.. 16350.. 16800., 

717220., 0., 419., 834., 1616.. 2339., 3032., 3654., 4515, 597S., 5759., 

86336 . . 6848 . . 7229 . . 7755 ..8154.. 8549 . . 8945 . . 9310 . . 9720 ..10119., 

919549.. 19940.. 11340.. 12639.. 13630.. 14400.. 15959.. 15619.. 16110,/ 
DATA (TTOC L),L-591,660)/16699., 17970., 17520., 17980., 0., 419., 834., 

11616.. 2339.. 3032.. 3654.. 4425.. 5110.. 5774.. 6386.. 6935.. 7440.. 7905., 

28352.. 8755.. 9155.. 9579.. 9999.. 19400.. 10846.. 11239.. 11639.. 13910., 

314090. . 14940 . . 15640. . 16249. . 16830. . 17320. . 17800. . 18300 . . 18820. .0. , 
4 419., 834., 1616., 2339., 3932., 3654., 4425., 51 10., 5790., 6440., 7030., 

57589.. 8999. .8599. .9049. .9499. .9905. . 10370. . 19809. . 11210. . 1 1640. , 

612069. . 13540. . 14719. . 15689. . 16489. . 17140. . 17770. . 18369. . 18910. , 

719459.. 19980,/ 


^ 4 * 


>M5 


1591 DATA <ZTO<L>, 1*1, 2491/6(1. 0,1. 007, 1.024, 1.049, 1.076,1 

1691 11. 156. 1. 189, 1. 291, 1.210, 1.216. 1.226.1. 238, 1.2S1, 1.265, 

1611 21. 356. 1.422, 1.486. 1.SS1. 1.617, 1.681, 1.746. 1.819, 1.974, 

1621 31.904,1 .929, 1.044, 1.971, 1.098. 1.124, 1.1S2, 1.176, 1.197, 

163i 41 .224, 1 .235, 1 .248, 1 .269, 1 . 275 , 1 .299, 1 . 352, 1.416,1. 489, 

164i M .674,1 .718, 1.892, 1.86S.1.92S, 6(1. 9. 1.993, 1.918, 1.949, 

1661 <.l . 119.1 .147,1.171.1.103.1.296.1.213.1.221,1.232,1.244, 

IM> 71 .,'f«7.1 .340. 1.41 1.1. 475.1.539. 1.694. 1.667. 1.732, 1.795, 

167i Rf.fl .9.1 .991. 1. 918, 1.037. 1.960. 1.088, 1.1 IS. 1.1 49,1. 16S, 

168i 91. 212, 1. 219, 1.229.1. 241. 1.2SS. 1.268, 1.283. 1.343, 1.406, 

169i 1 1. S97, 1. 661, 1. 725, 1. 788, 1.8S1.1. 999, 7«1. 0.1. 913, 1.932, 

1791 21.198,1.132,1.156,1.181.1.197,1.298.1.216.1.226,1.237, 

17H 31.277,1.337,1.398,1.461,1.526,1.589.1.652,1.715.1.778, 

172* 4781.0, 1 .012, 1.928, 1.952, 1.076, 1 . 190,1.128, 1.152,1. 175, 

173) 51. 214, 1.223, 1.234, 1.246, 1.259, 1.273, 1.332.1. 393, 1.45S, 

1741 61. .,44,1. 707. 1.770, 1.831. 1.891, 711.0,1.911, 1.924, 1.047 

1751 71.120,1.144,1.169,1.188,1.201.1.211,1.229,1.231,1.242 

176) 81.326,1.387,1.449.1.511,1.574,1.636,1.699,1.761.1.822 

177l 91.010.1.022,1.942,1.063,1.088,1.112,1.136,1.169,1.180 

178< DATA <ZT0<1 1,1*250, 5021/1.215, 1.226, 1.237, 1.249.1.263, 

179* 1 1.439, 1 .500, 1 .563,1.625, 1.687, 1.748, 1.898, 1.867,7(1.9, 

190) 21.037,1.958,1.980,1.104.1.128,1.152,1.172.1.190,1.291, 

181 i 31 .232, 1 .244, 1.258,1.312,1. 372, 1.431,1. 492. 1 . 554, 1.616, 

188) 41. 798. 1.8SS. 7(1. 0,1. 998. 1.015. 1.932, 1.052. 1.976, 1.098, 

1831 51.164,1.184,1.197,1.207,1.216,1.228,1.240,1.252.1.396, 

1841 61.484.1.545.1.606,1.667,1.727,1.787,1.844,7(1.0,1.007, 

185) 71.046,1.066,1.088.1.110,1.133.1.154,1.173.1.188,1.200, 

186) 81.232,1.245,1.297,1.353,1.412,1.472,1.531.1.592.1.653, 

187) 91 .826,7(1 .9. 1 .005, 1 .009, 1.022, 1 .940, 1.058. 1.080,1 . 101, 

188) 11.164,1.180.1.193,1.294.1.215.1.227.1.238,1.289,1.345, 

189< 21.522.1.581.1.649,1.699,1.757,1.812,7(1.9.1.094,1.907, 

1901 31.052.1.972,1.992.1.115.1.135.1.154,1.172,1.186,1.197, 

191> 41. 231. 1.281, 1.337, 1.394, 1.452, 1.519, 1.569,1. 627, 1.68S, 

192< 57(1.0,1.003,1.005,1.913.1.027,1.044,1.963.1.981,1.102, 

193* 61 . 159, 1 . 174, 1.188,1.199,1.210.1.221,1. 279, 1 . 324, 1 . 379, 

194) 71 .552, 1.609, 1.666, 1.721, 1.777,7(1.0. 1.002, 1.004, 1 .012, 

195) 81.055,1.072,1.092,1.112,1.129,1.148,1.164,1.177,1.189 

196> 91.261,1.314,1.368,1.424,1.480.1.537,1.594,1.650,1.705. 

197) DATA (ZT0(l>, 1*503, 6691/1. 001, 1.003. 1.008, 1.019, 1.030, 

198> 11.082,1.098.1.117,1.136.1.152,1.165,1.179,1.191,1.203, 

199) 21.355,1.410, 1 . 466, 1 .522, 1 .577, 1 .632. 1 .687, 1.739.8(1 .0. 

200) 31.014,1.024,1.037,1.052.1.069,1.085.1.103.1.119,1.136. 

201 ! 41.176,1.189.1 .235, 1 .285, 1.338, 1.390, 1 .445, 1 .500, 1.553, 

202) 51.712.8(1.0,1.001,1.003,1.010.1.019,1.030.1.044,1.059, 

203> 61.107.1.122,1.136.1. 150, 1 . 163, 1.175,1 .223, 1.272, 1 . 322, 

204) 71.480.1.534.1 .586. 1.639, 1.692. 9(1. 0.1. 001, 1.007, 1.015. 

205) 81 .050,1.064.1.079,1.093,1.108,1.122,1.136,1.148.1.161, 

206) 91.306.1.357.1.410.1.460,1.512.1.564.1.615.1.665,10(1.. 

207) 11.019.1.028,1.039,1.050,1.063,1.076,1.090,1.104,1.116, 

208) 21.187,1.234,1.282.1.332.1.382,1.431,1.480,1.531,1.581. 

209) DATA (GAMEOUl. 1*1, 2121/2(1. 4, 1.390, 1.344, 1.317,1. 295, 

210) 11.114,1.108,1.108,1.114,1.131,1.169.1.226,1.264,1.238, 

211) 21.121,1.110,1.103,1.099,1.094,1.092,1.093,1.094,1.996. 

212) 31.111,1.122,1.138.2(1.4,1.390,1.344.1.317.1.295,1.214. 

213) 41.113,1.113,1.118.1.131,1.162.1.212.1.262,1.248,1.190. 

214) 51.116.1. 109, 1 . 10S, 1.098, 1.096, 1.097, 1.098.1. 100. 1 . 192, 

215) 61.124,1.138.2(1.4.1.390.1.344.1.317,1.295,1.220.1.154, 

216) 71.116,1.121,1.132,1.159,1.203,1.258,1.255,1.198,1.158, 

217) 81.114,1.109.1.101.1.099,1.100.1.101.1.103,1.105.1.110. 

218) 91.139,2(1.4,1. 390, 1.344,1.317,1 .295, 1 .227, 1.160,1.139, 

219) 11.124,1.134,1.156.1.195,1.250.1.261.1.206.1.165.1.140, 

220) 21.114.1.105.1.102,1.103.1.104,1.106.1.108.1.112,1.118. 

221) 32(1.4,1.390.1.344,1.317.1.295.1.236.1.169.1.138,1.128, 

222) 41.137,1.154,1.186.1.237,1.261.1.218.1.174,1.148.1.134. 

223) 51.110.1.107.1.107.1.108,1.110,1.112,1.116,1.122,1.130, 

224) 61.390.1.344,1.317,1.295,1.244.1.176,1.144,1.133.1.130, 

225) 71.154.1.181,1.226.1.258,1.227.1.182.1.156,1.140,1.131. 

226) 81.111,1.111,1.112,1.113,1.116,1.120,1.125,1.133,1.145, 

227) 91. 344. 1. 317, 1.29S. 1.250, 1.183, 1.151, 1.138, 1.134, 1.1 37, 

228) DATA ( GAHEOl L >, L*213,276 >/ 1.179,1.218.1.252.1.236. 

229) 11.147.1.138,1.131,1.118.1.114.1.114,1.115,1.117.1.120, 

230) 21. 136, 1. 147. 2(1. 4. 1. 390, 1.344, 1.317.1. 295. 1.2S8, 1.193, 

231 1 31.140,1.142.1.148.1. 157, 1.177,1.209.1. 244, 1.244. 1 .207, 


103,1.130. 

1.289,1.296, 

1.934,6(1.0, 

1.208,1.214, 

1.545,1.609, 

1.067,1.994, 

1.259,1.272, 

1.858,1.918, 

1.189.1.204, 
1.470,1.534. 
1.0S6. 1.081, 
1.249,1.264, 
1.840,1.899, 

1.193.1.205, 
1.518,1.582. 
1.979,1.996, 
1.256,1.269. 
1.880,7(1.0, 
1.196,1.206/ 
1.318,1.378, 
1.009,1.017, 
1.211,1.229, 
1.677.1.738, 

1.121.1.144, 
1.364.1.424, 
1.012,1.926, 
1.210.1.229, 
1.711,1.770, 

1.124.1.144, 
1.493,1.462, 
1.918,1.034, 
1.298.1.219. 
1.742,1.798, 

1.122.1.140, 
1.437,1.494, 

1.923.1.037, 
1.209,1.213, 
1.759,7(1.0/ 
1.046,1.964, 
1.251.1.392, 
1.992.1.005, 
1.150,1.164, 
1.607,1.660, 
1.075,1.091, 
1.375,1.428, 

1.024.1.037. 
1.208,1.256, 
1.903,1.010, 

1.128.1.141, 
1.630/ 

1.206.1.143, 

1.180.1.143, 
1.099,1.104, 

1.149.1.119, 

1.152.1.128, 

1.107.1.114, 

1.124.1.117, 
1.133.1.121. 
1.116,1.126, 

1 . 121 . 1 . 120 . 

1 . 126 . 1 . 118 , 
1.128.1.149, 

1.125.1.129, 
1.125,1.120, 
1.142,2(1.4, 
1.133,1.140. 

1.126.1.114, 
2(1.4,1.390, 
1.143,1.155/ 
1.191,1.163, 
1.123, 1.128, 
1.161,1.146, 
1.177,1.159. 
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2381 31.148, 

2331 41. 138. 1.124, 1.128, 1.1 10, 1.129, 1.122, 1.124, 1.128, 1.133, 1.141, 1.158, 
8341 52*1. 4, 1.399, 1.344. 1.317. 1.295, 1.263. 1.291, 1.169. 1.152, 1.148, 1.147/ 
236t DATA (GAREOd J,L*27?,424 )/ 1.152, 

2381 61. 169, 1.177, 1.294, 1.237, 1.245, 1.219, 1.189, 1.169.1.156.1.146.1.128, 
2371 71 . 124, 1 . 123, 1 . 124, 1 . 126, 1 . 128, 1.132,1.137.1.143,1. 154.2*1.4, 1 . 399. 
238* 81.344,1.317,1.295,1.267.1.219,1.177,1.159,1.152.1.152.1.156,1.163, 
239< 91 . 178. 1.299, 1.239, 1.242, 1.226,1.291,1.179.1. 164. 1 .152. 1 . 134, 1 . 128, 
2491 11.127.1.128,1.139,1.132.1.136.1.149,1.147,1.156.2*1.4,1.399,1.344, 
2411 21. 317, 1.29S, 1.271, 1.223, 1.188. 1.168. 1.169, 1.159. 1.161, 1.169, 1.189. 
2421 31.197,1.221,1.236.1.231,1.214,1.192,1.174,1.162,1.141,1.135.1.133, 
2431 41.133,1.136.1.138,1.142,1.146,1.152,1.161,2*1.4,1.399.1.344,1.317, 
2441 51.295, 1.273, 1.232, 1.196,1.175.1. 167, 1 . 164, 1.166,1. 173, 1 . 182. 1.196, 
245* 61.217,1.231.1.232.1.222,1.292.1.183,1.179.1.147,1.149.1.138.1.138, 
246* 71.149,1.142,1.146,1.159,1.156.1.164,2*1 .4,1.399. 1.344.1.317.1.29S, 
247* 81.274,1.238,1.296,1.183,1.172,1.179,1.171,1.178.1.186,1.197,1.214, 
248* 91 .228, 1 .232, 1 .226, 1.219. 1.1 92, 1.179, 1.1 S3, 1.1 45, 1.143, 1.143,1. 145/ 
249* DATA<GAIC0<IM.*42S. 636 1/1.147, 1.15, 1.154, 1.169,1.168,2*1.4,1.399, 
258* 11.344,1.317,1.295,1.275,1.244,1.218,1.195.1.182,1.178,1.178,1.183, 
251* 21. 198, 1.299, 1.213.1.22S, 1.231. 1.228. 1.218, 1.294, 1.191, 1.162, 1.153. 
252* 31. 158,1. 158, 1.1S9, 1.153, 1.156, 1.169, 1.168, 1.173, 211. 4. 1.399 1.344, 
2S3* 41 .317, 1 .295, 1 .275. 1.247, 1.226,1.294,1.189,1. 184, 1.184,1. 188, 1 . 194, 
2S4* 51.292. 1.213.1.224,1 .239, 1.238, 1 .224. 1 .212, 1.299, 1 . 169, 1 . 159. 1 . 155, 
255* 61.155,1.156,1.158,1.161,1.165,1.171.1.177,2*1 .4.1.399,1.344,1.317, 
2S6* 71. 295, 1.275, 1.259, 1.232, 1.21 1,1. 196, 1.199, 1.198. 1.192, 1.198. 1.295, 
257* 81.214,1.224,1.239.1.231,1.227.1.219,1.298,1.176,1.165,1.161.1.169, 
258* 91.161.1.163.1.166.1.169,1.175.1.182,2*1.4,1.399.1.344,1.317,1.295. 
2S9* 11.276,1.253.1.237,1 219, • . 295. 1 . 199, 1 . 198, 1 . 199, 1 .294, 1 .219, 1 .216, 
269* 21. 224. 1.239, 1.232, 1.231, 1.226. 1.216. 1.166,1. 174, 1.169, 1.167, 1.168, 
261* 31. 1 79. 1. 173, 1. 1 77, 1.181, 1.1 87, 2t 1.4,1. 390. 1 . 344, 1.317, 1.295,1.276, 
262* 41. 254. 1.240, 1.223. 1.212, 1.205, 1.203, 1.204, 1.208,1 214, 1.229. 1.22S. 
263* 51.230, 1.233.1.232,1 .228,1.221.1.194,1.181,1.176,1.173.1.174.1.176. 
264* 61.179, J .1P2. 1.186. 1.192, 2*1. 4. 1.390, 1.344. 1.317, 1.295, 1.276, 1.255, 
265* 71.241,1.226,1.217,1.211,1.208.1.209.1.212.1.218,1.223,1.226,1.229, 
266* 61.233,1.232,1.229,1.224,1.203.1.188,1.182,1.179,1.180,1.182.1.184, 
267* 91. 187. 1.191, 1.196,2*1. 4. 1.398, 1.34-*, 1.317, 1.295. 1.276, 1.255, 1.242/ 
268* DATAICAAEOU ), L-637,668 >/l . 228, 1.223. 1.22, 1.214,1. 218.1. 219,1 .223. 
269* 11.228,1.228.1.238,1.233.1.233,1.229,1.227,1.214,1.198,1.191.1.187. 
278* 21.188,1 .189, 1.192. 1.19S, 1.198, 1.282/ 

271* DATA (ENTROOd >,1-1,244)/ 23.68.31 .56,33.97,36.43,37.92,39.81, 

272* 139.88,48.87,41 .78.42.65,43.58.44.33,45.16.45.94,46.67.47.33,47.94, 
273* 248.58,49.94,49.58,59.10.58.59.51.98,53.85,54.88,56.70.58.46.69.16. 
274* 361 .93,63.61,65.22,66.85,68.48,22.76,30.64,33.85,35.51,37.99,38.08, 
275* 438. 96, 39. 93. 48. 84, 41. 78, 42. 54, 43. 32, 44. 88, 44. 84, 45. S6. 46. 22. 46. 81, 
276* 547.39.48.01.48.41,48.95.49.43.49.93.51.84,53.58,55.32,57.08.58.72, 
277* 668. 3S.62. 00. 63. 58, 65. 10. 66. 66. 22. 87. 29. 95. 32. 36. 34. 82. 36. 31, 37. 40, 
278* 738. 25, 39. 21. 49. 11, 49. 99. 41. 79. 42. 57. 43. 39, 44. 01. 44. 75. 45. 40, 45. 99. 
279* 846.53.47.16,47.56.48.02.48.51,49.02,50.85.52.61,54.32.55.99.57.64, 
288* 9S9. 24. 60. 80. 62. 33, 63. 84, 65. 34, 21. 37, 29. 26, 31. 66. 34. 12, 35. 61, 36. 70. 
281* 137 .56, 38 .59, 39 .40, 40 .24. 41 01.41.79,42.52.43.22,43.92.44.56.45.15, 
282* 245.70,46.28.46.70,47.17,47.65.48. 12,49. 94, 51. 65, 53. 31, 64. 96. 56. 52. 
283* 358.06.59.61,61.10.62.57.63.99.20.46,28.34,30.75.33.21,34.70,35.79. 
284* 436.63,37.57.38.43.39.23.40.01.40.77.41.45,42.14.42.82,43.44,44.02. 
285* 544. 57. 45. 09. 45. 55. 46. 03, 46. 49, 46. 93. 48. 73,50. 3S, 51. 98. 53. 54, 55. 06. 
286* 656.56.58.04.59. J6. 60. 86, 62. 28. 19.77,27.65.30.06.32.52,34.01,35.10. 
287* 735.96,36.89,37.72,38.53.39.27.40.00,40.72,41.36.42.90.42.63.43.20, 
288* 843.74,44.23,44.79.45. 18.45.63,46.04.47.77.49.38,50.94.52.49,53.97. 
289: 955.42,56.87.58.23,59.60.60.98,19.07.26.96.29.36,31.82,33.31,34.40/ 
890* DATA (ENTROOd >.1-205, 408)/ 35. 27. 36. 17, 37 . 01 , 37. 80. 38. 51 , 39. 22. 
291* 139.90.40.57,41.17,41.80,42.37.42.90,43.39,43.86.44.32,44.76,45.17, 
292* 246 .80, 48. 42, 49. 93, 51. 40, 52. 87, 54. 31, 55. 66 , 5/. 00, 58. 36, 59. 68 , 18 . 16. 
293* 326.04,28.45.29.91 , 32. 40, 33 . 49, 34 . 35, 35 . 25, 36 .08, 36. 65, 37. 56, 38 . 25. 
294* 438.88,39.53,40. 12,40.64,41 .25.41 .89,42.30,42. 74,43.16. 43. 61, 44. 80, 
296* 545. 61, 47. 11, 48. 62. 50. 00, 51. 40. 52.80, 54.) 3, 55. -<0.56. 68, 57. 90, 17.46, 
296* 62S. 35. 27. 75. 30. 21. 31. 70. 32. 7g, 33 . G4, 34. SS. 35. 38. 36. 12. 36. 81. 37. 50. 
297* 738.13,38.75,39.33.39.90,40.43,40.95.41.42,41.89,42.32,42.72,43.13, 
298* 844.70,46.19,47.62.49.00,50.33,51 . 64. 52 . 97. 54 .24 , 55 . 42,56 . 60. 16.77. 
299* 924.65,27.06,29.52,31 .01,32.10,32.96.33.88,34.66,35.40.36.08.36.72. 
30# * 137.37,37.97.38.54,39. 11.39.63,40. 12,40.60.41 .02,41 .48.41 .88,42.27, 
341 * 243.79.45.25.46.64,47.96.49.28.50.52,51 80.53.00.54.20,55.32, 15.85, 
342* 323.74,26.14.28.60,30.09,31.18,32.03.32.97,33.75,34.45.35.11.35.76. 
343* 436. 35, 36. 93,37. 48, 38. 04, 38. 54. 39. 02. 39. 48. 39. 94, 40. 35. 40. 75. 41. 13, 
344* 542. 60. 43. 99. 45. 33. 46. 65. 47. 87, 49.06, 50. 2b. 51 .45.52.58.53.67,15. 16. 
345* 623.04,25.45,27.91,29.40,30.49,31.34,32.26,33.04.33.74.34.41,35.00, 
346* 735.61.36. 16, 36.69, 37. 23. 37. 73. 38 . 1 9, 38 .66, 39 . 10, 39 . 50, 39 . 8S. 40 . 28 . 
347' 841.70,43.04,44.35.45.60,46.80,47.98,49. 10,50.27,51.35,52.40,14.47, 
348* 922.35,24.76.27.22,28.71,29.80,30.67.31 .57,32.35.33.02,33.68,34.28/ 
349* DATA (ENTROOd ) ,1-409,612 >/ 34 . 8S, 35 . 40, 35. 92. 36 . 45, 36 .92. 37 . 40, 
314* 137.84,38.27,38.68.39.03,39.40,40.81.42.13,43.38,44.56,45.73,46.88, 


( • f 1 
«.*-*» / 


311 1 247.98.49.07,59. 15,51.15,13.55,21 .43,23.84,26.39,27.79,28.88,29.7s. 

3121 339.64,31.42,32.10.32.72,33.30,33.87,34.39.34.90,35.40,35.86.36.32. 

313i 436.75,37.17,37.57.37.95,38.30,39.65.40.93,42.10.43.27,44.34,45.45. 

3141 546. 51. 47. 53. 48. 53, 49. S0, 12. 86, 20. 74, 23. 15, 25. 61. 27. 10, 28. 19,29.01. 

315* 629.97,30.71,31.39,32.00,32.60.33.15,33.68,34.14,34.62.35.10.35.54, 

3161 735.95.36.35.36.77,37.11,37.44.38.77,40.02,41.16.42.27.43.32.44.37, 

317< 845.42,46.38.47.36.48.28.12.16.20.05.22.45.24.91,26.40.27.49.28.36, 

318> 929.25.30.93.39.68,31.30,31.86,32.41.32.93,33.40.33.88,34.31,34.73, 

3191 135.14,35.53,35.91,36.27,36.62.37.90.39.10.40.22,41 .31,42.34,43.31. 

3201 244.28,45.27,46.18,47.09,11.25.19.13,21 .54 .24 .00,25 . 49. 26 .58. 27 . 44, 

32H 328.34,29.11,29.77,30.37,30.95,31 .46.31 .96,32.41.32 88.33.30.33.72. 

3221 434.11.34.50,34.84,35.19.35.53.36.77,37.90,38.99,40.00,41.00,41.95, 

3231 542 . 87, 43 . 79.44.63, 45 . 49, 10.55, 18. 44, 20 . 8* ,23 . 30,24 . 79, 25 . 88. 26 . 75, 

3241 627.64,28.42,29.06,29.69,30.22,30.76.31 .22,31 .68,32.12,32.54,32.95, 

325i 733. 31, 33. 68. 34. 05. 34. 38. 34. 71. 35. 93. 37. 06. 38. 05. 39. 08, 40. 00, 40. 92, 

3261 841.80,42.64.43.51.44.30. 9.86.17.76,20.15,22.61,24.10,25.19.26.08, 

3271 926.95,27.72,28.37,28.98,29 . 53. 30 . 04, 30.50, 30 .95, 31 . 40, 31 . 80, 32 . 20/ 

328* DATA <ENTROO< L ) , L *61 3, 660 )/ 32 .SS. 32 . 91 . 33 .24. 33.59, 33 .90, 35.09, 

329* 136. 16, 37. 18, 38. 12, 39.0', 39 . 90, 40 .76, 41 .57, 42 . 37, 43 . 18, 8.95, 16.83, 

3301 219.24,21 . 70,23. 19,24 ,28. 25 . 15,26 .03,26 . 79, 27 . 44 , 28 . 03, 28 . 59, 29. 10, 

331* 329.56,30.00,30.41,30.80,31.18.31 .54.31 .89.32.23,32.54,32.84,33.94, 

332* 435. 01, 35. 96, 36. 90. 37. 77. 38. 57, 39. 33. 40. 12, 40. 92, 41. 68' 

333* DATA (HZO(L ),L«1, 33 )/0., 100., 200. .400. ,600., 800., 1000., 1250. , 

334 * 1 1 500 . , 1 750 . . 2000 . , 2250 . , 25 A0 . . 2750 . . 3000 . . 3250 . , 3500 . , 3750 . , 4000 . , 

33S * 24250 . , 4500 . , 4750 . . 5000 . , 6000 . . 7000 . . 8000 . , 9000 . , 1 0000 ..11000., 

336* 312000. ,13000. .14000. .16000./ 

337* END 




a- 48 


CL 

►4 % 


C 

m 


LU 

cu. 


K 

300 


»► 

o x 5 

00 


• UUI 

*-« 


n i-a 

« IT 

in 

in ^-4 

»- CO 

-j 

*oo in 

I 

K 

* • 

r 

a 

3 VOX 

•3 O) 

t 

«L E CJ 

W CO 


£ ■ K 

z in 

r 

a ■ ho 

M \ * 

M 

i- r -4 r 

*00 • 

1 

H 

~ '00 

4 - 

1 - * * • 

•> «ru 

-J 

inr 4 

<r %in 

M 

3 ® 3 — 

r ^ • 

■3 

▼ co •— 

n *n 



® 

® 

* 


® ® 
* * 


CD 9 
<9 <$• 

8? 


® 4 > 
** «* 

% * ru ru 


ru 

. . 


® *®® 

® ® 

* . 



. . 

OM 1 

® ® 


® ®rxs> 

® ® 





z 'in 

W M 


m •• m ® 

00 CO 

ru 

. 

. 


ttco • 

ruru 


*fur\jn 

ru 

in® 

® 

® ® ® 

run 

• ' . 

« « 

« 

• ' • 


com 

*4 

**oo-*a»oo 

® r 

~ oo ru 

. • 



• • 

(U *4 

m 

m 'c ru oo oj 

mm 

M **4 

« in in 

tt 

® 

• ® ® ® ® 
® ® ® m® 

® ® 
® ® 


% 

.... 

• • 

H ' » 

r* m 

M 

m in oo co in 

m ▼ 




% 


-> ru • 
w *irs 

ruruoo • 
u *fur- 

h * IA *•* 

2 no m m T 

* 

* V 


m ru — * ~ nj ru 


r- oo 
r- co 


» ! 


CO 

CS * 


• W CO - <T «(DMS 

• n \ k 
«cl^ «® -> ru 

hOH C • so * » % * « 

n »«Sia« .... % . 

'~ C 


SSSIS.SS 

ru®r-ru®**r-ri 
^ ru — ' ru ▼ 
*••••%•% ru ■* • -« 

«(£«(£ 

99999999 ^ (vnwn 
dd®®®dind~rU-4 l/w« 
* ini/>r^minruruou • • 

•-» W* r. <!»»*• 

CO 

“> cor- cn — 


* * • » s ® ® 

«x> run m m 

ru • • ® r- ru ru 

in go ci) run 


cn * oo 

r- go ru/' • co • 

(7> *4 CO C7> ®n- * • ® 

* ru oo ru run ® n #4 

* r> 


~«® m 
* 00 CO 
▼ 00 

® • r- 
*^in •< 

in 

r- r ^ cn® * » * * 
oo in oo -^r *• ® in 

rvi * -« r* ru ▼ n ▼ cu«~» 


— • vn r- co ▼ 

*r ru in — cr oo o> ▼ 

— corururumv 


r- 

00 •Ok 

• in ® • 
n • «i ^ 

cn^in -4 


in 

ru 


* % 
® ® 


• S 

o»ru 

Zt 


in • ® ® 
cu® *in 
. tt -*ru 
co ru ▼ . 


® ru® ® 

• ® m ® m 
> in • • o 

• • cn o -- 
"» * — **• • 


®® UO® 00 

r •u.ini/i 
• in • •* 
▼ • in co — 1 
rur in • 


28 


% % 

• I** (U 

• id 

JS3 


® ®ry ® 
ooinl® 

• in • *® 
cr> . ▼ ^ — 
rum *-«r* • 


ru< 
in c* - 


^ ® ® ® ® * *®®®®®®®® • ® m ® m . • ® in cti win oo *n 

mo>«*®«£ocru- 4 inr-i^m * 
— <1 * • in «_••«*’ ® m (V) C ® r\j 4 ) <*> f- — . nj ® *4 *» ® n in ru ru ruru — 

■ 

h “1 

*6 


~ r r-mco ® $ $ lj 

• VODOU) • • 00 4> ® ® oo 
a • r- «m® ® • m • • m 
h r* . r- • m m m • ~ co ® 
win oo *ii/» rn — co co co ru oo •-« w • 


« 

c 


N 

a 

x 

s 

o 


ISI 

a 

o 

% 

a 

o 


■ h * 4-7 . -4 i i ru 

~ ^ * 4*1 « « « 


» CO — ^ — -4 


m ru 
. — ▼ mmru 


coccwcti — ®oo®in ® 
- — r-i . . • m r- ru ▼ 


•>. — ® ® ® ® •®G>®®®®®®®^)Ooc 7 ico®ao« > r-r>u r-®®^® • • • oo — rn . •* > ® 

“-> cm *~>in®®®«®®®® -»r>Ln- 4 &cDii>r\ir .~*iD*in«~c 7 *~»*<ucMn— r\jruc*r^~* • 

two«wr\jn»(\j®Jiiflru/» m ru ru ®» 4 ^®ru — n'ururu^* - 

i l^ ru a. rg***%44^**** 

“KJSSSSSSS'iu, 



g *> in ' » j) ® ® ® oc • 

ago msBum oo n noaia • • n — -*in -tcoja) 

f* ••• _• • ot ® 

^J?U« -4*®*® ® • — *4^ ru-4®®<l-^ 

.OOQmO II III III l I I • o 


-lyn^muirawnHrunt -*rum mior oocti •• —rur 


vO ▼ n* ruru® ▼ nj 

— rum^in®r ^®®##-4 


• • . V 

N 

• * * h . 

£ ®S*J l- X 

• • • J w • X 

»-NOJl 

• * •«om®o 

8 -5 i • 

♦ • hJO- 


uuououuuu 


— rurn^ mcor^w®®— run^incor- 
^ 1 1 * * 'rtfii/uAinininmin 


38SS3SSS«» 


8-4»a 


u*b 

u. *ujh 

56 

'flLWXU 
ISIO I %0 
* • «JO 

a« Sujh 
K fttt -o 

CL 3 Q.|fct Ld 

ujuj *h33 

s *« « •< 

M- -LA&1- 

w28»5« 

P&gaoN 
% *5 wan 
CUN «&U.K 
OJN -X W 

;&!g 

HOJuiO 
^ZO « « « 

o««o 



cu 

N 

9 * 

CU 
- H 
O X 

«4 

-LA 

a 

n 


-IX 

. J -Ui 

— 9 0 

w 

X 

u. o 

u •"* - 

a. — cl 

«>4 

«* 

•» X 

-r 90 

•3 w O 

•SI 

1 

X Or Or O — -J 

• *n iah - 


a 

x 

(L 
• -» 
~ <T 
CU • 

wft) 
•“*“7 
O — 
O-I 

*a 
—» h 

ruul 

->«i 


& 

i 

a 


i 

* 


wutoi 
^acAi^ « 

^bjuTtn-S 

O * «)LUO 

KhMlu - 

• 4 fun*tfi 


o 

x 


LA w 
CUN X 


in 


3CU — t 
bl->~ S X 
ftwiu - c 
^*•41 rupr 

ru o ru — . 

2-2 

»-» (\) • «w 

UJ3 ~ ~ H 

c wru ru z 

Ui X — w -a 

I -xt-rl 

a cu % «ivjma 

ck *7^*^ * % * 
H • w fUAJfU ^ ~ 
w wio 
hi/tw-hn *w 

r • r -o^mu. 

a ru — 3 «qc w u. 

£*’.£*■ 2 

“5§S2§fg 


O % 

JW 

u 



t-jM mQKm, 
-OOQ O O 

— ru 


ox r -h 
a ^ u k a 

— * H -CL * 

u *oo * 
k la«-»« 


«Sfe“.f 2 . 

3 — ia uj -~ a a. — 
Or w ^oou.® 

• nailed -3 la 
o: 3 uj >-• u. — « ~ 


% w 


«.Ska=S66*H!2«*i 

LA 

S s .£a 

LAN 

a 

X Ui -Zlax -la - -^ *- Ui 

r 


ru 3 

X 

a aw-L. nornotH 


uj h az 

w *-« 

*— « 

LJHO • -9 -CL— <Z <1 LA 

H 

x - 


n 

— o 
la h 
« 

a> <i 

»- 

s 


3 

CL 

r. 


LA 


LA < _ _ 

cut- O 

w r “ 


IftOWH 
ruoH® 
w * <z 

g^?5 


. > 4 M W M M M 


LA NLA 

r lat 


wO~H— ~ 3 «X -U. 

3 w <* (/> u •- cn - 


KOI -® 

a •IhH 

**J -o — - 

OH ~ OiA 
jaoiDa 
at- * * 
ouj <x q ~ 
h -u 6 ® 
Q Aj 3 « — 
»aaz w 
▼ a • —3 
k a -ua 
• q. w * * 
m -~o~ 


iiiniJ ~ 
3o- - <i . j ~ J 9 
ZdhOJ *H LALA 
jcao.-co.w 
•f->- J®Z -H 
r - 3 la — uj ru 3 
►^dliI'' *i a 

or X 2 -nj > la 

—• x - ru — — ui u. 
uu^Ncaa - 
or -• ojuju. -N 

*~iA3 -3 CUH 

g <j»- w- a a -j 
»*k 5® • j r 


— HH l*<Xi®Ot — H «« *JC 

qw 3 dhm«® '-h^r 

Or * C - vA ui ru a ru — u. UJ 

ui m >wonawu.<lNQLO • 


la LA 

«* 


CL 

O 

9 

— ~ H 

U> • * 

CU LA - 
X w w LA 
ON- 

'2*2 _ 
~ X — ^ Dh 


O O 

h or 

3i 


W LA 

8p 


• won 
h -joc •« 

•WHO * 

a* 3 a .j ~ 
h • -U. ® 
Ci^Ou ^ 
*1-10'- 
(Ud. L -CU 
*- -3 « <r 

• h gr ® - 

— r 5 la ~ 

Hu J • — «► 
Q -ON — 

--J i D w 

— Xu. • — 

huiJ^a 

S -UJ % » 

r- -iA ^ 

o a -in - 
odu. -a 
v -a ^ h 
X -»H «> * 
O ui -IA — 

C *C wf- 

x xa n r 
o or o H uj 

o 

— run ▼ 


- h u — a »n 

h ><IJZ -ru a. 

3 ^ •« a n — -j 

Cl ® u - -Q. <r 

H — ui — ^-JQtpL 

3-oa «>ca«3o 
o ru — -j — Or a 
%h oc (A ^ • -ru m 
owaiNruHLAiA 
O <t • »*JCLAiQ.— 
u *— — u j O ^ 
H ^ N -y> u & — 

• — h — j ® »a -la 

-WOW— CU * O UJ 

(wcuaa^o — rua 
— -»u *n <r <l <r 

5^ - « E & ja b 

« U- -&d.<LCA 

• — 

— «cpoc — — — 

C— w • — <V — 0-H 

--no-a« — J cn 

LA tor- ODOM 


5 

X 

o 

t 
•— « 
CL 

o 

LA 

X 


QC 

2 

2 


0 

-J 

1 

X 


Di 

a 

0 

1 


LA ac 

3 


o 

“3 

3 


w o 

S 3 


it 

0 * LA 

1 O % 

a « h 

»1 s 

UJ « • • 3 

5o «o — 

¥ 

- « 
r x 

»H LAO 
UJ O 3H H ►-•X 
H X TLAQCX • 

-or -3 CL UJ LA o 

U!la3x J wS 2 

3 

S XXH -LAW 

-ooio h Z 

(kndHtt 
-XC NLA 

k -XH * 3H 

LA -HlJm - 

«^35=« fe 2 

k?^5 S N 

— a: 3 la i a ^ 

*Ja*- 5 is 

s 

— 3 CL LA C -LA O 

ax o«xi • 

- •flrtdaj a 
ljolx • *h a 
JP « -LA X * * 

raktfis h 

\Oh|/|Ihh S 

\Xh -O n H H S 
X -X 3 -X -LAX 

OOftIUOHM© 
OH ‘ O' -O - O 

— ruc*S — — cu 


a 

a 

a 

x 

% 

o 

LA 

« 

O 

3 

O 

6 


o o o o o o uooooo ooo 

- 5S; “ 5? " s 22222£2S22iBS8SSjc»SaRSS«SS««fcifl?^5i«n 


:,- 49 


$t& 4 . 


Q 

u 


»- 

i 

O 


i 8t 

*o «y 
►« ojoc 
h * 3»- 



X 

a 


o 

3 


ilatl 

S w «o 
% %«»- • 
oc 

-ms 

_ 

\-S&R ' j *J 


_ lowruduio 
» 4ty *« »Johh 

*- o< 5 £ - 


J5 «T1C3 

5^-ujo:d± * ~ 

Ear HOj« 0 (^K\t 


5 

CL 


o 

X 

at 


a 

c 


o • 

- OJ 


r 

3 


•S 3 

cd % at 
u. -I • 

. %u i * 

Ss £ 

r\l * x 

«Q • 
X l/> 3 

in * • 


x r u. 

^ u v * 

- .TX 
. % .(VI 
UJ • ^* * 

§ «• 
m • r- 
mac u. ^ 


r uj n 

►- qc m 

X 


>03 
X • 

<e r- 


*x-J r * 
x r T • • 
to r- 

% ■» 

' DX * 
k r j*xm 
u ' iOJI 


OJ 

ft • 

OJ VO 


<4 • 

. <L 


14 X * 

U <4 CO 


U. <t QC 

i h o 


“a 

N 

in x 


co 

ftUi 

n ■ 

• 14 

‘ft* 

WHQ. 


ft X 


X 

(D 


fi 


<n«4 


o 

X _ 

▼ 

5o* 

ftUJ 

n • 


o«n 

o 


S * i; :i 5 .& 


X 

OJ 


X 

OJ 


o ▼ c* ►- 
z «uMn 

Qt CO ■ * 

•o-* o 


OJ 

CM 


or- • 

u 


fS 


3« £-gf Sx» 

5 mu » C t « 
H OX 
ft xo •<* UJ 
u. X r X © • ' 

r 0) +ot $r * 

5 

• • wO 

U. 9 UJH 


58 C: 

5<r*fe X 

• • • «t 

:?i i-: 

uox x r 

t * ftift. 

033 3 * 


Ml4 5 
*u * 
u. *n 


U. 9 « 
ft' ftOX 

x o • * oi 

▼ X X * 

• or oj ▼ 

-» » • • 

x • ^ 9 

O X QC • UJ N 

d co a. 0 ) • 

c »xu • 

ft ft r* * x 

• a. %x 
xr xoj 

▼ LJ 0"» • 

• P 9 W ?- - - - 

ft * v -in co c* W ; *9 

> «TU)Z u-o qc oj • Q. 

m>-x *ur 


i l 

:x i*U 


r *■ 9 

• • UJ 

? 9 * 

<3 • 


i « ft *9 

T? S'*. 

9 « ftX 

If = “ 


•T JjE 


XU 



r imdci * 

;:** S A8.*H?5 


OX X 

~ r; qV Jj > ^ 

~ in x. uj cl » «> • 
-ft 9 - ~ <t 3-J *2 


■ ■ 

uj -> u. cd r- o 
r ^ t u u. uj 

9 m i • QC 


- - * 

Ifc x 

d ▼ u. 

O 

£ 


“£ 


fttii 
-?uF- 


• x • 

^ _ w o ao x 



5 9 
* UJ 

n • 

CD O) 
U. N 
• U. 
30 


s 

5* 


**.3 

* 9 d 

• uj r- 

OM « 
030 

0 . • 

o 

TO 

% 


9 UJ O 

UJ • x 


3 0 

*11 

x xo 
r wo 
* « • 
OJWO 
. • UJ 
CD (ft. QC 


sHiltiiilppiipi ppiis i 

Sio8S8omm0o3mow3u.! 5 *^, u ^,rvj n + 

S P s ■ 


m(v r* m »- K w 

x x oc o u oc aooj 

rtnzu. *3 !li 




oouuo 


o oooo 

I-fur) ▼m 


»-D0 


*01 

XI^U 

5 v* 

CX • 

j 

«k< 

o «« 


QO • 

•on 

u •Ui 

ax % 
•— f\l ' 

« « • 
ov 

X '3 


I 


Xh 

OJ' 


•MX 

Q • n 
CL (A • 

• U *• 

* *t 

O • • 
K *01 

x«aw 

•k w 

§ ‘1 

X 

-QftJ 

5x • 

(f V • 
£ 1 *% 

sf # 


3« 

'Ui 

n • 

OKh 

V. 

Mj 

xc*x 

• u * 

8o$ 

“i*. 

& u 

xfc? 

** 
▼ • % 
•AW 
0>Ui€ 

2-1 

tfUs* 

*S8 

• • • 

"!•¥ 

<ir x 
u u • 

• • Ui 

oo 5 

XL- ~ 

*i* 

h *w 

|*^ 

8CS 

k • a 


X 

(U 

<x • 

h » 

Ui • — 
o»r^ ru 
•u • 
« *r 

K XU 
Jf % 

ssv 
% • 

ujr- c 

■hfc 

*5“ 


* 


if 

x * 

X 
U* 
•II % 
U • ▼ 

• *u 

iT .* 

• A 

I«g 

• M X 
ui *r 
a vo • 

•u x 

Ui ««• 
UiX • 

fc ru - 

• *x 

xdv 

•k • 

S *2 

^UiOl 

Sx * 

•t X 
(|w(U 
Vh % 



x 

X 

(JO 



X 



X 



a 

? 

«, 





u 








X 

Of 

o 

% 

• 

▼ 


% 

V 


• f 


X 



»- 




ro 

n 

X • 

VOOV 

X 

VO 

X 

ru 

r> 


n z • 
• id at 


«« 

• 

3 



X 

% 

5 


U) 

VX) 

•Ui 

a 

% 

Ui 


03 Ui 



* 


lj 

5 

Ui 

u. 

u 

r> • 

r> 


■ 


u • 


J5 

n 




% 

a 

• • 

• • 

• « 

• 

• 

•# 


• <x 


• 



-i 

r 

? 

* 

•r 

VDK 

03 

r- 



U L- 


% 

VP 


K 


X 


• 

U. -J 

u 

u. 

— » 


Ui 



U. 


X 

X 

v 4 

5 

2 « 

2^ 


• 

• 

UI 


V 3 X 

o 

• ^ 

■ 


ui 

•4 

« 

3 

X 

o a 

»> 

m 

C* 03 

Qc 


— J 


X 

CL 

• u 

■ u 





L- 


“D 

Ui • 



2 

% 

(A 

a 



O 


• 


z 

</■ 

w V 0 

O 


D 

• 


J 9 

hO 

X 



X 

X 

ru u 



3 


$ 




00 vo 

• 

00 


• 00 

% 

». • 


X 

00 

U 

U 

X • 

« X 

r* 

• 

• 

•H % 

3 

X. 

X ^ 

ru 

K 

U. 



°i 

a i 

VO -T h vo 

u 

^ VO 

* 

Ui 

ru * 

VO ▼ 

• 

• 


•ru 


• Ui A • 

cx • -o» 
0 Ui 30 

“§1* 

a .o 

sth 

•t • • 
X • u. 

•013 • 
AWL Ol 
Ui • *Ui 
• Ui 3 • 

Ui X ui 

sfss 

£5 5 * 


to t" or ifl 
u u w 
■ • * • 
WWQO 

EL*3 

5x*.5 

W ft VO ~ 


g * 


fi 

r> 

vo 

u 


• =5 VO 
<*CU 

ujy • 
• -.j 
3 vox 
o \ 

3TL 

*& 

X 


? 

3 

* 

CL 

O 

K 

s 

CO 



• *<A 

o ^ 

• • 

T£V 

'V 0 

CA • • 

XU 

UVOCA 

*4 • 

■ d(ui 

• • 

J| • 

▼ • 

C M 

OI M 

X .Ui 

a 

X 

ui 

• _i 


i? X 

-i • 

oo ui r^ 


• a • 

X X 

« •« 

• 

• vO • 


ca 

ui 


o 

X 

a 


x ~ 

ifi 


$2 

ru 

{ x* 
4 r\J 




>uico 0 * 

.2U U 


3^Ui 

h-D 

r- 

cax 

x 'vr 

vo vo w 
*w»- 

M I 

(X 


r* 

u 


x 

K 

vo 

i 

X 


VO v/) 
u o 
• T 
X • 
fXJ * 


VO X • 

5 *£ 


r* -• 

u \ 

*n 
x • 

' 00 


or 
% 
u 

VA 
K 

? 

VO 

a 

i 

vr 

Ui 

a 

6 

vo 

w Kh 

P I* 

5 £* 


« 


i 


at 

Ui 

■ 

ct 


r 


- -5 x 

•* r* • F 


r* ■ 
u 

•vo 
ft 3 


o» X 
Ui lO 


z • 

ru 


r* • 
u 

<i2 

Ui • 

•* 


or 

O 


o 

s 

£ 


u 


01 

Ui 

ru 


x 

ru 


vo 

u 

5 


gS\' 

5 u,*t 5 

U. w #- • 

in*. 

- -fir 
utOu 
~3U. • 

M « 

s 


n £s 

* Ui~* 


X 

nor 


a « 

r 

d? 

JSf 


£5 

d£ 


% % 

•UI N 

xJn 

'i*- ’ 

r < 3oJ 


X • 
00 r) 


oi • oi 

x Jn ui 

O^VJ-l • 


#--j <r 


•a 

o 


U X 
-*£*• 
c 


~J • h 
W h- W 

g >f 


- a ru — ux^-j • vo m 

U c -k8h*^ 

u.5 

w J ff* 


•»Oh X 

nS.JiSr, 


u 


a 

$ 


-j vo 

5 ., 

-5 


CAW 

Ui Ui 

*5 


u 

a 


?*•*«»„* tfn 

•r»a .XX • X • 
— vt 

Lm 

fi* 

un 

• • *o • • 

^ 2 |ST-* 

V 0«4 •XhAU 

« .u UNJ »Vh • 

« —i 1 (u r> ♦ m 
XkWhO • • V * 
••k winuioB 
_ ^ • ~ 4 JVI <3 

o \ \ «VJ V < 

b w . w »u» 3 v ja H 

i- ••*-«* w »r < 

w & • wK • vft ••— UI U 

kOAkO 0 kOXftaOXX X 
MkkMkk^kttllXkQ A A 

— •< M « M « M M 

s 8 * i 


3d 

vr 


v> ru 
U.H 


h 3 
vl 


OUUUO UUUOOU u ooououo u o uoo u uuooooooo 


B-i,i 



SUBROUTINE NEUT (A, b , C, T. ICANT) 


in 





Ui 





a 




% 

3 



o 

• 

h 



<z 

3 

o <r 



3 

X 

r cl 



u 

cu 

M UJ 





q a 



K 

% 

r c 



3 

m 

H-» UJ 



in 

• 

3 K 




cn 




• 

UI 

O' c 

CO 


m 


O 3 



X 

% 

u. •-« 

i 


o 

• 

CL 



*— « 

<r 

O ou 

m 


K 

X 

O ►-« 

M 


a. 

t\j 

X 

H 


CL 


1- 3 

fr- 


UJ 

% ^ 

UJ *”-* 



»- 

x n 

E 3 

H 




O 




in 

X UJ 

<X 


> 

* 3 

o 



»~ 

3 

in a 




O % 

x o 



•— » 

<r * 

CL 

• 


3 

OQ 1- 

a 3 




3 <r 

nL 3 

1 X 


X 

X fr- 

1 <x 

W CL • 



<3 3 

x 3 

3 ru 



3 

o 

\ h- 


o 

x n 

fr- Ui 

~ 3 N 


3 

O X 

3 H 

O C* ~ 


X 

in r- 

UJ 3 

fr- 


M 

X 

r a 

♦ ~ <r 


<r 

a ■> 

x 

— fr- 


fr- 

<r n 

3 O 

»~ 0 3 


ee* 

3 • 

X O 

<S> 3 


O 

fr- cn 

fr- 

* • Q 



X X 3 

o 



7~ 

fr- O 

in k 

CO • 4 


3 

3 fr— % 

> 

fr- 


3 

3 3 • 

o in 

1 3 fr~ 


no 

Q U h 

3 fr- 

• 



* X X 

CL O 



»- 

h inj 

E O 

*> ^ « 


O 

* ru 

u a 

O i J- 


X 

O OJ * 


IDh h- 



- D 

UI 3 

\ 


in 

00 - • 

x a 

* M H ^ ^ 


<i. 

• x cr> 

M M 

^ <r a * 


X 

<x cn 3 

h* E 

a fr- H <s> 




3 O 

• 1 3 3 


3 • 

O * - 

o x 

^ 3 3 • 


O X 

<s> »• 

CL > 

<$ x nouj 


x a. 

fUIO H 

00 3 

• ^ J UJ 

3 3 

- -- X 

3 O 

» i in • 

3 

O fr~ 

3 ~ i\J • 

in a 

® h 00 

X 

Q: 3 

w h~ z. 


h®<Ih^h 

t— « 

3 3 

3 a - or 

in 3 

Z ^ h w w 

fr- 

O 

fr- E m x o 

^ a 

<1 3 ■ 

X 

X o 

Q • 'l h Q 

X 3 

0 0 3 U U 

o 

O X 

u o cn o uj x 

•- a 

Q Q fr— 

o 

O <1 

j u uj h a uj 



o 


c> 



<s> 


<s» 



•-4 


ru 



oooooo oooooo 

• • • •••»« •• »• •• M • • •• •• •• •• •• •• •• t 9 9# M Af • # • « ## ## # * 

Hf\jnv Ln«x>r-oocD<s>-*rurn^ inior-oo(7><s>-«rj' s >^io3 

«h «-« «-4 m •-< •— « •-< r\j ru ."\j nj f\j nj 



SUBROUT I lit OPTMYZ<H> 


II 
2iC 
3«C 

4 i C THIS SUBROUTINE PERFORMS THE OPTIMIZATION PROCEDURE. 

5*C 

SIC 

?l COMMON/ THICK/ DTEQ, TDOT, TRE, QC. OR. ENC , EMISS. MRECOU. 

8l 1 DT # TINE. PL. DTO, DTI. TEQ2 

9i COHHON/OPTT/ TOPT.PERCNT.DELO.TI, 

19* 1 SIGMA, tniS. NATL. NMTL, RMOM, CPN. ICAlST. I DI C 

111 COMMON/ FLUFLD' PHIL( S ), PE. PHILISU9). TU. PHILUPUS) 

12» DIMENSION TVMi 100 > . E NCS( 100 > , HRS C10O). TS< 100 ) .PS ( 109 ) , 

131 1 CPMAT<10>, TNAT(10) 

14*C1883 FORMAT ( 1M1.54X.2HBD, 14, SH RR.I3.7H CASE , I 4//20X . 
lSlC 1 39HAE RO D YNAM I C HEATING £ NUIRONMENT —ARY / 

lStC 29X.4MTIMC.4X, 13HHCAT TRANSFER . 4X, SMRti r RV . 6 X . SHLOCAl . 10/.4HWALL/ 
1?IC 317X. IlHC0EFFICIENT.6x.8HENTHALPY.SX.8HHh uFT , SX, 1 IMTEfIPE SATuRF / 

18 1 C 48X. SH< SEC ).4X. I 3H( LBN/FTP -SEC ),4X,9M<BTU u&M ) . 4X , 9M( LBF /F T2 ) . 6X . 
1»IC 57H( DEG F ) /) 

281 2888 F ORMAT ( 1H0, 49X. 29H0PT IN I ZAT I ON PROCEDURE RESULT/ 1H8. 8x. 

21i 1 1 SHAT ITERATION. 13. 

22 » 2 3SH MAXIMUM UALL TEMPERATURE REACHED 

23 1 3 F 11.3, 34H DEGPt£S F WITH A SKIN THICKNESS •. Ell. 5. 

24i 4 8M INCHES./ IH . 8X 27MC0MPLETE TRAJFCTORV UILL BE, 

25* 5 SHH RERUN USING THIS THICKNESS.) 

28* 2588 F ORMAT ( 1H0. 7X. 36HTRAJECT JRY UILL BE RERUN WITH A SKIN, 

27* 1 12H THICKNESS •, Ell.S) 

28* 3000 FORMAT* 1H0, 18X, 37HFAILURE TO OBTAIN CONVERGENCE OF WALL. 

29* l 47H AND OPTIMUM TEMPERATURES AFTER 10 TRIES. CASE, 

39* 2 11H ABANDONED. ) 

31* 4000 FORMAT ( IH0, 5X . 22H0PT I MUM TEMPERATURE OF. Ell. 5. 

32* l 5 1 H DECREES F IS GREATF R THAN THE MAXIMUM LOUILIBRIUM. 

33* 2 t SH TEMPERATURE OF, Lll.S.ltH DECREES F . / IN . 5X. 

34* 3 16H CASE ABANDONED. ) 

35* 5000 FORMAT* IH0, 7X, 23NUITH A SKIN THICKNESS •. Ell.S, 8H INCHES.. 

36* 1 3SH MAXIMUM UALL Tf MPERATURE REACHED F10.3, 

37* 2 1 0H DEGREES F/ IH . 7X, 1 5HUHICH IS WITHIN. F103. 

38* 3 36H PERCENT OF THE OPTIMUM TEMPERATURE.) 

39* 6000 F ORMAT * 1H0, 48X, 32H0PT ]M 1 ZAT I ON - SURPRISING RESULT/ 

40* 1 1H0, 13X. 29HUITH INITIAL SKIN THICKNESS -. Ell. 5. 

411 2 43H INCHES. MAXIMUM UALL TEMPERATURE REACHED •, Ell. 4. 

42* 3 11M DEGREES F./ ]H . 12X. J6H WHICH IS UITHIN, F7.3, 

43* 4 51 H PERCENT Of THE OPTIMUM TEMPERATURE. v<)UR COMPUTER. 

44* 5 43H CONGRATULATES YOU ON YOUR EXCELLENT GUESS.) 

45* 7000 FORMAT* 1H0, 50X. 27H0PT INI ZATION RESULT - RERUN. 13) 

46* 8000 FORMAT* 1M0. 7X. 39M0PT IN I ZATION PROCEDURE WILL BE REPEATED. 

47* l 24H TO COMPUTE A N£ U GUESS. ) 

48* 9000 FORMAT* 1H0. SIX. 25MFINAL OPTIMIZATION RESULT' 

49* 9500 FORMAT* IH . 7X, 44NV0UR FRIENDLY NEIGHBOR^ 0 TER SUGGESTS/ 

58* 1 IH . 12X. 43H( 1 ) LOOSENING VCUR CONUERGENCE CRITERIA, OR/ 

51* 2 IH , 12X , 4 3H* 2 ) IMPROVING. BY WHATEVER MEANS AVAILABLE,. 

52* 3 27H YOUR ORIGINAL LOUSY GUESS.) 

53 ! C 

54 x C THE FOLLOUING EXECUTABLE STATEMENT T(>pT. TOPT * 459.7 

55 * C UAS Rf MOOED FROM THE MAIN PROGRAM < H8M0 > TQ THIS SUBROUTINE 

56 » C 

57 i TOPT- TOPT ♦ 459.7 

58* IF (M - 1> 100. 200, 300 

59 *C 

60tC 

61 *C SET THE COUNTER. ERROR FLAG, AND FQUILl IB»IUN TEMPERATURE TO 0. 

62 *C 

63 *C 

4i 18» :♦ - 8 
8* TEONA* • 9. 

881 ICANT - 0 

•7i RETURN 

88 *C 

69 * C 

79 *C SAVE FLOUFIEID PARAMETERS AND STORE NEW EOUILLBRIUM TEMPERATURE 

71 »C IF IT EXCFEDS PRCUIOUS MAXIMUM. 

78«C 
73*C 

74* 299 N • N ♦ S 

75* TVM(N) - TIME 

76 i ENCS(N) • LNC 

77* MRSfM) - MRECOU 

78* TS(N) • TU 

79* PS ( N ) • PE 

89* IF * TRE .GT. TFOMAX ) TfQMAX • TRE 

81* RETURN 

82 * C 
83*C 




' /> 

% 2 Q *- 


() n 


•4«C IF THE OPTIMUM TEMPERATURE IS LESS THAN THE MAXIMUM 

9SlC EQUILLIIRIUM TEMPERATURE, ABANDON TRV. 

•SIC 

•7IC 

Ml 3## NMAX-N 

•91 IF (M.EQ. 4 ) CO TO 1000 

Ml IF (TOPT . GE. TEQMAX ) GO TO 950 

9|tC 

oaic 

93 i C DETERMINE MAXIMUM UALL TEMPERATURE REACHED. 

94 1 C 
95 1 C 

Ml TMAX • 0. 

971 DO 40# N • t » MM AX 

9Ri IF (TS(M) .CT. TMAX) TMAX • TS(N) 

991 4M CONTINUE 
lMt IF (M .EQ. 3) GO TO 850 

1011 I TRV • 0 

1021 ICANT • 0 

1031 DELHEU • DELO 

1041 5M TDIFF - (TMAX - TOPT)/ TOPT 

1051 TDIF2 - (TMAX - TOPT )/ C TOPT-TI ) 

IMiC 

107IC 

10«»C CHECK IF COMUERCENCE HAS BEEN OBTAINED UITH THIS THICKNESS. 

IRfttC 

110IC 

111* IF (ABS(TDIFF) .LT. PERCNT) CO TO 700 

112» IF ( I DID .EQ. 2) CO TO 750 

113> IF ( IDID .EQ. 3) GO TO 800 

1141 ITRY - ITRY ♦ 1 

1151 IF (ITRY .CT. 25) GO TO 900 

USiC 

1 17 *C 

1 18 * C COMPUTE NEU THICKNESS AND CO THROUGH SIMULATED TRAJECTORY TO 

119IC GET NEU MAXIMUM UALL TEMPERATURE. 

120 1 C 

i2i>c 

1221 DELNEU • DELNEU * TDIF2 ♦ DELNEU 

1231 Tg - T I 

1241 TMAX • TU 

125* NM - NMAX - 1 

1261 550 DO 600 N • 1, NM 

127* TUO • TU 

lit* HU - 0.24 t TU 

1291 QM • ENCS(N) t (HRS ( N ) - HU) * SIGMA X EMIS X 7UIX4 

1301 CALL MATRESCHATL, QN, RMOM, CPM, DELNEU, TDOT, TU, THAT, CPMAT, 

1311 1 NMTL) 

1321 TU - TDOT f <TVH(MM> - TVM(N)) ♦ TuO 

1331 IF < TU .CT. TMAX) TMAX • TU 

1341 600 CONTINUE 

135* GO r 0 500 

1301 700 IF (ITRY .EQ. 0 ) GO TO 860 

1371 TMUR • TMAX - 459.7 

1301 URITE (6,2000 > ITRY, TMUR, DELNEU 

1391 DELO • DELNEU 

1401 RETURN 

1411 750 T1 • TMAX 

1421 D1 • DELNEU 

1431 DELO • DELO ♦ TDIF2X DELO 

1441 GO TO 82S 

1451 800 T2 - TMAX 

1461 D2 • DELNEU 

1471 DELO - 02 ♦ (TOPT - T2 ) f ( D2 - D1 )/ ( T2 - Ti ) 

1481 825 URITE (6,2500 ) DELO 

149> DELNEU- DEL 0 

1501 RETURN 

tsne 

152IC 

153 iC If THIS UAS RERUN OF ORIGINAL TRAJECTORY PRINT RESULTS AND 

154 l C RETURN 

156 iC 
ISSiC 


tS7t 8S0 TDIFF • ABSCTNAX - TOPT)/ TORT 
1581 TDCENT - TDIFF t 180. 

trot fun . thaw _ 4 CQ 7 

1601 IF (TDIFF . LE . PERCNT .OR. IDID .GE. 3) GO TO 855 

1611 UR l TE ( 6 • 7000 ) IDID 

162t WRITE (6, 5000 ) DELNEW, TflUR, TDCENT 

1631 WR l TE ( 6 , 8000 ) 

1641 ICANT - 1 

165i RETURN 

1661 855 I CANT • 0 

1671 URI TE (8.9000 ) 

1681 WRITE (6, 5080 ) DELNEW, TNWR, TDCENT 

1601 IF (TDIFF .GT. PERCNT) UR ITC ( 6. W500 ) 

1701 RETURN 

1711 860 TDCENT • AB5( TDIFF ) t 100. 

1721 TNUR - TNAX - 459.7 

1731 WRITE (6, 6000 > DELO. TNWR, TDCENT 

1741 ICANT • -1 

1751 RETURN 

176IC 

177IC 

178 1 C IF CONVERGENCE HAS NOT BEEN OBTAINED IN 25 TRIES ABANDON CASE. 

1791C 

180IC 

18H 900 WRITE ( 6 , 3000) 

182i ICANT • 1 

1831 RETURN 

184 (C 
185 1 C 

186 1 C IF OPTinUfl TENPERATURI IS GREATER THAN NAXinun EQUILLIBRIUN 

187i C TEMPERATURE, ABANDON CASE. 

188iC 

1891C 

190) 950 TOUR • TOPT - 459.7 

19H TEQUR • TEQNAX - 459.7 

'Iff* WRITE ( 8 , 4000 ) TOWR, TEOUR 

iWl ICANT - 1 

l|4i RETURN 

190IC WRITE OUT ENUIRONHENTAL SUNNARV TABLE 

197IC 

198H000 CONTINUE 

109 IC WRITE (6, 1003 ) < WU J >, J-l ,3 ) 

0001 DO 1002 K*1,NNAX 

0011 TWF.TS(K)-459.7 

8021 WRITE ( 6 , 1001 ) TYN<K ), ENCSOC J,HRS( K ), PS<>C>. TWF 

8031 1001 FOR»»AT<6X,F8.2,3X,Ell.S,6X,F8.2,3X,Eli.S,4X,F8.2) 

0041 1002 CONTINUE 
206i WRITE (6,1 004 ) 

206 » 1004 FORMAT! 1M1 ) 

2071 RETURN 

END 


S 


PCSU0601 


11 SUWOUTIHE PCSU<X,Y, ITABLE, OFFTBL,Z) 

3IC 

3iC 

4iC 

5lC THIS SUBROUTINE CONPUTES PC AND S.U. ANCLE FOR A REAL CAS 
6iC 

7tc mtiiiitiimuiiiimmmmnmiiiiiiiiitiiimim 

8*C 


9* 


101 


ill 


12t 


13t 


141 


15t 


16i 


171 


I8i 


19 1 


201 


21t 

16 

22t 


23t 


24t 


25i 


26t 


27i 


28t 


29i 

30 

301 


3H 


32i 


33i 


341 

50 

3St 


36! 


37t 


38! 


39* 


401 


411 


421 


431 


441 

76 

451 


46t 


471 


48i 


49 1 


56 * 


51 1 


521 

100 

531 


541 


55f 


56 1 


57i 


681 

101 

69i 


00t 

102 

6li 


62i 

163 

63i 


64i 

104 

65 > 


66* 

105 

67t 


66* 

166 

69t 


761 

107 

711 


72i 

700 

73i 


741 


761 


76i 


771 



COAAON/ERRFLC/NERROR 
DIMENSION A ( 6. 6, 28 ) , BC10O8 ) 

EQUIVALENCE (A< 1,1*1), Bd)) 

OFFTBL-0. 

X3-X 

IFCX. GT. 26. )X-26. 

IPATCH-0 

IF ( X .LT. 1.) CO TO see 
IF ( V . LE. 0. ) CO TO 516 

IF< Y.GT.60. .AND. ( ( ITABLE .EO. 2 ) .OR . ( ITABLE .EQ.4)))G0 TO 526 

IF ( Y.GT.S5. .AND. ((ITABLE .EO. 1 ). OR. ( ITABLE .EO. 3 ) ) )GO TO 

CO TO (10,36,56,76), ITABLE 

IF( (X. LE • 1.7). AND • ( V • LE . 1 6 . ) )GO TO 101 

IF ( (X.GT. 1.7). AND • ( X . LE . 2 . 8 ) . AND . ( V . LE . 32.))G0 TO 102 

IF (X.LE.2.8) CO TO S20 

IFCY. LE. 16.) CO TO 163 
IF( V. LE . 35. ) CO TO 104 
IF( (X.GT. 3.4) .AND. (Y.LE. 

IF (X .GT. 8.) CO TO 106 
CO TO 526 

IF( (X.LE. 1.5). AND. (V.LE. 


45. ) )GO TO 165 


28. )>GO TO 161 


IF( (X.LE .2.8 )• AND. (X.GT. 1 .5 ). AND. ( V.LE .48. ) )GO TO 102 
IFCX.LE. 2 .8 )GO TO 520 
IF< Y.LE . 35. )CO TO 103 
GO TO 164 

IFCX.LE. 1 .05 )GO TO 500 

IT( (X.LE. 1.45). AND. ( Y. LE .8 . ) )GO TO 101 
IF( (X. LE 2.0) .AND. (X.GT. 1 .45 ) .AND. (Y.LE. 20. ) )GO TO 102 
IFCX.LE . 2.0) GO TO 520 
IF( (X.LE. 3. 5). AND. (V .LE. 12.))GO TO 
IF ( (X.LE. 3. 5). AND. (V.LE. 36.)) GO TO 

IFCX. LE .3.5) GO TO 520 
IF(Y.LE. 26. )GO TO 105 

IFCY. LE. 35. >GO TO 106 
GO TO 167 

IF ((X.LE. 1.5). AND. (Y.LE. 32. )) GO TO 
IFCX.LE. 1.5) GO TO 520 
IF ( ( V . LE . 20 • ) . AND .(X.LE.2.8)) GO TO 102 
IFUV. LE. 48. ). AND. (X.LE.2.8)) CO TO 103 

IFCX. LE. 2.8) CO TO 520 

IFCY. LE. 24. ) GO TO 104 
IF (Y.LE. 40. ) GO TO 105 
GO TO 106 

C-<X*X-t )tS INC .0174533BY) 

C-1.4 

Cl - ( ( G+ 1 )SC+2)/(CG-l >*C*2> 

C2-AL0G(1.2fi./C) 

Z-d+Cl*C2 )*(SIN( .0174533*V) )tt2 
RETURN 


103 

104 


100 


PCSU0615 

PCSU0026 

PCSU002S 

PCSU0030 

PCSU0035 

520PCSU0640 

PCSU00S0 
PCSU0055 
PCSU0060 
PCSU0065 
PCSU0070 
PCSU0075 
PCSU00E0 
PCSU0085 
PCSU0696 
PCSU0695 
PCSU01 00 
PCSU6105 
PCSU0110 
PCSU01 15 
PCSU6120 
PCSU6125 
PCSU0130 
PCSU0135 
PCSU0140 
PCSU0145 
PCSU0150 
PCSU0155 
PCSU6160 

PC5U0170 

PCSU0175 

PCSU0180 

PCSU0185 

PCSU0190 

PCSU6195 

PCSU6200 


IPATCH-1 
GO TO 70# 

IPATCM-2 
CO TO 700 
IPATCH-3 
GO TO 700 
1PATCH-4 
GO TO 700 
IPATCH-5 
GO TO 700 
IPATCM-6 
GO TO 700 
IPATCH-? 

IF (X.GT .20. ) X-20 
CONTINUE 

IPATAB • 7 t (ITABLE -1) ♦ IPATCH 
Z-6. 

Xl-l./X 
DO 740 1-1,6 

IF (A(l, I, IPATAB) .EQ. 6.) GO TO 760 


PCSUO205 

PCSU0210 

PCSU6215 

PCSU6220 

P C5U6225 

PCSU6230 

PCSU0235 

PCSU6240 

PCSU0245 

PCSU6250 

PCSU0255 

PCSU0260 

PCSU0265 


PCS J0275 
PCSU0280 
PCSU0285 


m- :> r 


PCSU0295 


0*f$ it p 


70* 
79* 
89* 
81 1 
82* 
83 1 
84* 
8St 
86 * 
87* 
88 * 
89* 
98* 
91* 
92* 
93* 
94* 
95* 
96* 
97* 
98* 
99* 
lift 
191* 
192t 
193* 
1941 
195* 
196> 
1871 
198* 
199* 
1191 
Hit 
112t 
1131 
1141 
1151 


116* 
117* 
118* 
119t 
129* 
121 1 
122 * 
123» 
124* 
125t 
126* 
127 1 
128* 
129* 
139* 
131* 
132t 
133* 
134t 
135* 
136t 
137* 
138* 
1391 
149* 
141* 
142* 
143* 
144* 
145* 
146* 
147* 
148* 
149* 
159* 
151* 
152* 
153* 
154t 


DO 739 J-1,6 

IF ( A( J, I, IPATAi ) .EO. 0.) CO TO 735 
2 • 2 ♦ A( J, I # IPATA8 ) * (Xlttd-m * (Vtt(J-l)) 

739 CONTINUE 
735 CONTINUE 
749 CONTINUE 
769 X-X3 

IF < C ( ITABL6 .60. 1 ) .OR. ( ITABLE .6Q.2 ) ) . AND. (Z.GE.90. ) )0FFTIL*2 
IF( ( ITAilE •EO. 3 .OR. ITABLE .EQ.4).AND.Z . CE . 1 .8 )0FFTIL-2 
RETURN 

599 OFFTiL-1. 

RETURN 

519 UR1TE (6,511 ) V 

511 F0RNAT(53H0 ftttt UEDGE OR CONE ANGLE NEGATIUE OR ZERO, ANGLE-, 
1E12.6,9H DEGREES. ) 

NERROR-1 

RETURN 

529 OFFTBL-2. 

X-X3 

RETURN 

DATA (1(1), I-l, 36) / - .455022956+3, .41710586+2. 

1 - . 25260026+1 , 3t0., .19546736+4, - . 14418886+3, 

2. 11371636+2,3*0. , - .2612456+4, . 15487736+3, -. 17671716+2,3*0. . 

3 . 12012416+4, - . 46531366+2, . 967671 EM , 15»0. / 

DATA (i(I), 1-37, 72)/ 

4 . 15622486+2, - . 36101596+1 , -.39188205E 1, -.24032116-2. 

5 2*0., - . 4861 1566+2, .32985646+2, -.3057881. 

6 .4272396E-1. 2*9., .23161276+3. - . 79822076+2. 

7 . 18081336+1, -.1520378, 2*0., - . 15352256+3, 

8 .63321716+2, - . 21045456+1 . .16028342. 14*0./ 

DATA (id), 1-73.108)/ 

9 .2105304E-1. . 1 1458246+1 . 1576662E-2. .4145022E-3, . 1464391E-3, 

1 .11686756-4, .53774586+2, - .6947946E+1 , . 16045016+1 ,-. 33449885 , 

2 . 2550358E- 1 , - .58811416-3, .22928466+2, .22100306+1, .6336085E-2, 

3 .73411788 . - . 7901796E-1 , . 20228376-2. -. 28333216+2, .30640776+2, 

4- • 48633066+1 , - . 10267036+1 , .14236532 , - . 4066452E-2, 12*9 . / 

DATA (i(I), 1-199,144)/ . 12642S3E+3. 1376618E+2, 

1 .6486186 , - . 1357706E-1 . .121729E-3 , .0 ,-.1220996+4 , 

2 .81416116+2, - .26783191 ,-.3160076E-l, .1795490E-4, .0 

3 .54199496+4, -.9224745E+2,-. 17773076+2, .35320902 , .3728U6E-2, 

4 .0 , -. 13022856+5, .30718976+3, .75589366+1. .16417966+1. 

5- .5539958E- 1 , .9 , . 1 368806E+S, - . 99499716+3. .8745465642, 

6- .60934416+1 , • 1231502, .9,6*0./ 

DATA (i(I), 1-145, 180)/ 

7 . 37671996+2, - . 3462785E+1 , .13918317 , - . 1230641E-2, 2*.0 

8 . 32089776+3, . 3631055E+2. -. 17861 106+1 . .1748266E-1. 2*. 9 

9 . 3019779E+4, - . 8006098E + 3, . 28700866+2. -. 25491071 , 2*.0 

1 .59647136+5, - . 3S07169E+4, . 79661 756+2, -. 75989497 ,14*.0/ 

DATA (B(I), 1-181, 216)/ 


2 

331.0 , -.22767006 +4, 

4- • 15731 35E +6, -. 34145156+3, 

5- • 12288706+6, .3798787E+3, 3t.0 
6 .638S858E+4, 9*0./ 

DATA (i(I), 1-253,288)/ 


. 5636795E+2, - . 482661 6E + 1 , 
. 3522476E +3.-.6491 722E ♦ 1 , 
.81019996+2, 3* .0 


.10074503 . 

3*.0 

.48485756+7, 


. - . 27780846+7, - . 23230706+6, 


.149132436+3 

.386450466+4 

2S.0 

.60628930 
. 101993406+2 


DATA (i(I), 1-289,324)/ 


- . 38690958E+1 
-.570543526+3 
-.488114906+4 

2 * .0 

-.31032062 


6 -.469351756+3 

7 .189994266-3 

8 -.395097976+1 

9 - . 1 1617541E+5 

1 .756084946-2 

2 - • 163000936+2 

3 -.743599306+4 

4 • 10012938E-1 
DATA 

5 .96082449 

6 . 45937638E+2 , 
7.377231316-4,0.. 

8 -.580373826-1 , 

9 - . 1 1717956E+2 , 

9-. 197605406+3. 

1 .203719016+2. 


749991316+1 

, .0 

, .538656446-2 

, .208522096+3 

, .0 

, .22614131 

, .123322026+3 

,7*0./ 

(id ), 

, .267373286-1 

-.900428126+1 
.134058826+3 
12922046-2,0 


.786742366-1 

.160572336+2 

.722894986+3 

.209112686+4 

14*0./ 

.62472434 
.390636726+4 
-.179986106-2 
.113726426+2 
.153193426+5 
- . 143596216-1 
.941923766+1 

1-325,360)/ 

-.143712806-2 

.50227665 

.284566536+1 

,-.532710556+3 


-.985437436+3 

2 * . 0 

- . 38328927 
- . 22178436E+2 
-.303387826+3 


-.311417536-2 

-.678842656+2 

.0 

-.572565086-1 
- . 26893756E+3 
.0 

-.23014992 

.11409437 
.19988716-4, 
-.846836186-2 
.31682443 
. 10812230E+3 


. 56917652 , - . 936372686-2,0 . , .718740486+3, 
.89913739, .141461786-1,7*0./ 


0 . 


PCSU0310 

PCSU0315 

PCSU0320 


PCSU0335 

PCSU0340 

PC5U0345 


PCSU0365 


PCSU0415 

PCSU0420 

PCSU0425 

PCSU0430 

PCSU0435 

PCSU0440 

PCSU0445 

PCSU0450 

PCSU0455 

PCSU0460 

PCSU0465 


P CSU0475 

PCSU0480 

PCSU0485 

PCSU049 

PCSU0495 

PCSU0500 

PC5U0505 

PCSU0510 

PCSU0515 

PCSU0520 

PCSU0525 

PC5U0530 

PCSU0535 

PCSJ0540 


PCSU0S50 

PCSU0555 

PCSU0560 

PCSU0565 

PCSU0570 

PCSU0575 

PCSU0580 

PCSU0585 
, PCSU0590 
PCSU0595 
PCSU0600 
PCSU0605 


PCSU0615 


o-57 


ISSt 

DATA (1(1), 

1-361,396)/ 


-.611783226+2, 

.618233766+1 

A 

PCSU6620 

I5«i 

1 -.15594696 

, .196466236-2 


-.728684426-5 

» 

.0 

0 

PCSU6625 

15?i 

2 -.398561456+3 

, .331470916+2 


-.113486486+1 

A 

•24518360E-1 

A 

PCSU8636 

158t 

3 -.258841 136-3 

, .6 


-.491688136+3 

A 

.297952966+3 

A 

PCSU8635 

159i 

4 -.782548826+1 

, .366344636-2 


.169785896-2 

A 

.6 

A 

PCSU8646 

1681 

5 -.243996876+5 

, .619139516+3 


.731811166+1 

A 

.796314236-1 

A 

PCSU6645 

16U 

6 -.63625578 E-2 

, .9 


.444713766+5 

A 

-.818468666+3 

A 

PCSU0650 

162t 

7 -.393124546+2 

, .67556079 


.318719826-2, 

788./ 


PCSU6655 

163t 

DATA (1(1 >, 

1-595,546)/ 


-.3321322 


.17676547 

A 

PC$y#666 

1641 

l .76334979E-1 

, - .288S861S6-1 


. 45590 165E-2 


.0 

A 

PCSU966S 

165t 

2 -.17324872 

, .57856695 


-.42637694 


.523978586-1 

A 

PC SU 06 70 

1S£: 

3 -. 38552S17E-2 

# .6 


.372171196+1 


-.216564696+1 

A 

PCSU0675 

1671 

4 .55374581 

, .411982576-1 


-.109116886-1 


.0 

A 

PCSy9686 

168t 

5 -.583737986+1 

, .144510096+1 


-.11471143 


-.763282196-1 

A 

PC SU 0685 

169t 

6 .26758654E-2 

, .0 


.263247906+1 


.16728179 

A 

PC$y0690 

1781 

7 -.57681632E-1 

. -.388211786-2 


.164532506-1 


780./ 



1711 

DATA <0(1 ), 








1721 

81-541,576)/ 

-.780035256+1 


.139877306+1 


-.25321884 

A 

PCSU0706 

1731 

9 .322487346-* 

, 28.0 


.396958176+2 


- . 37004653E+1 

A 

PCSy0705 

1741 

1 .16677517 

, -.776469986-1 


28.6 


-.668498566+2 

A 

PCSU8710 

1751 

2 -.63689274 

, .286423516+1 


-.526615606-1 


28.0 

A 

PCSU0715 

176t 

3 .372579886+2 

, .589497796+1 


-.277375756+1 


.15813941,1489./ 

PCSW0720 

1771 

DATA (8(1 >, 1*577,612 >/ 







I78t 

4 



-.153783656-2 


-.992527236-1 

A 

PCSU0725 

1791 

5 .28554818E-1 

, -.943996296-3 


28.9 


.867509046-2 

A 

PCSU0730 

1881 

6 .68662619 

, -.11797733 


.559924066-2 


28.0 

A 

PCSU6735 

18H 

7 - . 16218554E-1 

, -.73896671 


.15437942 


-.730387776-2 

A 

PCSy0?46 

1821 

82686./ 







PCSU0745 

183i 

DATA <B(I), 

1-613,648)/ 


-.187538866+1 


.23580759 

A 

PCSy07S9 

184t 

1 -. 72532 363E-2 

, .858166986-4 


28.9 


.165196396+2 

A 

PCSy075S 

185t 

2 -.13176M66+1 

, .466236616-1 


-.581801646-3 


28.0 

A 

PCSU8766 

1861 

3 - • 159153656+2 

, .282567986+1 


-.765986676-1 


.112762576-2 

A 

PCSy0765 

1871 

42886./ 








1881 

DATA <»<I >.1*649,684)/ 







189i 

4 



-.864159876-4 


.640886276-2 

A 

PCSy0770 

198i 

5 -• 17674787E-2 

, .345587886-3 

• 

-.284473156-4 


.397718626-6 

A 

PC Sy 0775 

19H 

6 - .62376658E-2 

, -.244188866-1 

$ 

.168289316-1 


-.385344666-2 

A 

PCSy0789 

1921 

7 .22488754E-3 

, -.582537986-5 

A 

.821952526-1 


.17951874 

A 

PCSU678S 

1931 

8 -.42795977E-1 

, .141668826-1 

A 

-.144394386-2 


•39476388C-4 

A 

PCSU0790 

194t 

9 -.21585754 

, .23418848 

A 

-.18898115 


-.92288878E-2 

A 

PCSU0795 

1951 

1 .2321 1899E-2 

, -.787648866-4, 

1289./ 





1961 

DATA <»<I>, 1-685. 72D)/ 







1971 

2 -.92395825 

, .841246996-1 


-.175689136-2 

A 

.196940846-4 

A 

PCSU0895 

198t 

3 28.6 

, .138675846+2 


-.104168666+1 

A 

.224281946-1 

A 

PCSU9816 

199i 

4 - • 163566366-3 

, 21.8 


-.721416656+2 

A 

.527662836+1 

A 

PCSU981 

288t 

5 -.10193683 

, .381639866-3 


28.6 

A 

.185173226+3 

A 

PCSU082 

2811 

6 - .689656616+1 

, .932354396-1 


.298733496-3 


1488./ 


PCSU082 

2821 

DATA (0(1), 

1-721,756)/ 


-.169836226+1 

A 

.922421166-1 

A 

PC$y0836 

2831 

1 -.5668193SE-3 

, -.268624196-5 


28.8 

A 

.218211986+2 

A 

PCSy083 

2841 

2 -.76306746 

, -.742024366-2 


.288618866-3 

A 

28.8 

A 

PCSU684 

285i 

3 .244599796+3 

, -.157275726+2 


.42597236 

A 

-.416836946-2 

A 

PCSU084 

286i 

4 28.6 

, -.134067216+4 


.582272566+1 

A 

.123488186+1 

A 

PCSU0856 

287i 

5 -. 151326S6E-1 

, 21.0 


.775253646+4 

A 

-.154299096+3 

A 

PCSU685 

288t 

6 -.157723256+1 

, .168193686-1 


28.8 

A 

.423317666+3 

A 

PCSU086 

2#9i 

7 .563138796+2 

, -.146212896+2 


.32861514 

A 

28.6 / 


PCSU6865 

2181 

DATA <8(I) # 

1-757,792 )/ 


.36263253 

A 

-.133709426-1 

A 

PCSy#876 

2111 

1 - .25370804E-2 

. .135391886-3 


.186183526-4 

A 

-.454187446-6 

A 

PCSU8875 

2121 

2 - • 134883816+1 

# .654157116-1 


.638769866-2 

A 

.153242556-3 

A 

PCSU088 

2131 

3 -.341250796-4 

, .816572546-6 


.165358066+1 

A 

- .721645806-1 

A 

PCSU688 

2141 

4 -.695389676-2 

, -.173468486-3 


-.281157816-4 

A 

.119469296-5 

A 

PCSy689 

215t 

5 -.66826358 

, .164882136-1 


.782216566-2 

A 

-.524504156-3 

A 

PCSy089 

2161 

6 .674871656-4 

. -.261172426-5 


1286./ 





2171 

DATA (8(1), 1-793.828)/ 







2l8t 

7 -.368639246-2 

, .281474936-1 


.673418556-3 

A 

-.894661736-4 

A 

pcsy096 

2191 

8 .461451626-5 

, -.288529816-6 


.266611586-1 

A 

-.17289645 

A 

PCSU691 

228t 

9 .23968585E-2 

, -.386882666-4 


.228624736-4 

A 

.496698066-6 

A 

PCSU091 

22lt 

1 -.666999026-1 

, .35755448 


-.236110526-3 

A 

-.534498806-3 

A 

PCSU092 

222t 

2 -.209727656-4 

, -.120360296-5 


.446168286-1 

A 

-.244173 47 

A 

PC5W092 

2231 

3 .512455636-3 

, .641236996-3 


-.982651696-5 

A 

.128965986-5 

A 

PCSU893 

2241 

41280./ 








225t 

DATA ( B( 1 ) , 

1-829,864)/ 


.63215967 

A 

- . 70699932c- 1 

A 

PCSy0935 

226i 

5 .205898356-2 

, .243295826-4 


-.149667516-5 

A 

.969904476-0 

A 

PCSU0940 

227i 

6 -.10601885 

, -.321463656-1 


.688155436-2 

A 

-.123878186-3 

A 

PCSU0945 

228i 

7 .113726756-5 

# .362616296-7 


.124375136+1 

A 

-.13636679 

A 

PCSy#956 

229t 

8 .618783876-2 

$ -.459165966-3 


.433658376-5 

A 

-.602688536-7 

A 

PC5U0955 

2361 

9 -.311919636+1 

# .31189343 


-.546599836-2 

A 

-.228638826-3 

A 

PCSU0960 

23H 

1 .188617846-4 

, -.168796816-6 


1288 / 




PCSU096 


8-58 










II SUBROUTINE PMEXPNCPMANGL.HT.HL, AL,SL,PL,TL,UL«DL,ML, CANAL ) 

2* DIMENSION F<4,13>,G<4,t3>,A(4,13),ENTHC4.13) 

3 I REAL NL 

4i COMNOH/FLAC/IDEAL 

Si COMNON/ERRFLC/NERROR 

Cl IFINL.LT.l.) RETURN 

7l H-12 

St IFtPMANGL.GT. 1*3.2) GO TO 990 

9> IF(P«ANGL .LE. 86.) N-N-2 

101 IFCPNANGL .LE. 68.8) N-N-2 

lit IF (PNANCL .LE. SI. 6) N-N-2 

121 IF(PNANCL .LE. 34.4) N-N-2 

131 IF (PNANCL .LE. 17.2) N-N-2 

141 EN-N 

1SIC 

lSiCtMtt DIUIDE P-R ANCLE INTO INCREMENTS *»8* 

17IC 

181 PMANG -PfMNGL/57.3 

18i H-PMANG/EN 

28 i I-l 

2H J-l 

22t A( I , J )-AL 

23i F(I,J)-HT-HL 

241 10 G(I > J)-2.*F<I.J)*H/SORT(SO103.8F<I,J)/(A(I,J)M2>-1. ) 

25 ■ IF ( I .EQ.3 ) FC4,J)-F(1,J)»G(3,J) 

26t ir(I.E0.4) GO TO 30 

27i IFd.LT.3) F(I + 1.J)-F(1,J)*.5*C(I,J) 

28i I-IM 

291 20 ENTH( I,J)-HT-F(I,J) 

30t GO TO 40 

311 30 F(l,J+l)-F(l,J)+(G(l,J)»2.*G<2,J)»2.*GC3.J)+G<4,J))/6. 

32i I-l 

33i J-J+l 

34i GO TO 20 

35 1 40 CALL NOLIER(ENTH<I,J),P.3.T,Z.SL,RHO,GANA) 

36* IF( IDEAL. GT.0) GO TO 999 

371 A(I.J)-S0RT(GAMA«1716.48*Z*T) 

38i 90 IF(J.LE.N) CO TO 10 

39* HL-ENTH(I.J) 

401 AL-A(I.J) 

4H UL-SQRT<FU.J >*50103.) 

42 i PL-P 

431 DL-RHO 

44i TL-T 

45* NL-UL/AL 

46* GAMAL-GANA 

47* GO TO 999 

48* 990 URITE(6,930) PNANCL 

49* 930 FORMAT (S8H0 ***** PRANDTL-MEVER ANGLE CREATER THAN 103.2 PN ANGLPNEX 325 

50* 1E-.E12.6 ) 

51* NERROR-l 

52* 999 RETURN 

53* END 


PNEX 055 

PNEX 065 
PNEX 070 
PNEX 075 
PNEX 080 
PNEX 085 
PNEX 090 


M-fiU 


II SUBROUTINE PN I D ( PMANGL , HLOC AL # A LOCAL » PLOCAL « TLOC AL » ULOCA L # DLOC A L # 

8i lALOCAL) 

Sic 

4tc iititmtttmmiimttiiitttitttmmmtitiittttitittmttttiiittt 
ftc 

•IC THIS SUBROUTINE COMPUTES LOCAL CONDITIONS USING THE PRANDTL-AEVER EXPN 
7IC 

tic ittiittitiitttmimtmiimtttmitittittttttmmtitmtiiiittttt 


OiC 

REAL NUSTRA, MLOCAL 
I DIMENSION X(26) 


PAID O05 
PAID BIB 


..IC 

19ICSSSSS SET UP CONSTANT TERNS BEFORE ITERATATION ttttt 
I4IC 

.SI OANMA-1.4 

.Bl NUSTRN-NLOCAL 

l?l IF(NUSTRM.LE.l. > RETURN 

.81 PUSTRM-PLOCAL 

»0> TUSTRN-TLOCAL 

■Bt Cl-S0RT((GANMA4l. )/( GANNA-1. )) 

It C2-PMANGLI .0174533 

■81 C31 "SORT ( NU$TRMSt2- I . ) 

■31 C4*C1 SATAN (C31 /Cl ) 

_il C3-ATAWC31) 

S i CS-C3-C4-C2 

l C6-C1M2 

871 X(1>-NU$TRM 

SttC 

88tCSSStS PERFORM NEUTON-RAPHSON ERROR IS LESS THAN ONE PERCENT SSSSt 
3RtC 

311 DO 20 1-2,26 

38 1 F-C1*ATAH(X( I-l )/Cl )-ATAN(X( 1-1 ) )*C5 

331 XX-X(I-l)St2 

341 FPR IAE-C6/(C€4XX >-l . /( 1 •♦XX ) 

3St X( I )-X< I-l )-F/FPRINE 

381 IF(ABS( (X(I)-X(I-l) )/X( I-l ) ) .LE. .81) GO TO 4B 

37t 20 CONTINUE 

381 GO TO SB 

30 l 4B XXX-XCI) 

40IC 

41 tCSSSSS COMPUTE LOCAL CONDITIONS SSSSt 
48iC 

431 SB NLOCAL -SWT ( XXXtS2* 1. ) 

441 GRINDS- (GANNA-1. )/2. 

481 C- ( 1 . ♦GMINUSSMUSTRNf S2 ) / ( 1 . +GMINUStNL0CALtt2 ) 

48l PLOCAL-PU€TRMSCSS( GAMMA/ ( GAMMA- 1 • ) ) 

471 TLOCAL-TUSTRMSC 

48 I ALOCAL-49 . B2SSGRT ( TLOCAL ) 

40 i DLOC AL-PLOCAL/ (1716. 48STL0CAL ) 

8# l HL0CAL-B.24 STLOCAL 

81 I ULOCAL-NLOC ALIA LOCAL 

5£i RETURN 

63t END 


PAID BIS 

PAID B25 
PAID B35 
PAID B3B 

PAID B4S 


PAID 060 
PAID B6S 
PAID 070 

PAID 08B 
PAID B8S 
PAID 09B 
PAID IBS 
PAID 11B 


PAID 120 
PAID 125 
PAID 13B 
PAID 135 


PAID 1SB 
PAID 16B 
PAID 170 


ORIGINAL PMW IS 
OF POOR Q[ “ 


8-61 


•USROUTINC PR INTI 


4lC THIS SUlROUTINi PRINTS THICK-SKIN SOLUTION RESULTS. 

p 

71 COANQH/ THICK/ DTEO, TOOT, TEO, QC, OR, HCOUT, EHISS. MREC. 

•I t DT,CTIRE, PL. DTO, DTI. TE02 

•l C0NN0M/MANDV/INAX,DX<25).T(2S ),HPFLAC<25). NTHB.TSYNK.EHISIN.HCIU 

1*1 1,T*SF,TIUZ(S*>,HCIU2<S*>.TCASZ<S*>,T$INK2<S*>. TUI. TSINK, TOPS, HCIN. 

Ill 2 IEIRST.DCLT 

111 COMMON/ FLUFID/ philuo. tu, REST OS) 

13i DIMENSION TPRINT<25> 

141 DAT* STEVIE/ .47S*E-12/ 

1SI MM POMATUM*. 44X. 1RMOUTER UALL, 1*X. P2I.2/ 

1*1 1 (IN . 44X, 13, 17X, F21 .2 ) ) 

171 3*** FOMATUM , 44X, 1 SHIMMER UALL, 1*X. F21.2) 

111 4*** FORMAT! OH MTCUU)-E11.S.9H T CAS • F9.2.12H T SINK -F9.2.4X. 
1*1 17N9 <IU)«E11.S.9H TEO IU-F11.2.UH ST OUTER'EU .5./193X. 

Ml X9MST IWCR*E11.S./S2X. 

Ill 224MTENPEAATURE 9ISTRI9UTI0M ) 

Ml DO 5*4 I • 1. IRAK 

Ml TPRINT(I) - T(I) - 4S9.7 

Ml SS* CONTINUE 
Ml TUPNNT • TM - 459.7 

Mt OIN • - ENISIN X STEVIE X TUI *1 I ♦ ENISIN X STEVIE X TSINK XX 

I7t 14. HCIN X (TCAS - TUI) 

Ml TUIPPT • TUI - 459.7 

Mi TS • TSINK - 459.7 

3*1 TO • TCAS - 459.7 

3U TEQ8F*TE02-4S9.7 

Ml URITE(S,4***> HCIN.TC.TS.0IN.TE02F.DT0.DTI 

331 URITE(6.2**4) TUPRNT. (I, TPRINT'I), I • 1. INAX) 

341 URITE<6.3*M) TUI ART 

35 • RETURN 

Ml EMI 


H-62 


1* SUBROUTINE RADEGT ( ENC , HR , EMIS, TU, PE , TRE ) 

2*C CALCULATES RADIATION EQUILIBRIUM TEMPERATURE 

3 * COMMON/FLAG/ IDEAL 

41 SIGMA-4. 75892E-13 

5* TG-TU 

6 * 1 CALL M0LIER(HG,PE,2,TG,ZZ,SS,RR,GG ) 

7 * IF (IDEAL .EG. 0) GO TO 10 

8 < HG- .234S*TG+9.786E-6*TG*TG+943.o/TG-l .57 

9 J 10 DHDT-. 2345+1 .9S72E-S*TG-943.6/(TG*TG ) 

101 IF (TG .LT. 400.) DHDT«. 23366 

11* RESID-ENC*(HR-HG)-EMIS*SIGMA*TG**4 

12* IF (ABS(RESID/(ENC*(HR-HG ) ) ) .LT. .001) GO TO 30 

13* DERIU--ENC*DHDT-4.*EniS*SIGMA*TG*TG*TG 

14* RATE-RESID/DERIU 

15* TG-TG-RATE 

16* GO TO 1 

17* 30 TRE-TG 

18* RETURN 

19* END 


8-63 


iUftffOUTtNC SCTnu»(TSINKl,riUl,HC:ui.T6ASl.ftTIJZl,N0D£$l, 

1 Tl«Zl.TCAS2l.T»Cl.DX.Tli.H»fUO» 

XIS SUMOUTINC SETS L# THICK skin NMnrCM KND UNITES OJT If#uT 


f ^ 
2 SG 

«|w 

i 5 

ss 

t * 

•l'' 

V 5 C 

H W 

• lA. 
HO 


1 N 

»- 3 « 
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<X 

r 

(I 


• —* * 


MMA 
N «4L X 

3 ^OO 
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X -Hl- 


X tn« IAM« 
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Q 

r n ia — «t t 

o o ~ft x 

O XN*«X«00 

« — 'N^wf C 

•a ^ei-cnhh 
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II SUBROUTINE STABLKFIN, FOUT, GIN, GOUT, QI, TREC J 

aiCinilMMIItllllllllMIIIIIMIMIIIIIIIIIIMMMIMIIIIIIMIMIIIItlMI 
3iC • 

4lC THIS SUBROUTINE ESTABLISHES THC STABLE TINE INCREMENT FOR • 

SIC THC THICK SKIN OPTIONi THC LESSER OF THE TIMES RCOUIRED FORI 

«iC THE IHHER AMD OUTER UALLS TO REACH THEIR RESPECTIVE I 

7iC RADIATION EQUILIBRIUM TEMPERATURES FOR THE GIVEN NET HEAT I 

QIC FLUX AMD HEAT CAPACITY. I 

SIC t 

lUC IF THE STABLE TIME INCREMENT IS LESS THAN 1 SECOND BECAUSE I 

UK THE HALL TEMPERATURE AND UALL EQUILIBRIUM TERMPERATURE ARE I 

l2lC U I THIN 1 DEGREE. THEN THE STABLE TIME INCREMENT IS SET TO 1 I 

1 3*C SECOND. B 

14 iC I 

lS)C888l88lfltm888lf888l88888l8IS8ff8ll888f888llf 188188188188818811818888 
IS • COMMON/ THICK/ DTEO, TDOT , TEQ, OC. OR. HCO, EMISS, HREC, 

171 1 DT,CTIME. PL. DTO. DTI, TEQ2 

It I C0NN0N/HAMDV/IHAX,DX(2$).T(2S>,NPFLAC<2S>, NTIUZ,TSYNK,EMISIN,HCIU 

19 i 1 , TGSF . TIUZC Si I ,HCIUZ (S# >, TGASZ ( SO ) , TSINKZ( SO ), TUI . TSINK , TGAS,HCIN, 

28* 2 IF IRST,DCLT 

21 l COMMOM/ERRFLG/NCRROR 

221 DATA STEVIE/ .4766E 12/ 

23i 3666 FORMATUHO, 13X, 3€HI IINEUTON-RAPHSON FAILED TO CONVERGE, 

241 1 S3H ON A VALUE F OR THE INNER EQUILIBRIUM TEMPERATURE AT, 

25 1 2 FIB. 3, 1 IN SECONDS!!! > 

2S iciiiiiiiiii iiiiiiiiiiii iiiiiii i linn i mu iiiiiiiiiiii iiiiiiiiimiiiii t 

27 1 C I 

28 iC STABLE TIME INCREMENT I OUTER UALL I 

29 * C | 

3# t ciiiiiii inn if iiiiiif if iiimiif iiiiiiiiiiiiiii i min mu i mi in in t 

3H CALL RADEQTCHCO, HREC. EMISS, T(1 ), PL, TEQ) 

321 00 • QC - OR 

331 DTO • ABSCGOUT 1 (TEO - T<1))/<Q0 - FOUT IGOUTi C Tc 1 > - T(2)))) 

341 IF < DTO. LT . 1 . .AND. C ABS(TEQ-T( 1 ) ) ).LT. 1 . ) DTO -1.8 

35 1 ciiiii ii ii iii iiiiiiiiiiiiiii mil mil iiiiiiii iiiiiii in mi in miiti if 

36 1 C I 

37 l C STABLE TIME INCREMENTI INNER UALL I 

38 1 C | 

39«C FOR ADIABATIC INNER UALL, SET RADIATION UIEU FACTOR AND I 

48lC HEAT TRANSFER COEFFICIENT EQUAL TO ZERO. I 

41IC | 

42 tCIIIIIIIIIIIIIIIII II IIIIIIIIIIII || llllllllllllllllllllllllllll llimi || | 
431 IF (TMISIN .NE. i. .OR. HCIN .NF. 8.) GO TO 10 

441 DTI • 1 . 8E+66 

45* TE02 • 8. 

461 CO TO 28 

47 1 18 CONTINUE 

48 1 TEQ2 • T ( I MAX ) 

491 AEQ2 - EMISIN I STEVIE 

581 BE 02 • HCIN 

51 i C102 • AE02 I TSINKII4 ♦ MCIht TGAS 

52IC88II8IIIIII8IIIII88IIIIIIIIIII888I8II8I8IIIIIIIIIIIIIIIIIIIIIIIIIIIIHI 

53 * C I 

54 iC FOR INNER UALL, IF SUM OF RADIATIVE AND CONVECTIVE HEAT I 

55 iC FLUX IS ZERO. SET RADIATION EQUILIBRIUM TEMPERATURE TO I 

56 i C -459.7 F. | 

57 1 C I 

§9108888811 IBIIIf 1811 I llllf I I I I I If If If If Mil If I If I Ilf If If tit I I I IIIIIII If III 
fit IF (CE02 .CO. 8. ) TE02 • 8. 

68lCI68l8limill888ll888l8llll88lllll8ll8lllllll8llllllll8l8lllll8lll888ll 

61 >C l 

6§lC THE NEUTON-RAPHSON METHOD IS USED to SOLVE FOR THE I 

63 1 C RADIATION EQUILIBRIUM TEMPERATURE. I 

84 1 C I 

66 i C888IIIII IIIIIIIIIIIIIII I IIIIIII till I III! II II I I IIIIIIIIIIII I III III Hill I 

66 1 IF ( CE02 .NE. 8.) CALL NCUT<A£02. BE 02 . CE02, TE02, ICANT ) 

671 IF (ICANT .NE. 6) GO TO 918 

68 1 DTI • ABS ( GIN! ( T( I MAX )-TEQ2 )/(FlNIGlNS(T( I MAX - 1 )-T< 1NAX ) ) ♦ 0 1 > ) 

69 i IF(DTI.LT.1..AND.(ABS(TE02-T(1MAX) n.LT.l. ) DTI *1.8 

761 28 CONTINUE 

711 DTEQ • AMINKDTO.DTI) 

72 • RETURN 

73i 918 UR I TE ( 6 , 3668 ) CTIME 

741 NERRORM 

76 i RETURN 

76i END 
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II SUBROUTINE SUCVL2<RN,PHI,ENCL.ENCT> 

2* COARON/FLUFLD'TINF,PINF,RHOINF.AINF.RHOE,PE.TE.ME.XRUE,UE.RHOO,PO, 

31 1TO,HO,XHUO,RHOU,TU,HU,XRUU»RHOS,TS»HS,XRUS,RMOR,TR,HR,XRUR,V,XRE 

41 2,REL,HIHF,XRUINF,PR.PU,DU.TU.XRUU.XHACHU.U , J.HU. alpha. 

S 1 3RH0STL , HSTL. TSTL , XHUSTL . RHOSTT, HSTT , TSTT, XMUSTT . HRECL . HRECT . S 

Cl *.C INF, CAHAU, CARAS, GARAE.GAHAO.GU.GARAC.PRL.PRT 

7l CORRON/CVL/ HSL.PSL.RHOSL.XRUSL.TSL.WSL 

•I CORRON/STAC/PSTAG. RHOSTC, TSTG. XRUSTC 

9 1 CORRON/ERRF LC/NERROR 

19i CORRON/FLAC/ IDEAL 

111 RED* DUIUUS2 • 4RN/XRUU 

18i REL-RED 

131 1 PSU-PO 

141 XSU*XHUO 

1SI RSU*RHOO 

1C< TSU-TO 

171 RUSU'RHOU 

lit XRUUSU*XHUU 

191 PO-PSTAC 

29l RHOO-RHOSTG 

211 XRUO*XRUSTG 

22t TO* TSTG 

231 CALL ROLIER(HU,PO,9,TX,ZX,SX,RHOU,GX ) 

241 CALL HANSEH(XRUU,PO,TU) 

2S< CALL FAYRID(RN.ENF) 

261 PRL* .71 

271 PRT*.71 

28> PO-PSU 

29 i RHOO*RSU 

39i XRUO*XSU 

3H TO-TSU 

32i RHOU-RUSU 

33i XRUU-XRUUSU 

341 ENCL-9.797*ENF*C0S(PHI*.«174 >«*1 . 1 

3S« 19 DUDX*SORT(2.t(PSL-PU)/RHOSL)/RN 

3£t RRR*RHOSTTtXRUSTT/ (RHOOIXRUO ) 

37i ENCT • 1 . 94SRHOOK . 8IXRU0M . 2*USLM . 6*RRRM . 81DUDXI* . 2/PRTt I . 667 

38t S-PRLM.667 

38 i 99 RETURN 
491 END 
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70lC I TRAN IS A FLAG TO INDICATE POINT OF LAHINAR/TRAHSITION OR 
71 IC TRANSITION /TURRULENT FLOW FOR PRINTOUT AMD PLOTTING. 
72 «C IT IS RESET TO ZERO WITH EACH PASS THROUGH THIS SUDR 
73iC ITRAN-0 NO CHANGE IN TYPE OF FLOW 
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tl SUtROUTINE TINTS (A,AT-l,IT,C»CT > O.OT,E # ET,r,FT,6,8T > IL.MX ) 

itC LINEAR INTERPOLATION FOR S VARIABLES - IF IL'I.TX INDEPENDENT 
3tC UARIAILE IS DECREASING 


4t 


DIMENSION ATC 10 ),8TC 10 ),CTC 10 ),DTC 10),ETC 10),FTC 10),QTC 10) 

Si 


IF(XL) 1,1,10 

Si 

1 

IFCA-ATCl) >2,3,4 

7i 

2 

I-.222E+30 

li 


RETURN 

Oi 

3 

•-•TCI) 

III 


C-CTCI) 

1H 


D-DTC1 ) 

lit 


E-ETC 1) 

13i 


F-FTC1) 

141 


G-QTCl) 

15* 


RETURN 

1SI 

10 

IFCATCl >-*>2,3,11 

171 

4 

DO 9 1-2. NAX 

lit 


IFC A-ATC I > >20, 18,9 

10* 

0 

CONTINUE 

2#t 


CO TO 2 

2H 

18 

•■•T<1> 

22* 


C-CTCI) 

231 


D-DTCI) 

24t 


E-ETC 1) 

25* 


F-FTCl) 

25* 


C-CTCI) 

27* 


RETURN 

28> 

20 

RATIO- C A-ATC 1-1 > )/C ATC I >-ATC 1-1 > > 

29* 


• -•TCI-1) ♦ RATIOtCiTC I )-IT C 1-1 ) ) 

30t 


C-CTCI-l) ♦ RATIOS ( CT ( I >-CT ( I- 1 ) ) 

31i 


D-DTCI-1) ♦ RATIOSCDTC I >-DT C 1-1 ) > 

32* 


E-ETC 1-1 > ♦ RATIOS! ETC I >-ETC 1-1 ) ) 

33t 


F-FTCI-1) ♦ RATIOSCFTC I )-FTC I-I > > 

34i 


C-CTCI-l) ♦ RATIOSCCTC I >-CTC 1-1 ) ) 

35t 


RETURN 

36i 

11 

DO 16 1-2, NAX 

37t 


IFCATC I >-*>20,18,16 

31* 

IS 

CONTINUE 

39* 


CO TO 2 

401 


END 
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II SUBROUTINE URUNL<TRFLAC.ELTRAH,ELL.ELT.ELTP.PARA1,PARA2,ENL,EL.UX, 

81 1 NCFLG ) 

3i DINENSION ENEX(7>,RETRX(?),XX(7) 

41 CONNON/FLUFLD/TINF ,PINF,RH01NF ,ainf , RHOE, PE, TE, HE.XHUE , UE. RHOO.PO, 

fit 1TO,HO,XNUO.RHOU.TU,HU.XHUU,RHOS.TS.HS,XNUS,RHOR.TR HR.XNUR.U.XNE 

Si 2, REL, MINF.XNUINF. PR, PU. DU, TU, XNUU.XNRCHU.UU.MU, RLPMA. 

71 3RHOSTL.HSTL.TSTL.XHUSTL.RHOSTT,HSTT.TSTT.XNUSTT,HRECL.HRECT.S 

• I F.GINF, CANAU, CAN AS . OANAE , GAN AO , GU, GAN AG . PRL.PRT 

9l DATA (EREX(K>,K*l,7)/9.5. 1.0, 1.5,2. 0. 3.5,4. S.S.S / 

l#i DATA (RETRX(IC>.r*1.7>/ S . 30193, S. 54497. 5. 81291 .6. 0,6. 0,6. #7918. 

Ill 1 6,25527/ 

12l DATA NX2/7/ 

13i NTRANS*TRFLAC*.901 

Hi IF(NTRANS.E0,6) CO TO 6* 

15i IF(NTRAN$.CT.6> GO TO 70 

161 PARA-PARA1 

17l IF ( PARA2. LT .PARAl ) PARA* PARAS 

111 IF(NTRANS-4 >39,49,59 

191 39 ELTRAN*PARAFXNUE/(RHOEFUE ) 

29t CO TO 190 

211 49 ELTRAN*<PARAFXNUE/.664)FF2/<RH0STLFXNUSTLFUE> 

22 ■ CO TO 199 

231 59 ELTRAN*2 . 26FPARAFPARAFXNEFXNE FRHOE <! , 4FXNUEFF 1 . 6/ ( RHOSTLFXNUSTLF 

241 1UEFF.61FENL 

25i IF(NCFLC.LE,9) CO TO 109 

26t IF(NCFLC. LT . 3 ) GO TO 51 

27i ELTRAN-ELTRANFEL/ELL 

29 i GO TO 199 

291 51 IF<ELTRAN.GT.(UX/2. >> CO TO 52 

391 ELTRAN— ,5<UX*AL0C( 1.-2. FFLTRAN/UX ) 

311 GO TO 199 

321 52 ELTRAN* 1090. FELTRAN 

33i CO TO 199 

341 69 CALL EDPARN(ALPHA.PARA) 

36 i GO TO 59 

Ml 79 CALL TBLIN(XHE,ENEX,RELG,RETRX,X,XX,NX2) 

37 1 ELTRAN* ( 19. FFRELC )FXNUE/ ( RHOEFUE ) 

3Bi 199 ELTP-ELT 

39 i IF (EL .LT. ELTRAN) CO TO 999 

49i IF (NCFLG. LE. 9 ) GO TO 119 

41 1 ELT • ( EL-ELTRAN JFELT/EL 

421 GO TO 999 

43l 119 ELT-ELTP-ELTRAH 

441 999 RETURN 

45 i END 
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660 


II SUiROUTINE PL0TL0<TINE,NHFLAG,ARID2F,ATRE,TZ,ZZ,UZ,ALFA9T, 

21 I DELTAT, ITHICK, LNGPLT, ITINIT,MAXTME,DEUICE, IMCOPV ) 

3«C 

4iC DECLARATIVE STATEMENTS 

SIC 

Si INTEGER 9DVPNT , TRANNE , DEU I CE 

7i PARAMETER Jl-17, J2-S98, J3-32 

ft DIMENSION TIMEH J2 ),XNE1 ( J2),REL1 ( J2),0C1 ( J2 ),0C1T< J2 ), 

9l ITREF1 ( J2 ),TU1 ( J2 ),PE1 < J2 ),CPl( J8 ),9ETA1 ( J2 ),ALPHA1 ( J2 ), AKLZ1 ( J2 )# 

19i I HFACZ1 ( J2 ) . OT < 2 ) , TT ( 2 ) , TU2 ( J2 ) 

111 DIMENSION U2(S9 ),TZ(69 ), ALFA9T(59 ),ZZ< 59), DELTAT (59 ) 

121 DIMENSION 9DVPMT<4 ), TRANNE (6), I»UF(2#9),LAIELC29,6), IPAKC49) 

13»C 

14IC 

1SIC 

161 COMMON/LAIS/LAiEl 

1? I COMMON/DI ftU JO/IDVPNT, TRANNE , ICASE 

18 * COMMON/ AXDAT/XOR, XSTEP ,XAXIS 

19t COMMON/ ARRAV/QC1,T I NE1 , BETA 1 , ALPHA 1 ,AKLZ1 , 

2Si BHFACZ1 ,XME1 , REL1 ,PE1,CP1 , TU1,TREF1,0T,TT,TU2,9C1T 

21 1 COMMON/MAX/OH, BETAM, AKLZN, HFACZM, PEN, CPPM, TREFM 

221 CONNON/PIOTT/JRCD, IT,NT,NT2,NFF,PCT 

23i COMMON IPT , IBUF, IPAK 

24 1 C 

26tC 

2StC SET X-AXIS LABEL UNITS AND INCREANENTS 

27tC 

29t OMAX-MAXTME 

29 tC DELTX-NAXTME/19 

3#iC 

31 tC INITIALIZE PLOTTING DEVICE 

32IC 

331 DEVICE-1 

341 GO TO (31 9, 329 ,338), DEUICE 

jc* c 

36iC INITIALIZE TEKTRONIX CRT 

37iC 

32i 319 CONTINUE 
39 i CALL ERASE 

49>C 

41 IC START PLOTTING 

42 * C 

43 i CALL TEKEGNC489) 

441 CALL BASALFC'L/CSTD ) 

46 i CALL MIXALF ( 'STANDARD ' ) 

48 1 GO TO 349 

47 * C 

48 *C INITIALIZE MOPS TERMINAL 

49tC 

691 329 CONTINUE 

sue 

62 tC INSERT CODING HERE 

63tC 

541 GO TO 349 

55 1 C 

56 *C INITIALIZE NICROFILN HERE 


STIC 
691 339 
69 1 C 
69 1 C 

sue 

62i 349 
63 1 C 
64 1 C 
66t 
66 1 C 
67 1 C 
68iC 
69i 5 

79 1 C 
71 IC 
72*C 
73t 


CONTINUE 

CONTINUE 

PLOT REFERENCE PLOTS 
IF ( ICASE-2 ) 5, 19, 29 
REFERENCE SPHERE PLOTS 
CONTINUE 
TIME US ALTITUDE 

CALL AXSPLT(9.9,0NAX,8.9,X0R,X$TEP,XAXIS> 
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3# 


74*C 

75i CALL HNDVLOC INCOPY, TZ,ZZ, NT, 1,1 ) 

75 1 C 

77IC TIRE US VELOCITY 

7SiC 

79* CALL NNDYLOC IMCOPV,TZ,UZ,NT,2, 1 ) 

MIC 

IHC TIRE US ANGLE OF ATTACK 

M*C 

93* IF (ALFAOT . LE. 0) GO TO 10 

me 

•SI CALL NNDYLOC IHCOPY, TZ,ALFA9T, NT, 3, l > 

MIC 

•7i CO TO 49 

MIC 

MIC TIRE US ODOT REF 

MiC 

911 19 CONTINUE 

92l CALL DSMRL0<IT,QC1T,VM1N,VRAX) 

931 CALL AXSPLT(9.9,YAAX,6.9,V0R,VSTEP,VAXIS> 

941 CALL PPL0T1 ( 11 , ICASE ) 

96iC CALL 9LNK1(9.9,2.S,5.9,6.2,3> 

Mi CALL DRAULO ( XOR , XSTEP , XAX IS, YOR , YSTEP , VAX IS#TIRE1,QC1,IT,3,3 

97t 1 IHCOPY ) 

Ml CALL DRAULOCXOR,XSTEP,XAXIS,YOR,VSTEP,VAXIS,TIRE1,OC1T,IT,2, 


99 


99i 

lMi 
191* 

1921 
193i 
194t 
196i 
IMS 
197t 
lMt 
1991 
1191 

line 

U2tC 

113*C 

114IC 

1151 29 

119*C 

117*C 

line 

119IC 

1M*C 

121 iC 

1M*C 

183iC 

184 tC 

185IC 

lMi 

187iC 

1M* 

129IC 

139IC 

13HC 

132i 

133IC 

134 iC 

13Sf 26 

139*C 

137*C 

139IC 

139i 

149IC 

141 IC 

1481 

143IC 

144IC 

145 *C 

1451 

147 *C 


1 IHCOPY) 

CALL LINES! 9 IT )U (R )AD (E)0 (T >EMP.f' , IPAIC, 1 ) 
CALL LINES(MT)U ( D )EC. < F >• ' , IPAK , 2 > 

IPT — l 

CALL LECENDUPAK,2,2.S,6.2> 

CALL INTN0CITINIT,3.4,6.2> 

IF (II COPY .LE. 9) GO TO 99 

CALL MDCOPY 

CALL ERASE 

CONTINUE 

CALL ENDPL(9 ) 

GO TO 49 


CONTINUE 

IFCLNGPLT .EO. 9) GO TO 39 

LONG PLOT 
LAMINAR FACTORS 

LAMINAR FACTORS 

CALL HNDVLOC IMC0PV,TIRE1 , AKLZ1 ,IT,4,3) 
IF( ARIDEF .ME . 1 ) GO TO 26 
DEFLECTION FACTORS 

CALL HNDVLOC IHCOPY, TINE1,HFACZ1, IT, 5. 3 > 
CONTINUF 

LOCAL MACH NUMIER 

CALL HNDVLOC IHCOPY, TIRE1,XME1, IT, 6, 3) 

IFCNHFLAG .EO. 1) GO TO 27 
REYNOLDS NUMiER LOCAL/ 1.9E6 
CALL *HNDYLO<IMCOPY,TIRE1,REL1,IT,7,3) 
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27 


28 


I4HC 
1481 
189*0 
181 1 
168*0 
183IC 
164tC 
188IC 
18ft 
167IC 
1881 
ISOiC 
169*0 
161 iC 
1621 
163*0 
164IC 
166IC 
1669C 
1671 
1*>3*C 
169*0 
17#iC 
17H 36 

178*0 
173*0 
174* 

178*0 
178*0 
177*0 
178*0 
179*0 
188* 216 
181*0 
188*0 
183*0 
184* 

186*0 

186*0 

187*0 

188* 

189*0 

196* 

191*0 

192*0 

193*0 

194* 

196*0 
198* 236 

197*0 
198* 

199* 

298*0 
881* 37 
282*0 
283*0 
284*0 
286* 

286*0 
287* 49 

288* 

289* 


GO TO 28 
CONTINUE 

REFAC 

CALL HNDVLO < I HCOPY , T I HE 1 , HF AC Z 1 , I T , 8 , 3 ) 
CONTINUE 

LOCAL PRESSURE LIF/FT-SQ 
CALL HNDYL0(IHC0PY,TIREl,PEl,IT,9,3> 

C SUI P 

CALL HNDVLO* IHCOPV - TIRE1,CP1,IT,16,3) 
SHORT PLOTS 
CONTINUE 

1F( ARIDEF .NE. 1 ) GO TO 216 
CONTROL SURFACE DEFLECTION PLOT 
CALL HNDVLO < IHCOPV,TINEl,DELTAT, IT, 13, 3 ) 
CONTINUE 

Q DOT CITU/FT-SO-SEC) PLOT 
CALL HNDVLO* IKC0PV.TINE1.0C1, IT, 11,3) 
TITLE TUAL (DEG F) PLOT 
CALL HNDVLO ( IHCOPV, TINE1 , TREF1 , IT, 12,3) 
IF(ATRE .GT. 6.881) CO TO 236 
T UALL (DEG F) PLOT 
CALL HNDVLO* IHCOPV, T INEl , TUI , IT, 12,3 ) 
CONTINUE 

IFC ITHICK .NE. 9) CO TO 37 
GO TO 46 

CONTINUE 

TUALL (DEG F) PLOT TU2 

CALL HNDVLO* IHCOPV. TINE1.TU2, IT, 12, 3) 

CONTINUE 

RETURN 

END 
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SUBROUTINE DSnHLO<NP, A, YORI, YLAST ) 
DIMENSION A (S00) 

VORI-A( 1 ) 
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10 . 


OFF-LINE PLOTTING 


Off-line plotting is basically the same as explained in section 4, except 
that the data is ready o be plotted and stored in a secured file. The 
name of the file contain ng the runstream is different; it is OFFLPT 
(off-line plotting) and resides in secured file ES35-N06516*L0REN. 

The runstream listing for OFFLPT is shown in figure 2; the program listing 
is shown in figure 3. 


ItiFREE TPFi. 

2»«ASC,A ES32-L78771 tOLHEDO 
3lfUSE OL. , ES32-L78771 tOLHEDO 
4»#ASG,T 9.,F/1/TRK/5O0 
5 1 0ASG, T TPFf.,F/l/TRK/800 
6>«ERS TPFf . 

7 1 f COPY ES35-N06216tLOREN.,TPFf 
8 » tED, I TPFf. HAP 
9HIB TPFf 
10HIB DISSPLAtl RV 
1USEG HAIN 

12 J NOT DISSPLAtTRV . QOTKEG 

13 1 IN LOREN. QOTKEG 

14« IN DEHONA,PLOTDE,HENUPT 

15 < END 

16 * f PREP 

17tfHAP TPFf . HAP, TPFf . ABS 
18 1 fXQT ABS 

19 1 f ADD OL.RTLS-EX/32779 
20 i f ADD OL . C 1234/DATA 
21 t tADD 0L.C12345/DATA 
22 > f ADD 0L.BD3RD1/DATA 
23*fADD OL. DIRECTORY 


Figure 2. - Runstream listing OFFLPT. 
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11. SUBROUTINES USED BY MINIVER 


Routine name 

Symbolics available 

Source system 
or library 

H800 

Yes 

MINIVER 

AXSPLT 

No 

DISSPLA 

AIR62 

Yes 

MINIVER 

ATMS4 

Yes 

MINIVER 

BASAL F 

No 

DISSPLA 

NINTRP 

Yes 

MINIVER 

CHEEVY 

Yes 

MINIVER 

CRSFLW 

Yes 

MINIVER 

CURVE 

No 

DISSPLA 

DETRAL 

Yes 

MINIVER 

DINT 

Yes 

MINIVER 

DINT1 

Yes 

MINIVER 

DOWN ID 

Yes 

MINIVER 

DRAWLO 

Yes 

MINIVER 

DRIVEL 

Yes 

MINIVER 

DSMMLO 

Yes 

MINIVER 

EDPARM 

Yes 

MINIVER 

ENDPL 

No 

DISSPLA 

ERASE 

No 

PLOTIO.LOCALIB 

ERTRAN 

No 

LOCAL IB 

FINALT 

Yes 

MINIVER 

FAYRID 

Yes 

MINIVER 

FLOW 

Yes 

MINIVER 

GAUSS 

Yes 

MINIVER 

GRACE 

No 

DISSPLA 

GRAPH 

No 

DISSPLA 

GRID 

No 

DISSPLA 

HANSEN 

Yes 

MINIVER 

HDCOPY 

No 

PLOT1 0 .LOCAL I B 

HNDYLO 

Yes 

MINIVER 
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Routine name 

Symbolics available 

Source system 
or library 

HEIGHT 

No 

DISSPLA 

INTNO 

No 

DISSPLA 

IOWAIT 

No 

PLOTIO, LOCAL IB 

LEGEND 

No 

DISSPLA 

LINES 

No 

DISSPLA 

MATRES 

Yes 

MINIVER 

MESSAG 

No 

DISSPLA 

MIXALF 

No 

DISSPLA 

MOLIER 

Yes 

MINIVER 

NEWOUT 

Yes 

MINIVER 

NEWT 

Yes 

MINIVER 

NOCHEK 

No 

DISSPLA 

OPTMYZ 

Yes 

MINIVER 

PCSW 

Yes 

MINIVER 

PLOTIO 

Yes 

MINIVER 

PMEXPN 

Yes 

MINIVER 

PMIO 

Yes 

MINIVER 

PPLOT1 

Yes 

MINIVER 

PRINT1 

Yes 

MINIVER 

RADEQ 

Yes 

MINIVER 

SETMUP 

Yes 

MINIVER 

STABLE 

Yes 

MINIVER 

STOCK 

Yes 

MINIVER 

STORED 

Yes 

MINIVER 

SWCYL2 

Yes 

MINIVER 

TBLIN 

Yes 

MINIVER 

TEKEGM 

No 

DISSPLA 

TINT6 

Yes 

MINIVER 

TRANS 

Yes 

MINIVER 

VRUNL 

Yes 

MINIVER 

WRINP 

Yes 

MINIVER 

XINTAX 

No 

DISSPLA 

YAXANG 

No 

DISSPLA 

YINTAX 

No 

DISSPLA 
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